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Preface 



The main reason to study B meson decays is their sensitivity to the flavor structure of 
nature. Indeed, the fact that the b quark is so heavy makes B physics a rich source of very 
different processes, and leads to a very rich phenomenology. An additional consequence 
of the large mass of the b quark, that makes B decays interesting on the theoretical side 
is the fact that the always troublesome strong interaction effects can be handled within 
a heavy-quark expansion, allowing for theoretical predictions of acceptable accuracy. In 
fact, one of the differences between B- and D-physics is that decays of D mesons are 
theoretically much "dirtier". B decays have therefore triggered intensive research on the 
QCD side, that has witnessed a huge progress in the last decade. 

One of the main points in the B physics program is the search for physics beyond 
the Standard Model (SM). To that end, a number of dedicated facilities have been taking 
data for many years, achieving a long list of new discoveries. Specifically, the B-factories 
Babar and Belle, and the hadronic machine at Tevatron, with its experiments CDF and 
DO, have made this possible. Joint work by theorists and experimentalists has led to the 
appearance of several puzzles in B decays which could be interpreted to be due to New 
Physics (NP). However, these may as well disappear as more accurate measurements are 
made, or new understanding is accomplished on the theoretical side. The imminent start 
up of the Large Hadron Collider (LHC), with its B physics experiment LHCb, and the 
possibility of a super-B factory will certainly play a very exciting role in the search for 
physics beyond the SM, and once NP is found, in understanding its nature. 

The motivations that guide us in the search for NP -that so many times have triggered 
false alarms- are many. They can be classified in three major classes. Motivations belong- 
ing to the first class are related to observational facts that our theory (the SM) cannot 
reproduce. The baryon asymmetry in the universe, as understood by now, requires the 
three so-called Sakharov conditions. One of such conditions is an amount of CP violation 
that the SM cannot account for. The solution most probably requires new sources of 
CP violation. The presence of dark matter and dark energy is also an observational fact 
that has not found its solution within the standard theory, and which points towards the 
existence of new particle content. 

Motivations belonging to the second class are related to observations that can in prin- 



2 



Preface 



ciple fit into the theory but which would then seem extremely unnatural. The SM predicts 
two independent sources of CP violation in strong interactions. One is driven by the 
presence of instantons and the other comes from the diagonalization of the quark mass 
matrix. The measurement of the electric dipole moment of the neutron tells us that these 
two -in principle unrelated- contributions must cancel to one part in 10 9 . The classical 
solution to this strong CP-problem is to postulate a U(1)a (Peccei-Quinn) symmetry that 
dynamically sets this small number to zero. This solution requires the existence of (at 
least) a new particle, the axion. The hypothesis of inflation, introduced originally as a 
solution to the horizon and flatness problems, is pretty much accepted by the physics 
community nowadays. However, whatever drives inflation is still unknown, and again the 
many proposed mechanisms involve physics beyond the SM. The most striking fine-tuning 
problem is maybe the one related to the cosmological constant, which states that several 
unrelated contributions to the vacuum energy and the bare cosmological constant must 
cancel to one part in 10 120 in order to agree with observations, which is preposterous. The 
Higgs fine-tunning problem arises from the instability of the mass of a scalar particle to 
radiative corrections. The difference between the electroweak (E W) and the planck scales 
requires a fine-tunning of one part in 10 30 . Supersymmetry solves this problem protecting 
the mass of the scalars, which are not protected by gauge invariance, but the solution 
might as well be ultimately of different nature. In any case it should manifest itself as NP 
at ~ ITeV. 

Motivations belonging to the third class are not related to problems, but to unanswered 
issues. The inclusion of gravity in the standard picture is an issue that has led to the 
study of extra dimensions and string theory. A thorough theoretical investigation over 
decades indicates that quantum gravity can only be merged with the SM together with a 
significant amount of NP. The SM suffers by itself from its own unanswered issues, more 
related to this thesis than the problem of gravity. The SM contains 28 free independent 
parameters, a feature that goes directly against the quest for unification that motivated its 
foundations. An attempt to reduce the number of gauge parameters, through the gauge 
coupling unification, requires Supersymmetry and a long list of new particles, introducing 
even more parameters. From the 28 free parameters in the SM, 22 are directly related to 
the flavor sector. The hierarchy of the CKM matrix and the hierarchy in the masses of 
the quarks and leptons of different families should have an explanation beyond the SM, 
and constitutes part of the SM flavor problem. The existence of 3 generations of quarks 
and leptons is intriguing, and the existence of a fourth family is neither theoretically nor 
experimentally excluded. The fact that neutrinos are not massless provides new puzzles 
concerning the mixings, hierarchies and flavor violations in the leptonic sector to which 
any NP would have a definite impact. 

So it is clear that the search for new physics is strongly motivated. Past experience 
has also taught us that purely theoretical arguments are a powerful tool that leads to 
actual discovery. The positron was predicted by Dirac as a product of the unification of 
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Figure 1: The progress in the determination of the Unitarity Triangle. 



relativity and quantum mechanics. The charm quark was postulated in order to provide a 
GIM mechanism to suppress FCNC's. The third family of quarks was postulated in order 
to allow CP violation through the KM mechanism. The weak gauge bosons W, Z were 
predicted with the right masses as a realization of the GWS theory of weak interactions. 
These and other successes are the proof that the theoretical method is in the right track, 
and that research on the problems mentioned above will bring discoveries as impressive as 
old problems did in the past. However, as important as theoretical elucidation might be, 
experiment is the only way we get to know the world. Both must be connected through 
the ever crucial link: phenomenology. 

The role of flavor physics is central to this enterprise at the phenomenological level. 
The search for new physics begins with the understanding of the SM itself, and the precise 
determination of its parameters. The majority of these parameters are flavor parameters 
that are not very well determined in the present, but where impressive progress has been 
made (see Fig{T|). In fact the determination of the flavor parameters of the SM will soon 
become of the •precision type. Flavor parameters can only be determined through flavor 
physics, and the determination of the flavor parameters is necessary in order to infer not 
only quantitative but also qualitative aspects of the physics above the EW scale, and of 
the mechanism responsible for the flavor structure that we see at low energies. 

The importance of flavor physics cannot be overestimated. If there is really new physics 
at the TeV scale, it will be most probably detected at the imminent LHC at CERN within 
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2 to 5 years from now. However, revealing the true nature of those new particles is a 
much more uncertain task. The relevance of the discovery of NP is certainly to find the 
answers to the problems raised above. Will the discovered new particles provide a clue 
about the origin of flavor? Or about the mechanism of electroweak symmetry breaking 
(EWSB)? Or a sensible reason for the existence of large hierarchies? As far as we know all 
these problems may be related, or they may be as well completely disconnected. Will it 
be clear from the beginning that we are discovering Supersymmetry, or extra dimensions, 
or something else? Indeed, the direct search for new particles is an important but not 
the only part of the New Physics program. And this is where flavor physics comes at 
hand. Flavor physics constitutes a powerful arena on which to investigate detailed aspects 
and properties of the new physics and its possible role on the observed phenomenology. 
Indeed, flavor physics has already provided important constraints on the properties of 
the new physics concerning its CP and flavor violating structure. It has even raised 
its own problems. The suppression of flavor-changing neutral currents (FCNC's) is an 
observational fact that introduces a strong requirement in the flavor nature of the new 
physics. Without any mechanism that suppresses FCNC processes, any generic scenarios 
for NP should be suppressed by a scale larger than ~ 10 3 TeV. Therefore, any TeV NP 
must have a very specific flavor structure in order to satisfy the flavor bounds. This is 
the NP flavor problem. The non standard CP violating phases that appear, for example, 
in supersymmetry, are constrained to be of order 10~ 3 . If supersymmetry is invoked to 
tame the hierarchy problem of the Higgs mass, it should provide a reasonable explanation 
for the smallness of these phases. Without a mechanism of this sort, such scenarios of 
supersymmetry are problematic. 

In the Standard Model, the description of CP violation is given by the CKM mecha- 
nism. This mechanism is extremely economic, allowing for one single CP violating weak 
phase. Therefore, within the SM CP is violated "minimally" . This provides us with an in- 
credibly predictive framework, which can be tested meticulously by over-constraining the 
four parameters of the CKM matrix. It is somehow surprising, and certainly remarkable, 
that this simple picture is so far in quantitative agreement with all laboratory experiments 
made up to now. The combination of all observables that constrain the CKM matrix is 
usually carried out through a fit to the apex of the Unitarity Triangle (UT). Today, this 
fit defines a consistent and very much constrained Unitarity Triangle (see Figj2|), giving 
the following results for the real (p) and imaginary (fj) parts of its apex: 

p = 0.164 ±0.029 , fj = 0.340 ± 0.017 (1) 

The general hope, however, is that more precise measurements and more precise theo- 
retical predictions will at some point reveal some inconsistency. Once a number of well 
identified observables have been observed to deviate from the SM expectations, they can 
be used to study the nature of the New Physics. If these are flavor observables, their 
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Fi gure 2: Fit to the apex of the standard Unitarity Triangle as of October 2007. 

deviations will be a holy grail for the understanding of the flavor structure of the new 
physics that might be directly observed at the LHC. 

There is a major theoretical difficulty, however, that affects the determination of the 
SM parameters and the search for New Physics. This difficulty is related to the non- 
perturbative nature of QCD at low energies. Without it, the question of whether the SM 
is in agreement with all the measured observables would be a strict matter of experimental 
precision (combined possibly with a perturbative calculation up to a large enough number 
of loops). 

There are several features of strong interactions that are well known. First, at high 
energies the hadrons show a "partonic" structure, as can be inferred from deep inelastic 
scattering experiments, revealing the existence of quarks. This picture of hadrons in 
terms of "partons" unravels one of the characteristic features of strong interactions, which 
is asymptotic freedom. Second, the mere existence of hadrons means that for some reason 
the elementary quarks from which they are composed, are confined. This is verified by 
the fact that quarks are not observed in isolation, and the fact that colored "objects" 
are not observed either. Although quark confinement and color confinement are slightly 
different concepts, it is clear that evidence favors both. Both features, asymptotic freedom 
and confinement, can be understood qualitatively in terms of the scaling of the strong 
interaction strength. In this picture, strong interactions are very strong at long-distances, 
and they decrease with increasing energy to become very weak at short distances. 

The reasons to believe that QCD is the theory of strong interactions are compelling. 
First, the theory is formulated exactly as it is understood nowadays that a relativistic 
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quantum theory of particle interactions should be. It is a renormalizable quantum field 
theory based on a local gauge principle. Its gauge group is SU(3), where the fact that 
there are 3 colors is known from e + e~ — > hadrons. From the non abelian nature of the 
gauge group it follows that, at least for a reasonable number of fermion families, the 
theory is asymptotically free, and that the coupling constant increases at long distances. 
Moreover, perturbative computations in QCD are very successful at large energies, where 
perturbation theory applies nicely. So, at least in the perturbative regime, QCD shows the 
qualitative features of strong interactions and reproduces quantitatively the experimental 
results. 

In the non-perturbative regime, however, things are less clear. First, it has not been 
proven that QCD implies confinement. Second, it is not known how to extract the hadronic 
spectrum from QCD. For example, the Bethe-Salpeter equation, as a dictionary that 
translates from quarks to hadrons and back, it is extremely difficult to solve. Despite 
these formal deficiencies, however, a great progress has been made since the development 
of QCD in the 70's, and many qualitative features of strong interactions at low energies 
can be connected to features of "QCD-like" theories. For example, the large-N limit of 
QCD is able to give qualitative explanations for the Zweig's rule, the dominance of the 
leading fock states in the mesons, or the success of Regge phenomenology. Other successes 
of QCD itself in the non perturbative regime concern predictions derived for example from 
QCD sum-rules or from lattice simulations. 

Concerning the strong interactions in B decays, the progress has been driven by the ob- 
servation that one can perform a perturbative expansion in the small parameter AQCD/mb- 
The first step was the development of an effective theory for mesons containing heavy 
quarks, called Heavy Quark Effective Theory (HQET). This theory manifests explicitly 
the symmetries that arise in the heavy quark limit and allows, for example, to relate dif- 
ferent heavy-heavy form factors to a single (Isgur-Wise) function. For inclusive B decays, 
the large value of the b mass allows to use the so called Heavy Quark Expansion (HQE), 
which predicts, for example, that the inclusive decay of a B meson is dominated by the 
partonic decay of the b quark alone. In the case of exclusive decays, the establishment of 
a power counting in terms of Aqc d / frib led to the development of the QCD factorization 
approach that has been used extensively to make predictions for all two body hadronic, 
and radiative B decays. Finally, a consistent effective field theory of exclusive and inclu- 
sive heavy meson decays at large momentum transfer has been formulated under the name 
of Soft Collinear Effective Theory (SCET). However, much has to be done before we are 
able to give theoretical predictions for exclusive hadronic decays with uncertainties at the 
percent level. 

In this thesis I present some topics related to my work in non-leptonic decays of B 
mesons. It is divided it two major parts. The first part is an overview of the basic 
matters that constitute the background on which the original work is based. Such a 
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general overview is not an easy task in this field, since there are many topics involved, 
some of them very well established, and some of them not so much. Therefore I might 
have been too extensive on very well known areas, maybe short in topics that require more 
explanation, I may have been repetitive on some parts and on the contrary have skipped 
completely things that some people will miss. However, I believe that the whole text 
is quite self-contained, and that this first part constitutes a rather complete link to the 
second part. I have divided this part in three chapters: first I discuss the weak effective 
hamiltonian in B physics, both in the context of the SM and for the study of NP. Then I 
present the concept of factorization for the computation of matrix elements, and also the 
use of flavor symmetries as an alternative to direct computations. Finally I present the 
theory of CP violation in meson decays. All these three topics are necessary ingredients 
in the work presented in the second part. 

The second part of this thesis contains a series of applications of the theory presented 
in the first part. It is based on some papers that have been published within the last 
three years. These chapters are intimately related, but constitute different "topics" in 
the field; this is what motivates the title of this thesis. The central issue is the study of 
B s — > KK modes within and beyond the SM, but it requires an interconnection with the 
whole grand compact field of flavor physics in general, and in particular of B physics. I 
hope the concepts are written clearly and the structure of this thesis makes it easy to read, 
and easy to pose questions, and comments. 

Finally, I would like to acknowledge my collaborators, David London, Seungwon Baek, 
Sebastien Descotes-Genon and Rafel Escribano for nice collaborations and thousands of 
mails of enlightening content. I am also indebted to the theory groups at Dortmund and 
Rome for the hospitality during my stays there, where also some part of this work was 
done. Many thanks also to Felix Schwab, who so kindly made a thorough reading of the 
first part of this thesis and contributed with intelligent comments and suggestions. I am 
particularly grateful to Guido Martinelli, Bernardo Adeva, Sebastien Descotes, Santi Peris, 
Joan Soto, Gudrun Hiller and Ramon Miquel for reading the thesis, and participating as 
members of the committee. Some of them, specially Guido, contributed with important 
suggestions, improving further the quality of this writing. Very special thanks go to my 
thesis advisor, Joaquim Matias, with whom Fve had endless conversations about physics 
and life since April 2004, and has gone through painful revisions of this text. I can hardly 
imagine a better tutor. Personal acknowledgements have been kept personal. 
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B physics describes weak decays of B mesons. This means that it studies physical processes 
which involve simultaneously -at least- three different energy scales. First, the fact that 
they are weak decays implies that they involve the scale of weak interactions, given by the 
mass of the W boson, M\y. Second, since the energy of the process is that of the decaying 
meson, a second scale involved is the mass of the B meson tub- Third, the fact that we 
are dealing with mesons implies that the physics of strong interactions of bound states is 
also important. The scale introduced in this case is the hadronic scale Aqcd- Moreover, 
since we are looking for new physics, it is reasonable to allow for the existence of at least 
another energy scale, Ajyp, the scale at which the SM breaks down as an effective theory. 

The mere existence of several energy scales does not by itself recall for the use of 
effective hamiltonians. The utility of an effective theory arises when two different physical 
phenomena are mostly independent of one another because they operate at completely 
different scales. For example, the vibration of the atoms in a macroscopic object is too 
fast to affect at all the movement of the object, and the movement of the object is too 
slow to affect at all the vibration of its atoms, so both physical processes can be studied 
independently. Weak decays of B mesons are most suitably studied within an effective 
theory approach because of the large separation between the energy scales involved: 



In fact, Aqcd ~ 0.2 - lGeV, m B ~ 5GeV, M w ~ 100 GeV and A NP > few TeV. The 
large scale of weak interactions and New Physics with respect to the mass of the B meson 
motivates the use of the weak effective Hamiltonian, which is introduced in this chapter. 
The low scale of the strong interaction inside hadrons with respect to the mass of the B 
meson (or the b quark in this case) is what motivates the use of Heavy Quark Effective 



Aqcd < m B < M w , A NP 



(1.1) 
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Theory (HQET). A classical exhaustive reference for HQET is the review by Neubert pQ. 
Other scales that appear in exclusive decays of B mesons are related to the collinear and 
hard-collinear degrees of freedom, that introduce a scale of order mb^-QCD which motivates 
the use of Soft Collinear Effective Theory (SCET). The theoretical and phenomenological 
importance of this recent development calls at least for a set of references linked to its 
formulation [21 HJ U E] . 

1.1 The Weak Effective Hamiltonian 

The study of hadronic weak decays is rooted to the concept of the effective weak Hamil- 
tonian. It describes weak interactions at low energy, relevant for energies below Myy, and 
it has the structure of an expansion of local effective operators (OPE): 



It can be regarded as a generalization of the Fermi theory to include all quarks and 
leptons and the electroweak and strong interactions described by the SM. Furthermore, it 
constitutes a suitable framework for the inclusion of physics beyond the SM. The effective 
Hamiltonian is defined so that the amplitude is given by 



When dealing with such an effective description, there are two ways to proceed. The 
first one is to treat the coefficients in the OPE as unknown phenomenological parame- 
ters, to be measured or constrained through experiment. Once these are measured, the 
effective theory can be tested by checking its predictions over different observables. This 
approach was the one followed in Fermi's theory, and it is relevant when the fundamental 
theory is unknown. Yet the SM itself is often seen as an effective theory in this context. 
Alternatively, if the underlying theory is known (or assumed) then these coefficients can 
be calculated in terms of the fundamental parameters. 

When following the later approach, one must compute the amplitudes in the full theory. 
One will, in general, encounter the typical divergencies that can be absorbed through 
conventional renormalization, and that introduces a renormalization scale \x. When dealing 
with the effective theory, the inclusion of QCD effects also introduces divergencies. Some of 
them can be reduced through field renormalization; however, the resulting expressions are 
still divergent. Hence, one is forced to introduce an operator renormalization, to remove 
these divergencies. This process generically mixes different operators in the OPE, in such 
a way that new operators can arise that were not present at tree level (without QCD 




(1.2) 




(1.3) 
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corrections). It is therefore important to work with a basis of operators which is "closed" 
under renormalization. 

There are several very important issues related with the appearance of a renormaliza- 
tion scale in the OPE. First, the full amplitude cannot depend on the arbitrary scale //. 
As a result, the //-dependence of the coefficients Cj(/x) has to cancel the //-dependence of 
the matrix elements {Qi (//)). This cancellation involves in general several terms in the ex- 
pansion in (|1.3p . Since this scale can be chosen freely, it is a matter of choice what exactly 
belongs to Cj(/i) and what to (Qi (//)). In general, by running the value of fi one assigns 
to these two quantities different energy ranges, such that Cj(/t) contains short distance ef- 
fects above \x and (Qi([i)) contains the long distance non-perturbative contributions with 
energies below /t. 

A second comment on the renormalization scale has to do with the perturbative regime 
of QCD and the suitable scale to "compute" the matrix elements of the operators. As the 
coefficients are calculated perturbatively, they should be computed at a scale in which the 
QCD coupling constant is small. This means a scale ~ O(Mjv). However, the scale 
at which the matrix elements can be factorized in a meaningful way is much lower, for 
instance /// ~ Oiyn^, m c , 1 GeV) for B, D, and K decays respectively. The problem that 
arises here is that the large difference of the two scales involved spoils perturbation theory 
through large logarithms. Fortunately, the Renormalization Group (RG) provides a tool 
to recover the validity of the perturbative series by a resummation to all orders of large 
logarithms. 

In this chapter we will review schematically the specific details that allow to realize 
these ideas in a quantitative way. A thorough but pedagogical treatment of these matters 
can be found, for instance, in [U [9]. They will play a central role in the forthcoming 
chapters. 

1.1.1 Operator Product Expansion in Weak decays 

The basic idea that motivates the use of low energy effective theories is that there is a 
lower limit on the distances that can be resolved through a process of given energy. This 
means that short-distance effects below this limit ~ 1/fc, with k an external momentum, 
can be treated as local. In particular, one can "integrate out" heavy modes, in such a way 
that non-local interactions mediated by heavy particles are reduced to local interactions. 
A simple example of this process is shown in Figure [LT1 where the W boson is integrated 
out to give a local four-quark operator, Qi- 

To be specific, we show how this works in the case of weak charged current interactions. 
The relevant part of the weak interaction Lagrangian density is 

Ccc = ~(d„W+- d v w+)(d»w- v - d v W~») + M^W+W~^+ J?-(J+W + »+ j-w-*) 
2 2v 2 

(1.4) 
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Figure 1.1: (a) Leading order contribution to u — > u'dd! in the full theory, (b) Contribution in 
the effective theory. Integrating out the W boson in the full theory generates the local operator 
Q2 = (du) v _ A (u'd') v _ A in the effective theory. Here, u and u' denote u-type quarks and d, d! 
denote d-type quarks. 

which, integrating the W fields in the generating functional, leads to a non-local action of 
the form 

S cc = ~jj d 4 xd 4 yJ;(x)A^(x,y)J+(y) (1.5) 
where is the W propagator. In the unitary gauge, 

The expansion in powers of 1/Af^ is meaningful for energies well below M^, so this step 
identifies the low energy expansion. Inserting this propagator in (jl.5p and reading off the 
lagrangian density we have, keeping only O (1/Mjy) terms, 

£ cff = -^J-J+^) = -^| Y, V^V u , d ,{du) v _ A {u>d>)y_ A (1.7) 

u,u' ,d,d' 

where Gf = fff/8M^. This la grangian no longer contains the W degrees of freedom, as 
they have been integrated out. It has the structure (jl.2p of a set of local four quark opera- 
tors, with coefficients which contain information on the short-distance physics, mainly the 
mass of the W boson. This expansion is called an Operator Product Expansion (OPE), 
and the coefficients of the operators are called Wilson coefficients. Often, the top quark 
appears in loops and must be integrated out too, introducing in the Wilson coefficients a 
rrit dependence. In Chapter [7] we shall be integrating out supersymmetric particles, and 
the coefficients in the OPE will contain the masses of these particles. 

This is an example of how to calculate the Wilson coefficients in (|1.2j) from the under- 
lying theory. For theories beyond the SM one shall begin with a more general Lagrangian 
than (|1.4p . which may in general introduce a different set of local operators, and whose 
coefficients will contain new parameters, couplings and masses. 
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Figure 1.2: Renormalization of the local operator Q^. The vertex corrections such as the third 
diagram cancel the contributions from the field-strength renormalization. The counterterms in the 
second line are thus defined as to cancel the divergencies from the other diagrams, and to satisfy 
the renormalization condition introducing the scale (i. Note that it is necessary to include the 
counterterm of the operator Q\. This is a consequence of operator mixing. 

Up to some point, however, we must introduce QCD corrections to go beyond tree level. 
In the effective theory, this corrections will renormalize the local operators (see Figure [L2]) . 
Since these operators have dimension greater than four, the field renormalization does 
not cancel all the divergencies, and an additional operator renormalization is necessary. 
Furthermore, the inclusion of QCD corrections generates new operators with new color 
structure (or new flavor structure due to Fierz rearranging, if preferred). This operator 
mixing is responsible for the renormalization constant being a matrix, 

QT = ZijQj (1.8) 

The important fact about the operator renormalization is that it does not depend on 
the full theory; once we have renormalized a complete set of operators to find Zij, this 
will be used for the effective theory of any fundamental theory that generates this set of 
operators. 

Since the renormalization scale dependence introduced in this step will play a central 
role, it is rather more convenient to introduce the anomalous dimension matrix of the 
operators, defined as 

l = Z^-^—Z (1.9) 
am /i 

Now we give the basic guidelines to calculate the leading order anomalous dimension 
matrix for a general set of local operators. 

1.1.2 Calculation of the Anomalous Dimension Matrix 

In order to compute the i'th row of the matrix Z^ to leading order we must compute the 
one- loop QCD corrections to the operator Qi, as shown in Figure [L2l For the sake of 
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generality (focusing on four-quark operators for simplicity), we denote this operator as 

Qi = {d^vfup) ® (n' 7 r 2 y/4) (1.10) 

where V" 13 (8) V^ 6 refers to the color structure, and T\ (g) r 2 to the Dirac structure [S]. 
We must calculate the last three diagrams in the fist row of Figure 11.21 (together with 
the symmetric ones not shown) with the insertion of this general operator. We denote 
these (pairs of) diagrams by V a , V^ and V c . Because we are interested in computing 
the Z constants, we just need the divergent parts of these diagrams. A straightforward 
calculation (in dimensional regularization) gives 

Va = £i7 ( c - 1)r i» ® r 2 7 p 7 At + 4 2) 7M7pri ® yyr 2 ) 
Vb = Sl7 { c b ) wFi'f'f ® r 2 + cf ri ® 7M7pr 27 p 7 M ) 
Pc = Si7 (ci 1) r 1 7 P 7^7Vr 2 + ci 2 ) 7/ , 7p r 1 ®r 27 v) (1.11) 

where the color factors are given by 

4 1 ) = V\T a ® v 2 T a ci 2) = T a v x ® r°y 2 

c (i) = T « yi T a ® y 2 } = Vi <8> T a y 2 T a 

= V\T a ® T a V 2 C^ 2) = T a V x ® V 2 T a (1.12) 

The operators appearing in (jl.lip can be transformed into the operators (Qj) if the set 
closes under renormalization. The sum of these divergent parts of the diagrams give the 
counterterms of the operators plus the counterterm corresponding to the renormalization 
of the quark fields, 

]T Vt = {Z- 2 - 1) (Qi) + (Zij - <%) (Qj) (1.13) 

i 

Since the dependence of Z on [i comes only from a s , the leading order anomalous dimension 
matrix (jl.9p is 

7<°>g = -2(Zy-dy) (1.14) 

with (Zij — 5ij) obtained from (|1.13p and (jl.lip . 

For example, the leading order anomalous dimension matrix for the set of operators 
p.28p is given in Table [TT1 The computation of the anomalous dimension matrix beyond 
leading order is a much more formidable task, especially when, besides QCD corrections, 
the desired precision requires electromagnetic or even electroweak corrections. 
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Table 1.1: The anomalous dimension matrix at leading logarithms for the set of twelve operators 
0i,...,io, 77 , °&g defined in eq. (flT28|) [TO]. 



1.1.3 Renormalization Group Evolution 

Up to this point, let us recall the basic procedure discussed so far. The basic steps in the 
calculation of amplitudes in weak decays are the following: 

• Calculation of the amplitudes within the full theory. This is done at a high scale 
{nw) suitable for a perturbative treatment. In this step one can see which are the 
operators generated at this scale, and write the relevant OPE for this theory. 

• Renormalization of the effective operators and computation of the anomalous di- 
mensions. Once this is done, the full theory can be matched into the effective theory 
to find the Wilson coefficients. These will be free from divergencies, but will contain 
logarithms of the form log (Myy/Hyy)- This /i dependence is cancelled by the logs of 
the form log(//jL/s) in the matrix elements (Qi). 

• Calculation of the matrix elements. As we will see, this has to be done at a low 
scale, and thus the \x scale has to run down from Mw to the lower scale, which for 
the five flavour effective theory will be m&. 

The last step introduces a subtlety. The renormalization scale is arbitrary and therefore 
one could think that after computing the Wilson coefficients it can be set to without 
further care. However, this introduces the large logarithms log (My^/m^), which appear 
in the perturbative expansion as [a s log (M^//i 2 )] n , a s [a s log (M^ / //x 2 )] n , etc., and spoil 
the convergence of the perturbative series. Fortunately, the Renormalization Group allows 
to sum the terms [a s log (M^/pi 2 )] n to all orders in perturbation theory at leading order, 
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the terms a s [a s log (M^r / /f 2 )] n at next-to-leading order, and so on, yielding what is called 
the RG improved perturbation theory. 

Consider the necessary /i-independence of the amplitudes, in the form of a RG equation, 

i 

Recalling the ^-independence of the bare matrix elements (Qi)^ ' and the definition (jl.9p . 
equation (|1.15|) can be written as a RG equation for the Wilson coefficients, 

^CV)=7 T (/^(M) (1-16) 

where for simplicity we use vector notation. This differential equation can be solved rather 
trivially at leading order, to give 

fT(0) 

'a s {M w y 



C(ji) 



C(M W ) (1.17) 



a 8 (fj,) 

where we have used the definition of the /3-function at leading order, given by da s /d\n fj, = 
— 2/3o(as/47r). In general, beyond leading order, one can still write 

C(ji) = U 5 (jJL,Mw)C(M w ) (1.18) 

but the /U-evolution matrix Us((jL,M\y) takes a more complicated form than in (I1.17|) . 
Furthermore, when running the renormalization scale below the threshold of nib, one 
should start integrating out the b quark to go to the four-flavor effective theory, and the 
charm quark to go beyond m c , and so on. 

By looking at the solution (11.17|) of the RG equation, one can immediately see that an 
order a\ matching of the Wilson Coefficients provides in fact terms of all orders in a s . A 
more careful inspection shows that the terms are of the type [a s log fi] n with n any positive 
power. It is, however, far from evident that higher loop computations do not contribute 
terms of this type. This is indeed the case, and that is why after a one loop matching a 
resummation of all the terms [a s log//] n is possible. Here I give a recursive proof. 

Write the Wilson Coefficients as an analytic expansion in powers of a s and £ = log /i, 

^) = E%(^)V (i.i9) 

with arbitrary coefficients Sij. Consider also the perturbative expansion for the beta 
function and the anomalous dimension matrix, 

dota ( a<t\ 2 



-8vr 



(s)5>(£)' ^ 
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It is convenient to let i, j G Z and take a™ = = y r ^ i > = for i,j < 0. This takes into 
account that (da s /d£) appears first at second order in a s and the anomalous dimension 
at first order. Now, because the Wilson Coefficients depend on i explicitly and implicitly 
through a s , one usually writes the RG equation as 

d +^4-\e=?e (i.22) 



di dt da 



Now, by plugging into this equation the expansions for (da s /d£), 7 T and C in powers of a s 
and £, and requiring equality at each order, one arrives to the following recursive relation 



for the coefficients Oy, 

(J + 1) " 8vr + 1) a i+1J (3 k _i_ 2 = £ a itj ^ k ~^ V j, k (1.23) 

i>0 i>0 

From this recursive relation one can extract all the coefficients a^- from initial conditions 
a n Q (that is, the Wilson coefficients at the matching scale). The important thing is to 
see what coefficients Pi and are required for each a^. All the necessary results to 
understand how these coefficients are related can be proved by induction from fll.23H . For 
example, the first result is that 

a i:j = Vj" > i (1.24) 

which tells us that the power of the logarithms cannot exceed the order in a s . The next 
thing one can prove is that 

a k k = \ (t T(0) + 8vr(A: - l)/3 ) 4-i,fe-i (1-25) 

This is a very important result. It tells us that all the terms of the form [a s log/i] n in 
the Wilson Coefficients are known from the matching condition C(M\y) and the leading 
order coefficients 7 T (°) and (3q. This is exactly what allows for a resummation of leading 
logarithms (LL). It is essential for the convergence of the series in f 1 1 . 1 9 1) when the difference 
between the two scales Myy and fi makes [a s log (Myy/fj,)] a number of order one. 

Next to leading logarithms (NLL) can also be resummed in a similar fashion. For 
example, all the terms of the type a s [a s log /i] n are known once one has computed the beta 
function and the anomalous dimension matrix at next to leading order. The corresponding 
recursion relation from which this fact follows is 



1 

ak,k-l 



k - 1 



1.26) 



and thus requires also a one-loop matching of the Wilson coefficients (notice that besides 
aoo also aio is needed). 

One can proceed in this way proving more complicated recursion relations. In general, 
for a resummation of terms of order ct"[a s log ji] k for all k one needs: 
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• Sqo, aiOi • • ■ , S n o (n-loop matching of the Wilson Coefficients). 

• 7 T (°\ A f T ^ 1 \ • • • , j T( - n ^ ((n+l)-loop anomalous dimension matrices). 

• Po> Pi, ■ ■ ■ i Pn (beta function at n+f loops). 

Of course, in practice the resummations are not done recursively. One just solves the RG 
equation differentially as in (|1.17p . In my opinion, however, the recursive proof is more 
transparent and allows to keep control on the required orders in perturbation theory. 



1.2 The Weak Effective Hamiltonian in the SM 

In the SM, only a reduced set of operators are generated in the weak effective Hamilto- 
nian. Here we show only the effective Hamiltonian relevant for |AJB| = 1 and |A£7| = 0, 1 
processes. In this case, the operators are usually classified as current-current operators, 
QCD-penguin operators, electroweak penguin operators, and electromagnetic and chromo- 
magnetic operators. The Wilson coefficients always contain the same products of CKM 
elements, so these are factored out and written explicitly in the effective Hamiltonian, 
denoted by Xp = V p bV* D . Here D = d for AS = and D = s for AS = 1 decays. The 
effective Hamiltonian is given by 

Hcff = % S X v ] ( Cl <2i + C 2 Ql+ Ci & + C ^ ^ + C *9 Qs 9 ) + h.c. , (1.27) 

v ^ p=u,c ^ i=3,...,10 ' 

where Q\i are the left-handed current-current operators, Q3...6 and Q7...10 are the QCD 
and electroweak penguin operators and and Q$ g are the electromagnetic and chromo- 
magnetic dipole operators. Their explicit form reads 



Qi 


= (pb)v- 


A (Dp) V -A , 


Ql 


= (pibj) V -A(DjPi)v-A, 


Qs 


= {Db)v- 


-A^2q{qq)v~A, 


Qa 


= (qjqi)v-A , 


Q 5 


= (Db)v~ 


-A^2 Q (qq)v+A , 


Qe 


= {Dibj) v _ A (qjq^v+A , 


Qi 


= {Db)v- 


-A^ q ^e q (qq)v+A, 


Qs 


= {Dibj) V -A %e q (qjqi) v+A 


Q9 


= {Db)v- 


-A^2q^e q (qq) V - A , 


Q10 


= {Dibj) v _ A %e q (qjqi) V -A 


Q77 




Q8g 


= _il mb Da flu (l+ l5 )G^b, 



(1.28) 

where (g~i q<z)v±A = 917/^(1 ± 75)^2 and the sum runs over all active quark flavours in the 
effective theory, i.e. q = u,d, s,c,b. If no colour index i, j is given, the two operators are 
assumed to be in a colour singlet state. 

The Wilson coefficients at NLO are given in Table [L2l at six different scales of interest 
(we will neglect electroweak penguins). These will be used in latter chapters. 
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Ci 


c 2 


c 3 


c A 


c 5 


c 6 






A* = 


2rrib 


1.045 


-0.113 


0.009 


-0.025 


0.007 


-0.027 


-0.281 


-0.136 


A* = 


m b 


1.082 


-0.191 


0.013 


-0.036 


0.009 


-0.042 


-0.318 


-0.151 


A* = 


vn b /2 


1.139 


-0.296 


0.021 


-0.051 


0.010 


-0.066 


-0.364 


-0.169 


A 4 = 


y/2m b A h 


1.141 


-0.230 


0.021 


-0.051 


0.011 


-0.067 


-0.366 


-0.170 


A* = 


y/m b A h 


1.184 


-0.371 


0.027 


-0.062 


0.011 


-0.086 


-0.395 


-0.182 


A* = 


y/m b A h /2 


1.243 


-0.464 


0.035 


-0.076 


0.010 


-0.115 


-0.432 


-0.198 



Table 1.2: NLO Wilson coefficients of the operators 0\ 6> 07 7 and 0% g at six different scales of 

interest. The scale Ah = 0.5 GeV. 



1.3 The Weak Effective Hamiltonian beyond the SM 

In the presence of NP, many more operators can be generated other than those in (jl.28p . 
If NP exists at a scale A^rp, the exchange of these new particles between the low energy 
SM fields will in general generate standard and non-standard operators, with Wilson 
coefficients suppressed by the new physics scale Ajyp: 

oo [D] 

n « = E E TT^-M D] (1-29) 

D=5 i A NP 

Indeed, the effective theory approach is the most general and model independent way of 
dealing with NP at low energies. It relies on the following assumptions: 

• The low energy degrees of freedom are the SM particles with masses below Myy, so 
these fields are the only building blocks in the operators in the OPE. 

• All the operators must be gauge invariant, the gauge group being that of the SM. 
Other global symmetries might be imposed in order to avoid proton decay and other 
phenomenological requirements. 

• A prescription is taken to cut the infinite set of operators. This prescription consists 
in dropping the operators with dimension higher than a given dimension d. Operators 
of dimension D are suppressed by A^ffp, so if the scale Atvp is relatively high, the 
"cutting" dimension d can be low, leaving a reasonable number of operators in the 
Hamiltonian. In the case of B decays this means to keep only operators up to 
dimension six. 

Following this procedure, however, one finds that the most general effective Hamiltonian 
contains too many operators. And since the Wilson coefficients are just phenomenological 
unknown parameters of the effective theory, the theory has too many unknowns. The 
general set of operators was derived in [TT]. The concept of a general effective Hamiltonian 
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will be used in Chapter [J] to write a general parametrization of the NP amplitudes for 
b — > sqq transitions, where the issue of the large number of operators is circumvented with 
an argument concerning the NP strong phases. 

A general structure of the new physics is, as mentioned before, difficult to reconcile with 
a relatively low value for Anp- Flavor physics constrains severely many NP contributions 
to flavor violating operators and CP violating phases in the NP Wilson coefficients. This 
leads to effective scenarios of new physics that reduce considerably the complexity of the 
general effective Hamiltonian and fulfill the experimental bounds almost virtually. An 
example is given by the minimal flavor violation (MFV) hypothesis, which assures that 
the flavor violation in the effective theory comes entirely from the CKM matrix [12] . While 
there are several non-equivalent definitions of MFV in the literature, they all share this 
property. 

This formalism is also useful to study model independently concrete scenarios of new 
physics. For example, there are several models that describe the origin of supersymmetry 
breaking. While these models might be wrong, low energy supersymmetry might be right, 
and the study of supersymmetry at low energies should not rely on the specific scenario 
chosen, for example, at the Planck scale. For a derivation of a low energy supersymmetry 
effective Hamiltonian see |13j . 

For model dependent studies of specific new physics, the use of the weak effective 
Hamiltonian is much the same as in the SM. The only difference is the extended set of 
operators that is generated, and of course that the matching conditions lead to different 
Wilson coefficients. In Chapter [7] we will match gluino-squark box and penguin contribu- 
tions to b — ► sqq transitions into the NP weak effective Hamiltonian. 



Chapter 2 



Hadronic Matrix Elements 



The concept of the Operator Product Expansion, discussed in the previous chapter in the 
context of the weak effective Hamiltonian, is a powerful tool. It allows for a factorization 
of long and short distance physics; the short distance contained in the Wilson Coefficients 
and the long distance contained in the matrix elements of the operators. The Wilson Co- 
efficients are calculable within perturbation theory at a high energy scale (e.g Myy), and 
are process-independent. Moreover, this large scale can be dragged down to a phenomeno- 
logically sensible scale (e.g m&) by means of the RG resummation of large logarithms of 
the form log(Mw/m&) that would spoil the perturbative series. 

In this way, a general process-independent effective Hamiltonian is derived, which in- 
cludes QCD (and possibly electroweak) corrections to a given order (LL, NLL, NNLL,. . . ) 
from gluons of virtuality above the scale [i. The effective Hamiltonian depends on fi 
through the Wilson Coefficients. All this is process independent, and the model depen- 
dence comes solely from the matching conditions of the WC (this can very well account 
for the fact that many operators might have been dropped from the beginning: one just 
includes all the operators, but puts to zero the WC's of those operators that are not 
wanted). 

In order to compute the amplitudes, this effective Hamiltonian must be sandwiched 
between the initial and final states (c/. (|1.3j) ). 



The last step of the process is then the computation of the Hadronic Matrix Elements of 
the operators (Qi(fi)), between the initial and the final states. 

The scale dependence of the matrix elements is conceptually subtle but enlightening. 
The operators, by themselves, do not have any scale dependence. Neither do the initial or 
final states. However, it is clear that the matrix elements must have a scale dependence 




(2.1) 
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that cancels the scale dependence of the WC's. (The same is also true for the scheme 
dependence, which leads to a criticism to some generalized factorization approaches as 
will be discussed later.) The scale dependence of the matrix elements can be partially 
understood through the following physical reasoning. Since the effective Hamiltonian is 
process independent, the matrix elements are precisely what distinguishes between dif- 
ferent final states (lets say the initial state is always a B meson). But the final state is 
not characterized only by the quark fields of which is composed. Indeed, one can always 
rearrange the quarks in the final state to get different hadronic final states. This rearrange- 
ment is due to strong interaction rescattering. The point is that QCD corrections with 
virtualities above the scale \i (hard gluons) are accounted for in the Wilson Coefficients, 
which are process independent, so that these corrections have no impact on the final state 
hadronizatior0- Therefore, the final state hadronization takes place through interchange 
of soft gluons (QCD corrections with virtualities below //), which must then be contained 
inside the hadronic matrix elements. 

Up to now, however, there is no fully consistent way to compute hadronic matrix ele- 
ments from QCD (apart from, arguably, lattice methods). This is a longstanding problem 
that goes back to the foundation of QCD and that partially motivates this thesis. 

The concept of factorization of matrix elements, as a second step after the OPE, has 
proven a powerful tool to face this problem. The idea is to reduce the matrix elements to 
products of form factors and decay constants, which are process independent quantities. 
However, even in the cases where factorization can be strictly justified, it just solves the 
problem partially, since a systematic way of computing form factors and decay constants is 
also necessary. This is a less pathological problem, though, since their process-independent 
nature allows to extract them from data. Also, because they are intrinsically simpler 
objects, it is possible to extract them from lattice simulations, which is up to now the 
most pure "QCD-based" technique for non-perturbative computations. 

Factorization is strict in the case of leptonic and semi-leptonic B decays. The simplest 
ones -concerning QCD corrections- are the leptonic decays, which are those with only 
leptons in the final state. Because leptons do not interact strongly, QCD effects take 
place only within the B meson in the initial state. Therefore, when computing the matrix 
elements, factorization is strictly valid, and all the long-distance strong interaction effects 
are contained in 



which defines the decay constant fs- 

Semi-leptonic decays are those with leptons and hadrons in the final state. Strong 
interactions take place inside the initial and final hadrons, in the process of hadronization, 
and between them during the decay. But these corrections affect only one vertex of the 

lr The values of the WC definitely tell what are the most favorable final states, but that is a feature of 
the structure of the interactions themselves, not of the hadronization process. 
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weak current (since the other vertex is leptonic), and thus factorization is still exact. The 
non-perturbative physics is parameterized now in terms of a form factor, 

(M q \bTq\B)~F B ^ M (2.3) 

The case of non-leptonic decays is more complicated. Since the final state is composed 
purely of hadrons, factorizable and non-factorizable effects take place. Therefore, factor- 
ization is no longer justified. The observation that, surprisingly, strict (naive) factorization 
works reasonably well also for non-leptonic decays, led originally to what has become an 
intensive line of research. Now we understand the reason, since we have factorization 
theorems that tell about the limitations of naive factorization. 

An alternative approach to deal with non-leptonic B decays beyond factorization is to 
take advantage of approximate symmetries of QCD. Flavor symmetries arise in the limit 
in which the masses of the quarks are equal, and it turns out that for the three light 
quarks u,d,s this is a good approximation. This introduces a SU(3) V symmetry group 
under which the QCD lagrangian is approximately invariant. Symmetries always allow 
to establish relations between matrix elements, according to the Wigner-Eckart theorem. 
Therefore, up to symmetry breaking corrections, SU(3) relations between amplitudes of 
different processes hold, and allow to make predictions without the need of computing 
matrix elements. Of course, how big the symmetry breaking corrections are is a dynamical 
question. 

In Section 12.11 we introduce more formally the idea of factorization for two body non- 
leptonic B decays, and give the relevant formulae for the evaluation of amplitudes within 
the framework of QCD-Factorization. In Section 12.21 we explain the main features of the 
use of flavor symmetries in hadronic B decays. Both subjects will be important in the 
development of the following chapters. 

2.1 Factorization 

2.1.1 Beyond naive factorization 

Consider a decay of a B q meson into two mesons M\ and Mi- We would like to calculate 
the matrix element 

(MiM2|Oi|Sg> . (2.4) 

The meson M2 can be generated directly by a quark current containing the appropriate 
flavor and Lorentz quantum numbers, say <?i7/x75<?2 for a pseudoscalar. A dimension- 
six operator Oi containing this particular current will contribute to the decay through 
a factorized product of two current matrix elements if the other current has the correct 
quantum numbers to generate the transition B q — > M\. Let's say this current is q'^fj). So 
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factorization amounts to the following simplification, 

(MiM 2 1 (^7592) («?V)I^) __ C -U (M 2 \qi7^ 5 q2\0){M 1 \^b\B g ) ~ f M2 -F B ^ A,h 

(2.5) 

If this simplification is justified, then one can express the matrix elements such as (j2.4|) in 
terms of products of a decay constant and a form factor. This oversimplified presentation 
gives an idea of the concept of factorization. The strongest version of this process consists 
simply to promote (|2.5p to an equality, and it has been called naive factorization. 

In the weak effective Hamiltonian, one in general will have to deal with pairs of oper- 
ators of the following type 

Oi = (q a b a )r 1 {qipq2/3)r 2 , 

Oj = (qabp)r 1 (qipq2a)r 2 ■ ( 2 -6) 

where a, (3 are color indices and Ti, r 2 are Lorentz structures. These operators can be 
Fierz rearranged and written as 

O'i = {qipb a ) r/i (q a q 2 p) r > 2 , 

Oj = {qi a b a ) T > i {qpq 2 p)r' 2 • (2.7) 

where in general the Lorentz structures r' l5 T' 2 are different from T\, r 2 . Now, depending 
on the flavor structure of the decay, it is possible that the operators can be factorized in 
the form O.;, in the reordered form Oj, in both forms or in none of them. In the case 
that both forms contribute, both have to be taken into account since they correspond to 
different quantum mechanical paths (different ways to reorder the quarks to form the final 
mesons). The factorized matrix elements are defined as 

(Oi) F = (M 1 |(^)rjB 9 )(M 2 |feg 2 )rjO) + (M 2 |(g6) ri | J B ? )(M 1 |(g ig2 )r 2 |0) , 
(0'j) F = (Milfa&VjBgXMaK?^ (2.8) 

In order to introduce factorization it is useful to work in the singlet-octet basis for the 
operators. Using the identity for the SU (3) generators 

(2.9) 
one can write 

aOi + CjOj = (d + 1c,-) Oi + 2Cj of 

dOi + CjO'j = (Cj + ^C i )o' j + 2C j O'f (2.10) 



rpa rpa 
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where Of = (qT a b)r 1 (qiT a q2)r 2 i an d similar for O^. It's clear that in order for the 
octet operators to contribute through a factorized matrix element, strong interactions 
must somehow change the color structure, so it is reasonable to think that these operators 
contribute only at order a s . Be that as it may, in naive factorization, by definition, the 
matrix elements of octet operators are set to zero. At the end, in naive factorization (NF), 
the contribution to the amplitude of this pair of operators is given by 

CiQj.)(Oi) + Cji^iOj) = ai(^){0 % ) F + ajiriiO'^F (NF) (2.11) 

where 

OiOx) = Ci(fi) + iCfOx) , a» = C» + . (2.12) 

The hadronic matrix elements of the currents that result from factorization do not show 
any renormalization scale dependence that can cancel the scale dependence in the Wilson 
coefficients. Therefore, naive factorization cannot be correct. However, one may hope that 
a particular factorization scale [if exists, at which the Wilson coefficients are evaluated, 
for which this is a good approximation. As we will see, there is a major inconvenient to 
this point. 

The attempts to give a formulation for factorization that didn't have the problem of 
the renormalization scale dependence led to different generalized factorization approaches. 
We first present the formulation by Neubert and Stech |14] . 

The structure of factorization can be made exact by introducing process dependent 
non-perturbative parameters e\{\i) and £%(\i) that parametrize the non-factorizable contri- 
butions. For example, in the case of a decay with a flavor structure such that {Oj)f = 0, 
they are defined as 

so that Ci(Oi) + Cj(Oj) = [(Ci + Cj/N)(l + e\) + esCj\{Oi) F is exact. Moreover, naive 
factorization arises in the limit £i,£s 0. The formula (|2.1ip is then generalized to 

Ci(jj,)(Oi) + CjfjiXOj) = af(Oi) F + afiO'^p (GF) (2.14) 

where the effective parameters are given by 



eflf 



af 



Ci(») + ^Cjfjj)) (1 + eiGu)) + e^)Cj(ji) 
1 

N' 



+ Tt C ^) i 1 + £ 'i^)) + £ M C ii») ( 2 - 15 ) 



Since the factorized matrix elements do not have any scale dependence, all the fi- 
dependence is contained inside the Wilson coefficients and the parameters £j(/i). As 
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the full amplitude is //-independent, the hadronic parameters ej(/t) restore the correct 
//-dependence of the matrix elements, which is lost in the naive factorization approxima- 
tion. This observation is enough to extract the //-dependence of the hadronic parameters 
by means of the RG equations. The condition for this amplitude to be independent of the 
scale is that, for any particular scale fio, 

C i (//)+-C y (/i)) (l+ei(/i))+e 8 (/i)C j ( / u) = (c i ( Ato )+-^C J (/xo)) (l+£i(// )) W/^Guo) 

(2.16) 

holds. Now, using the RG equation (| 1 . 18[) we can write the Wilson coefficients at the // 
scale in terms of those evaluated at /io- Taking into account that the evolution matrix 
C/(/to,/u) cannot depend on the values of the Wilson coefficients, we find 



1 + eiO) = (ku + — «ji)(l + £i(Mo)) + Kjieg(n ) 



e 8 (p) = (/% - («« - Kjj) - ^ i 1 + £ i(^o)) + ( K Jj - K ii) £ s(^o) (2.17) 

where k = U(fj,o,fi). The main idea in generalized factorization is to extract from data 
the non-factorizable parameters £i(m{,), E%{m\,) and then run these parameters according 
to (|2.17p to find the factorization scale \x$ for which they vanish. 

There is, however, a major drawback to the generalized factorization procedure as 
presented in [Uj. While the scale dependence is properly taken into account, at next- 
to-leading order in the renormalization group improved perturbation theory the Wilson 
coefficients depend also on the renormalization scheme. This scheme dependence can 
only be compensated by non-factorizable scheme dependent contributions in the matrix 
elements. It has been proven [15] that for any chosen scale fj,f = 0(mb) it is possible to 
find a renormalization scheme for which the parameters £i ; s(///) vanish simultaneously. 
Therefore, the factorization scale is scheme dependent and can have no physical meaning 
as such. 

A different approach to generalized factorization that doesn't suffer from the scheme 
dependence problem is the one discussed in [HI [EH HHJ H9] . The idea is to calculate in 
perturbation theory the matrix elements of the operators between the quark states. In this 
way one can extract the scale and scheme dependence of the hadronic matrix elements. 
Combining these scale and scheme dependent contributions with the Wilson coefficients 
Ci{fx), one obtains effective coefficients Cf s that are scale and scheme independent. Then 
one can write 

Ci(Oi) + CjiOj) = Cf{Oi) Ucc + Cf(Oj) tree (2.18) 

where (Ojj) tree denote tree level matrix elements. Once this is done, factorization can 
be applied to the tree level matrix elements without the problem of the scale and scheme 
dependence. The result can be cast in the form of eq. (|2.14p . with a?^ calculable up to a 
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phenomenological parameter (an effective number of colors) that should in principle give 
information on the pattern of non-factorizable contributions. 

This approach, however, has also its own drawbacks. As noted in [15], the effective Wil- 
son coefficients are in general gauge dependent, and also depend on an infrared regulator. 
These dependencies originate in the perturbative evaluation of the scheme dependent finite 
contributions to the matrix elements, that are necessary to cancel the scheme dependence 
in the Wilson coefficients. 

Some papers have been written attempting to solve the problems that arise in gener- 
alized factorization approaches (sec for example |20j). Here wc will focus on the so called 
'BBNS' approach or QCD- factorization (QCDF), presented initially in [21] in the context 
of B -> 7T7T, and extended later to general B q -► M 1 M 2 decays [221 EH E3 ESI EH W\ ■ An 
overview is postponed until Sections 12.1.31 and 12.1.41 

We finish this section writing a general formula for the amplitudes in SM within any 
generalized factorization approach in terms of the coefficients a| ff . The SM operators in 
eq. (ll.28p can be rearranged to fit the form of eq. (j2.14| ). 



Qi 


= {pb)v-A{Dp) V -A, 


Qi = 


-- (Db) V - A (pp)v-A , 


Q 3 


= ^2,1 (Db)v-A(qq)v~A , 


Q'a = 


X// (qb) v ^ A (Dq) V -A , 


Q 5 


= E« {Db) V - A (qq)v+A , 


Q' 6 = 


-2^ q (qb) S -p(Dq)s+P, 


Q 7 


= ^Zql e q{Db) V - A {qq)v+A, 


Q's = 


- 3e q (qb)s-p(Dq) s+ p , 


Q9 


= ^2q^ e q(Db) V - A (qq)v-A, 


Qw = 


Y,q hi (qb)v~A(Dq) V -A , (2.19) 



Notice that the difference between the contributions from operators Q3 and Q5 in their 
factorized form is just a minus sign if the qq meson has odd parity. The same is true 
for Qj and Qg. The difference between the contributions from operators Q' 4 and Q' e is, 
besides the factor —2, a minus sign if the Dq meson has odd parity and a chiral factor 
r x from the scalar dirac structure of the currents (see below). Then, the amplitude for a 
B qs — ► M1M2 decay is given by 

(M 1 M 2 \H eS \B) = £ \P*{M 1 M 2 \TI\B). (2.20) 

p=u,c 

with the transition operator T V A given by 

T\ = afA([q s p]\pD]) + afA([q s D}\pp]) + (af T af ) £ A([g a Z>][?g]) 

q 

+(of ± r x af) A([q s q][qD}) + (of T of) £ I e « A (^D][qq]) 
q q 

+ (af ±r x af)Y J le q A{[q s q][qD}) . (2.21) 
q 
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The upper signs correspond to the case when the second meson has odd parity and the 
lower signs when it's even. The matrix elements of the operators A{[- ••][••• ]) are non-zero 
only if the flavor of the mesons match the quarks inside the [• • • ] in either order, and will 
be expressed in terms of form factors and decay constants in the following section. In 
naive factorization the coefficients af are given, for i odd, by 

af F = d + jjCi+i j a f+i = Ci+i + — Ci . (2.22) 

As a final comment, just note that contributions from weak annihilation or hard interaction 
with the spectator quark are not taken into account in this formulation. Naively they are 
expected to be small, but at some point they can have an impact in phenomenology. The 
inclusion of these contributions is an issue in the QCDF approach. 



2.1.2 Form factors, decay constants and meson distribution amplitudes 

In this section we give the definitions for the decay constants and form factors of pseu- 
doscalar and vector mesons, and give the expressions for the matrix elements of the oper- 
ators A{[- ••][■•• ]) that appear in the factorized amplitudes. We also give the definitions 
of the meson distribution amplitudes that appear in the hard scattering kernels in QCDF, 
and present their representation in terms of Gegenbauer polynomials . 

The decay constant fp of a pseudoscalar meson P with 4-momentum q is defined as 

(P(q)\qWq\0) = -ifpqp (2.23) 

For a vector meson V(q, £*) with 4-momentum q and polarization vector e* the longitu- 
dinal (fv) and transverse (fy) decay constants are defined as 

{V(q,e*)\qj^\0) = -i/ V my< , (V(q, e^qa^O) = -i/^(^e* - g„e*) • (2.24) 

Using the Dirac equation for the quark fields, the following identity follows 

d^q^q') = (MYW - qi^{d^q') = i(m q + m q ,) q l5 q' . (2.25) 

Therefore, we have (m q + rn q i){P{q)\q^q l \Q) = —id^—ifpq^) = ifpmp, and 

<P(<Z)I<Z75<Z'|0> = . (2.26) 

m q + m q > 



The scalar matrix element {V(q,e*)\qq'\0) is zero because it can only depend on q ■ e*, 
which is zero as can be easily seen going to the rest frame of the vector meson. 
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For the form factors we use the conventions in [28]. For a B — > P transition, the form 
factors Fq and F + are defined as 



(P(p')\q^b\B(p)) 



(p+P 



2 2 
f\fi _ m B ~ m P 



2 2 

mo — m 



F+(g2)+ ^ »*P q , Fo{q 2 h (2 . 27) 



where q M = (p — p')^ is the momentum transfer. The identity 

d^b) = {d^-fb + <tf(dpb) = i{m q - m b ) qb , 



(2.28) 



and the fact that the F + term in (|2.27p vanishes when contracted with q^, implies that 
for the scalar current 



(P(p')\qb\B(p)) 



m B — rrij 



m b - m q 

For aB-*y transition, we will need the form factors Aq, A%, A 2 and V defined as 

2iV(q 2 ^ 



(2.29) 



(V(p',s*)\q^b\B(p)) 
(V(p',e*)\q^ 5 b\B(p)) 



tub + my 



c livpa * / 



2m v A (q 2 )^^-q^ + (m B + m v )A 1 (q 2 



-Mr) 



TUB + my 



(p+p') 



_ m| m 2 p ^ 



(2.30) 



Now we are ready to write the matrix elements of the operators A{[. ••][... ]) that 
appear in the formula for the amplitudes. We define 



Gp 
V2 



(Mi\(qb) V -A\B q )(M 2 \(qiq2)v-A\0) 



(2.31) 



Using the definitions given above for the decay constants and form factors and neglecting 
terms of 0{jr? M jm 2 B \ we have 

— 1 ^2 m BJP2 P Q 



A T 



A, 
A 
A" 



-2i^m V £* V2 -p B f V2 F^ 



\\P-2 



-2i 



% m v£* Vl -PBfp 2 A^ 



^E-m 2 f A B ^ Vl 



(2.32) 



with all the form factors evaluated at q 2 ~ 0. The case for B — > VV with transversally 
polarized mesons can be found in [291 \27\ ; we omit it because we will be more concerned 
about longitudinal polarizations. Now, we finally have 



{M x M 2 \afA{[ 



B Qs ) = ci a? (MiM 2 ) A UlM2 + c 2 af{M 2 M 1 ) A 



(2.33) 
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where Cj = 0, l,±l/v2 are obvious constants related to the flavor composition of the 
mesons. The chiral factors r x , that account for the difference between (5-P)®(S+P) and 
(y—A)®(y+A) matrix elements are given by 



2m 2 p 



m b (n){m qi +m q2 )(n) 



2my /y(^) 
m b (n) fv 



(2.34) 



neglecting light quark masses with respect to m&. At leading order should be put to 
zero, since the scalar current gives no contribution to the vector decay constant. 

For the meson light-cone distribution amplitudes (LCDA's), the definitions used are 
the ones in [6U\ 131]. We will need only the two particle leading twist (twist-2) LCDA's <3?p 
and <3V, an d the twist-3 LCDA's & p and for pseudoscalar and longitudinally polarized 
vector mesons. 

The definitions for the leading-twist light-cone distribution amplitudes are 



(P(q)\[q(y)wq(x)]\0)\ {x _ yr=0 

(^(9)|[9(y)7^ , (^)]|0)| (a: _ 2/)2= o 



Jo 

ifyq^ [ due^+^y^viu^) (2.35) 



with u = 1 — u, and the brackets meaning that the fields at x and y are connected by 
a Wilson line that makes the matrix element gauge invariant. They are conventionaly 
represented by a Gegenbauer expansion, 



u — u) 



n=l 



(2.36) 



where a^f '(//) are the Gegenbauer moments, and Cn^ 2 \x) are the Gegenbauer polynomials. 
This expansion is usually truncated after n = 2. The relevant Gegenbauer polynomials 



are cf /2 \x) = 3x and Cf'' l} (x) = 3(5x 2 - l)/2. 

The definitions for the twist-3 light-cone distribution amplitudes are 



t(3/2). 



(P(q)\[q(yHq(x)]\0)\ (x _ y)2=0 = -zf Pf i P J due^+^^u, 



(V\\(q)\[q(y)*^(x)]\0)\ {x _ y)a=o = -i/jW / duS^ +u ™\u-u) 



q^Zv qvZfi 



(2.37) 



with z = y — x, and (ip = mp/(m q + m q /). When three-particle amplitudes are neglected, 
<$>p(x) = 1 and 



§ v (u,{i) = 3 



\u — u 



n=l 



(2.38) 
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Here P n (x) are Legendre polynomials. At second order, p2(x) = (3x 2 — l)/2 and Ps(x) = 
(5x 3 — 3x)/2. It should be emphasized that twist-3 distribution amplitudes only contribute 
formally at order 1/rrib to the decay amplitudes. However, they can appear in terms that 
are chirally enhanced by the factors r x , so they must be taken into account. 

An example of how meson distribution amplitudes arise in the computation of matrix 
elements in QCDF is the following. The outgoing meson in the decay has momentum q. 
At leading order, the meson can be described by it's leading fock state: it is composed of 
a quark with momentum uq and an antiquark with momentum uq, with u + u = 1. Then, 
when computing, for example, a vertex correction to the matrix element, one will come 
up with a term of the form 

u aa (uq)T(u, . . . )af3,ab,...Vf3b(uq) (2.39) 

sandwiched between hadronic states. Let's say for definiteness that the meson is pseu- 
doscalar. Then, the prescription is that this term between the vacuum and the meson 
state, must be changed to 



ifp f 1 



AN 



o 



du$ P (u)(4^)p a T(u, . . . ) a p, aa ,... (2.40) 



This prescription is a manifestation of factorization, and requires the proper disentangle- 
ment of long- and short-distance contributions. The fact that this disentanglement occurs 
in the heavy b-quark limit is discussed in the following section. 

The B meson distribution amplitudes arise in the computation of contributions where 
the spectator quark suffers a hard-scattering interaction (hard spectator-scattering terms). 
This is due to the fact that a hard interaction can resolve the inner structure of the B 
meson, probing the momentum distribution of the b and spectator quarks. 

At leading power in 1/rrib the B meson is described by two scalar wave functions. In 
the case in which the transverse momentum of the spectator quark can be neglected in 
the hard-scattering amplitude [22], these scalar wave functions can be defined through the 
following decomposition of the B meson LCDA, 

(0\[q a (z)bp(0)]\B(q)) = + m b H]fr C dCe~^+ z - [$ B1 (£) + 4-*B2{OU 

(2.41) 

where n=(l,0, 0, —1) and the subscripts (+, — , _L) refer to the usual light-cone decompo- 
sition of 4- vectors. Specifically, at the leading order the hard spectator scattering contri- 
bution depends only on the first LCDA $>bi- This dependence is of the form 

( 2 - 42 ) 

which defines a new hadronic parameter A# of order Aqcd [21] ■ 
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2.1.3 Factorization in the heavy quark limit 

Taking advantage of the fact that the b quark mass mj is large compared to the hadronic 
scale Aqcd, it is possible to tackle the problem of hadronic matrix elements in two-body 
B decays in a systematic manner. Indeed, once a consistent power counting in terms of 
A-QCD/irib has been established, it can be proven that factorization is a formal prediction 
of QCD in the heavy quark limit |32| . For this strong statement to be true, however, some 
comments must be added. The proof is valid for a certain type of decays, and the mass of 
the charm quark m c might have to be included in the power counting as a formally large 
parameter, as well as the energy of the ejected meson, which should therefore be light. 
Nevertheless, the importance of such a result is clearly monumental. 

The physical picture is always the same. The most clear example to illustrate this 
picture is the case of a B meson that decays into a heavy and a light meson, being the 
heavy meson who pics up the spectator quark. The b quark in the B meson decays into a 
set of very energetic partons. The charm and the spectator quarks form the final heavy 
meson with no difficulty: since the charm is quite heavy, its velocity will not be large. On 
the other hand, the two light quarks formed in the weak vertex will be very energetic, so 
that if they are going to form a meson they must be highly collinear and in a color-singlet 
configuration. This compact color-singlet object will leave the interaction region without 
interacting with the degrees of freedom that hadronize into the heavy meson, since the soft 
interactions decouple. This is the usual color transparency argument [33] due to Bjorken. 

The fast color-singlet pair of quarks will then hadronize with a probability given by the 
leading- twist light-cone distribution amplitude of the light meson depending on the 

momentum fraction of each quark. The transition from the B meson to the heavy meson 
is parameterized by a standard form factor. This is how the color transparency argument 
links to the concept of factorization. This physical picture was made quantitative with 
the development of the QCDF approach [2"Tj 122]. 

The mathematical formulation of this picture is powerful because it allows to compute 
corrections in a systematic way. For example, it is not required for the ejected pair of 
collinear quarks to be in a color-singlet configuration in order to form a meson. Indeed, a 
hard gluon exchange with the B - heavy meson system can put this pair in a color-singlet 
state before hadronization. The important point is that this correction is calculable in 
perturbation theory (see Fig l2.ip . One must prove, however, that this escaping object 
doesn't interact softly with the interaction region even if it is not in a color-singlet state. 
This is the sort of things that must be checked in order to prove factorization. 

The case of two light mesons in the final state is more complicated. The light meson 
that picks up the spectator quark, receives from the weak vertex a very energetic light 
quark. Therefore, since the spectator quark is slow, the meson is created in a very asym- 
metric configuration. The probability of hadronization is then given by the meson distri- 
bution amplitude near its end point, which leads to a suppression of order (Agc/j/mfc) 3 / 2 . 
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Figure 2.1: Order a s vertex corrections. 










Figure 2.2: Order a s hard spectator interactions. 



Figure 2.3: Penguin contractions. 



Hence, one should take into account the competing contribution from a process in which 
the spectator quark suffers a hard interaction. If this interaction is the exchange of a hard 
gluon with the b quark or with the other quark with who is forming the light meson, this is 
just a contribution to the heavy-to- light form factor. If the hard gluon is exchanged with 
the pair of quarks that form the ejected meson, then this is again a calculable correction, 
called "hard-spectator interaction" (see Fig ]2.2p . 

Other calculable corrections to naive factorization are given by penguin contractions 
(see Fig l2.3p . with the insertion of a penguin operator, or the chromomagnetic dipole 
operator. These contributions introduce complex phases in the hard-scattering kernel 
which account for perturbative strong-interaction phases, due to the rescattering phase 
of the penguin loop. This is basically the familiar Bander-Silverman-Soni mechanism for 
strong phases [Mj- However, this is not the only source of strong phases since the vertex 
corrections in Fig l2.ll also generate imaginary parts. These are correctly included in the 
QCDF approach. 

There is another set of contributions to the decay amplitudes that are completely miss- 
ing in the generalized factorization formula (|2.2ip • They consist of the processes in which 
the two quarks inside the B meson annihilate to form the final state partons (see Fig l2.4p . 
These contributions are formally suppressed in the heavy quark limit. However, some 
of the annihilation topologies related to the corresponding twist-3 light meson LCDA's 
are chirally enhanced by the factors r x , and for realistic b-quark masses the suppression 
Aqcd/to;, might not be sufficient to neglect these terms. Unfortunately, logarithmic end- 
point divergencies do not cancel properly within these terms, leading to the appearance 
of non-factorizable contributions. Some power suppressed terms contributing to hard- 
spectator graphs also suffer from this disease. These corrections constitute one of the 
weak points of QCDF at the phenomenological level. We shall come back to this issue 
later. 
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Figure 2.4: Weak annihilation contributions. 



The final formula for the amplitudes in QCDF is then exactly the same one as Eq. (12.211) . 
except that the coefficients a| ff are not phenomenological parameters anymore, but are 
completely calculable in perturbation theory. Moreover, one should add the hard-spectator 
and annihilation contributions that do not fit in the scheme of Eq. (l2.21j) . The final for- 
mulae for the amplitudes at order a s will be given in the following section. An important 
observation is that naive factorization arises formally from this picture as a prediction of 
QCD at leading order in a s and AQCD/ m b- This justifies mathematically the phenomeno- 
logical successes and limitations of the naive factorization approach. Moreover, beyond 
leading order in a s the scale and scheme dependences are correctly cancelled between the 
Wilson coefficients and the matrix elements. A proof of factorization in QCDF has how- 
ever been only given at order for heavy-light final states [22], and the all order proof is 
not straightforward. Towards this end a more systematic formalism based on an effective 
lagrangian (SCET) is more useful [32] . 

The significant property of the factorization approach in phenomenology is the fact 
that the non perturbative quantities that appear in the formulae for the amplitudes are 
much simpler objects than the original hadronic matrix elements. These simpler objects 
are either related to universal properties of a single meson state, as is the case for the 
light-cone distribution amplitudes, or describe the B — > meson transition matrix element 
of a local current, parameterized by a form factor. These are objects that appear in many 
different decay amplitudes, and can therefore be extracted from data and used to make 
predictions on other modes. Moreover, one can also use QCD sum-rule techniques to 
study them, or extract them from the lattice. Factorization is a consequence of the fact 
that in the heavy quark limit only hard interactions between the l B — ► meson' and the 
'ejected-meson' systems are relevant. 

2.1.4 B -> M 1 M 2 amplitudes in QCD Factorization 

Here we collect the relevant formulae for the B — > M\M 2 amplitudes in QCD Factorization 
according to [261 127] . The amplitudes are given as matrix elements of the transition 
operators and Tg, 

(MiMalWdlB) = Y, X p D> *(MiM 2 \TP + TP\B). (2.43) 

p=u,c 
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The transition operator contains the contribution from the leading order tree diagrams, 
and the next-to-leading order vertex corrections, penguin and hard spectator-scattering. 
It can be written as 



= <^ai(MiM 2 ) A([q s u][uD]) + S pu a 2 (M 1 M 2 ) A([q s D][uu]) 
+ a§(MiM 2 ) J] A([q s D][qq}) + a\{M x M 2 ) Y A([q s q][qD}) 



q 



q 



+ a 



3,EW 



(MiM 2 



^ 2 

q 



e q A([q s D][qq}) + ag EW (AfiAf 2 ) 



J2^e q A([q s q)[qD}), 



(2.44) 



where the operators A([- ••][••• ]) and their matrix elements are the same as introduced 
above. Notice that this is basically the same formula (|2.2ip . The sums extend over 
q = u,d,s and q s is the spectator antiquark. The cf- coefficients are given by 



ai(MiM 2 ) 
a 2 (M 1 M 2 ) 

o|(MiM 2 ) 
c?^M x M 2 ) 

aV{M x M 2 ) 



ai(MiM 2 ), 
a 2 (MiM 2 ), 

f a p 3 (MiM 2 ) - af {M X M 2 ) ; 
\ a|(MiM 2 ) + af(MiM 2 ); 

f ^(M!M 2 ) + a%{MxM 2 ) ; 
1 ^(MiM 2 ) - rf ^ ag(MiM 2 ) ; 



if MiM 2 = PP, VP, 
if MiM 2 = PV, VV, 

if MiM 2 = PP, PV, 
if MiM 2 = VP, VV , 



a 



p 

4,EW 



{M Y M 2 



al(M x M 2 ) - o?(MiM 2 ) ; if MiM 2 = PP, VP , 
al{M 1 M 2 ) + 4(MiM 2 ) ; if M X M 2 = PV, VV , 

< (MiM 2 ) + 2 a p 8 (MiM 2 ) ; if M X M 2 = PP, PV , 
a? (Mi M 2 ) - r^ 2 a\{M x M 2 ); if M x M 2 = VP, V V . 



(2.45) 



where, at next-to- leading order in a s , the coefficients a\ have the following general form 

Ci±i 



<(MiM 2 



+ 



C< + 

JV, 

Ci±i Cfq. 



N c 4vr 



Ni(M 2 ) 

4vr 2 

V(M 2 ) + — fli(MiM 2 ) 



Pf (M 2 ) . (2.46) 



The coefficient Ni{M 2 ) is the leading order term, Vj(M 2 ) contains the vertex terms, the Pf 
are the penguin terms and the Hi(M\M 2 ) are the hard spectator terms. The expressions 
in terms of meson distribution amplitudes and the various hard-scattering functions can 
be found in (23EZJ. 
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The transition operator 7g contains the contributions from weak annihilation and is 
written as 

1% = 5 pu b 1 (M 1 M 2 ) J2B([uq'][<?u}[Db]) + S pu b 2 (M 1 M 2 ) Y, B (WW D][ub\) 

q> q> 

+6g(MiM 2 ) J2 B (WWD][Qb]) + mMiM 2 ) J2 B (WW<l\[Db}) 

<?,<?' 9,9' 

+ 6^ EW (M 1 M 2 ) \ H B([qq'][q'D][qb]) + ^(MM) £ \ e q B([qq'][q ' q][Db]) 

q,q' q,q' 
+ 5 pu b S i(M 1 M 2 ) YB([uu}[q , q'}[Db])+5 pu b S2 (M 1 M 2 ) ^ B([uD][q'q')[ub)) 

q> q> 

+ ^ 3 (M 1 M 2 ) J2B([qD}Wq'}m) + ^S4(MiM 2 ) ^B([?g][gV][D6]) 
q,q' q,q' 

+ ^ 3 EW (M 1 M 2 ) ^ | e q B([qD][q'q'][qb]) + ^^(MiMa) £ | e g B([gg][^][D6]) , 

9,9' 9,9' 

(2.47) 

where the matrix elements of the operators ••][•••][••• ]) are given by 

(MiM 2 |S([. ..][.. .]][. . .])|B g > = cB J¥f2 , with B Mi m 2 =±i ( j= SbJmJm* , (2.48) 

whenever the quark flavors of the three brackets match those of M±, M 2 , and B q . The 
constant c is the same as for the A operators. The upper sign applies when both mesons 
are pseudoscalar and the lower sign otherwise. 

The coefficients with the subscript 'S' contribute only to final states containing flavor- 
singlet mesons or neutral vector mesons. The expressions for the bi coefficients can be 
found in [26|[27]. 

Up to this point it useful to look closely at the contributions that are dangerous in 
QCDF, mainly the hard-spectator and the annihilation terms with twist-3 LCDA insertion. 
The hard-spectator terms receive a contribution of this form: 

H { M lM2 ) D f\ fdy r ^M 2 ^_rnAy) x f {2M) 

A Ml M 2 *B JO JO X V Jo V 

where \b is defined in (|2.42p . This contribution is formally suppressed in the heavy quark 
limit, but the presence of the chiral factor r^ 1 makes it phenomenologically important. 
Moreover, as can be seen by plugging in <3? p (y) = 1 or $ v (y) as given in (|2.38j) . there is a 
divergence coming from the region y ~ 1. This infrared divergence is a manifestation of 
the fact that factorization does not hold in this case. 
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It is customary to extract this divergence defining a parameter X H X = J Q dy/y such 
that 

C (y) = * mi (l)X A H h + f 1 ^ [$ mi (y) - $ mi (1)] (2.50) 

jo y Jo y 

Because this divergence is associated with a soft interaction of the ejected meson with 
the spectator quark, the divergence arises specifically from the region y ~ AQco/mb, 
and therefore one expects that ~ log(m;,/ Aqqd) with an a priori arbitrary strong 
phase. The choice for the values of Xff introduces an unavoidable model dependence in 
the predictions. 

Consider now the chirally enhanced contributions to the annihilation terms with 3- 
twist LCDA insertions. For illustration we take the coefficient b\ in Eq. (|2.47p . which 
contains a term 

6i(MiM 2 ) D ^>Cnra s f dx f dy r^r^^ m2 (x)^ mi (y)^- . (2.51) 

iV Jo Jo x y 

This contribution is again infrared divergent from the regions y,x ~ AQCD/ m b ■ As 
before, the divergencies are parameterized by the quantities X^ 1 and X^ 2 that introduce 
an unavoidable model dependence. The rest of the coefficients in Eq. (|2.47|) are also affected 
by these divergencies, and signal the breakdown of factorization. 

A way to get around these inconveniences at the phenomenological level was proposed 
in [35] and used later in [361 133 ESI US] • The basic idea is that the so called GIM penguin is 
in many cases free from these divergencies when computed in QCDF, since it is dominated 
by short-distance physics. This quantity (that we call A), is given generically by 

A MlM2 = A MlM2 ^C x [G(ml/m 2 b ) - G(0)] (2.52) 

where G(x) are penguin functions. Using this infrared safe theoretical quantity together 
with data and some SU(3) relations, predictions can be made safer and more precise. This 
will be discussed in Chapters [5] and [6) 

2.2 Symmetries 

It should be clear by now that the strong interaction dynamics inside the matrix elements 
is a very big deal. An impressive progress has been observed in the past decade concern- 
ing the idea of factorization in B physics, but direct computations of amplitudes from 
QCD still suffer from two diseases. First, the computations involve a great amount of 
work, combining tedious perturbative calculations, lattice simulations and other techni- 
calities. Second, the theoretical validity of such approaches at the phenomenological level 
is threatened by corrections to the heavy quark limit and by potentially non-factorizable 
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long-distance leading effects, like the case of the long-distance charm loop issue ("charming 
penguins" ) 0QI HQ M H21 SH • 

Fortunately there is a complementary approach to B-decay phenomenology that relies 
on the approximate flavor symmetries of QCD, and incorporates naturally all the hadronic 
information that cannot be extracted from first principles. The drawback is of course the 
uncontrollable symmetry breaking effects, and the fact that only certain relations between 
observables can be obtained. However, these relations have been very useful in order to 
extract CKM parameters from experimental data on non-leptonic decays. In this section 
we review schematically how this flavor symmetry relations are obtained. 

Consider the lagrangian of massless QCD with nf quark flavors 

n f 

Cqcd = — A F a ^F« v + i q k Tp q k (2.53) 

fc=i 

This lagrangian has global symmetry SU(rif)v® SU (n f) a^U (l)v U (1) a- For example, 
for nf = 2 (u & d quarks), this symmetry is expected to hold with good precision due to 
the smallness of m u and m^. The first term SU(2)y is isospin, and it is seen to be a good 
symmetry. The second term SU(2)a is not observed as a symmetry of nature, because this 
symmetry is spontaneously broken by the quark condensate {qq), giving rise to goldstone 
bosons that we observe as pions. The third term U(l)y is nothing but baryon number 
conservation, which is obviously observed. The last term U(1)a led originally to the so 
called [/(l)-problem |45| . but later it was found that in fact this symmetry is broken 
by instantons [16]. So it is clear that, at least for two flavors, the approximate global 
symmetry is well understood and under control. In the case of isospin, the symmetry 
breaking effects are expected at the level of a few percent. 

For nf = 3 flavors, the flavor SU(3) symmetry is not so exact as isospin, since the 
strange quark mass is high compared with m u ^- However, the main feature of a SU(3) 
flavor symmetry is phenomenologically realized: we observe the (tt, K, rf) pseudoscalar 
octet, and the power of group theory can be applied to SU(3) in the same way as for isospin. 
The predictions, nonetheless, will be affected by symmetries breaking uncertainties of 
~ 20%, or maybe higher |47| . 

2.2.1 The Wigner-Eckart Theorem 

The Wigner-Eckart theorem is the mathematical prescription to extract the maximal ad- 
vantage from a symmetry group in a quantum mechanical system. Although the Wigner- 
Eckart theorem can be formulated with total generality, we just quote here the theorem for 
SU{2), as used in elementary quantum mechanics when adding angular momenta. In the 
next section we will use this theorem to write a simple flavor symmetry relation between 
B — ► 7r.fr amplitudes. 
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Given a symmetry group of the system under consideration, any operator can in general 
be decomposed in terms of irreducible tensor operators, denned as operators that transform 
as tensors of a given rank under the symmetry group. In the case of SU(2), an irreducible 
rank-/c operator is a (2k + l)-tuple of operators Q q (q = —k, —k + 1, . . . , k), satisfying 
the proper commutation relations, 

[J±,QW] = ±y/(k^q)(k±q+l) Q% 

[j z ,Q ( q k) ] = qQ { q k) (2.54) 

where J± are the ladder operators and J z is the Cartan operator of SU(2). The Wigner- 
Eckart theorem states that the matrix element (a'; f , m'\Q q \a; j, m) depends on the 
quantum numbers m, m! and q only through Clebsch-Gordan coefficients, specifically 

(a';j',m'\Q q V\a;j,m) = C(j' ,m'\j,k;m,q) (a';j'\\Q k \\a;j) , (2.55) 

where C(f ,m'\j,k;m,q) is the corresponding Clebsch-Gordan coefficient and the object 
( Q 'i j'llQ^II ; j) is a so called reduced matrix element. 

The power of this result is that, although the reduced matrix element cannot be com- 
puted from symmetry principles alone, it only depends on the representation on which the 
objets live in, and therefore one can establish relations between transitions that connect 
states in the same representations. 



2.2.2 A simple example: B — > tiK Amplitudes 

As a simple example of how the Wigner-Eckart theorem applies to flavor symmetries, we 
derive the famous B — > ttK amplitude relation |48^ |4^| [50] . 

We want to extract the isospin structure of the matrix elements 

(7rK\H cS \B) , (2.56) 

where TC e {[ is the weak effective Hamiltonian, given in Eqs. (ll.27| ). (jl.28p for the SM case. 
First we decompose the initial and final states in the isospin basis. (B°, B + ) and (K°, K + ) 
are isospin doublets, and (jr~ -,7r°,7r+ ) is an isospin triplet: 

\B°) = |l/2,-l/2) , \B+) = |l/2,l/2) , 
\K°) = |l/2,-l/2) , \K+) = |l/2,l/2) , 
IO = |1,-1) , K°) = |1,0) , |tt+) = |1,1) (2.57) 
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The isospin decomposition of the I7r.fr) states is given by the (1 <g) 1/2) SU(2) Clebsch- 
Gordan coefficients, 



\tt + K°) 



3 1 

2 '2 

3 1 

2 ' 2 



+ 



1 1 

2 ' 2 

1 1 

2 ' 2 



|vr"K" 



k°K°) 
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1 1 

2'~2 
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(2.58) 



The next step is to decompose the effective Hamiltonian in terms of irreducible SU(2) 
tensor operators. But this is easy: the current-current operators Q\ 2 ~ ubsu have com- 
ponents AI = ad AI = 1, and so do the electroweak penguin operators Q7,...,io ~ 
Db eqQQ', the current-current operators Q\ 2 ~ cbsc and the QCD-penguin operators 
Q3,...,6 ~ sb^qq only have AI = components. Therefore, the effective Hamiltonian can 
be decomposed as the sum of triplet and singlet tensor operators: 
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n 



(0) 
eff 



(2.59) 



Now the matrix elements can be expressed in terms of reduced matrix elements. First 
have (3/2,±l/2|wJ5|l/2,±l/2) = due tc 
rest, the Wigner-Eckart theorem implies that 



we have (3/2, ±l/2|ft^ 1 1/2, ±1/2) = due to the commutation relations (12331) . For the 
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(2.60) 



It is conventional to define the invariant isospin amplitudes A3/2, Ax/2 and Bx/2 , 
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^1/2 = 
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2 /3 
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so that combining Eqs. (|2.58p and (12.60P the amplitudes can be written as 

A(B + ^tt + K°) = B 1/2 + A 1/2 + A 3/2 , 
B 1 / 2 + A1/2 - 2A 3 / 2 , 



-V2A(B + -> 7r°ir + ) 



A(B° = B 1/2 -A 1/2 -A 3/2 , 
= B1/2 - Ax/2 + 2A 3 / 2 



(2.61) 
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From this representation of the amplitudes in terms of reduced matrix elements, the 
promised B — > ttK amplitude relation follows: 

A{B + -> 7T + K°) - A(B° ^~if + ) + V2A(5+ -> vr°K + ) - ^2^(5° vr°K ) = (2.62) 
2.2.3 Discussion 

Following the same procedure as for B — > ttK, it is possible to write other isospin relations. 
For example, for B — > tttt , isospin implies 

A(B° -> tt+tt^) - V2A(B + vr+vr ) + -> vrV) = . (2.63) 

Relations arising from SU(3) flavor symmetry are richer, since they mix many other decays 
that are not just connected by an isospin transformation. The complexity of the group 
SU(3) with respect to the simplicity of SU(2) is the main reason to avoid a detailed 
discussion of the technical aspects of the SU(3) analysis here, and the important idea is 
clear from the example above. The SU(3) Clebsch-Gordan coefficients were collected for 
the representations of interest by de Swart [51], and the complete SU(3) decomposition 
of the non-leptonic decay amplitudes of the triplet B®,B®) into pseudoscalar mesons 
was given by Grinstein and Lebed |52j . 

However, the general decomposition of SU(3) amplitudes introduces too many SU(3) 
invariant amplitudes, so that in order to extract amplitude relations it is often necessary 
to neglect some of them. But the relative sizes of the reduced matrix elements is a dy- 
namical question, so this process requires a dynamical assumption. Because the language 
of invariant amplitudes is often obscure in this sense, it is often convenient to use a dia- 
grammatic approach in which the invariant amplitudes are related to flavor-flow topologies 
(see for example [531). The amplitudes can then be written in terms of T (tree), P (pen- 
guin), C (color suppressed), E (exchange), A (annihilation), PA (penguin annihilation) 
and electroweak penguins Pew an d Pew Including b — > s transitions doubles the num- 
ber of parameters, which are usually denoted with primes. Then one can argue that the 
electroweak penguins are small, and also neglect the contributions E, A and PA on the 
basis of helicity or is/m^-type suppressions [53] . It is, however, far from clear whether 
these assumptions can be made safely. 

An analysis along these lines allows to find flavor symmetry relations that do not rely 
on any dynamical assumption, like for example |55[ [56] 

A(B° K°7T ) - V3A{B° -» K°7]s) = , 

A(B + K+7T ) + V2A(B+ -> K°tt + ) - V3A(B + -> K+rj 8 ) = , (2.64) 

A(B° -> K+7T-) + A{B + -> K tt+) - V6A(B+ -> K + V8 ) + \/6A(B° -> K%) = . 

From these relations the isospin B — ► irK relation f|2.62)) follows. 
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Up to this point it not obvious what is the usefulness of these amplitude relations. 
Indeed, a powerful aspect of flavor symmetry relations is to be able to write relations 
between observables and not just amplitudes. These relations can be then tested exper- 
imentally to probe the size of the SU(3) breaking. Since they are specially sensitive to, 
for instance, non standard isospin violating contributions, they provide a window for New 
Physics (see for example |57| for a discussion concerning B — > ttK). 

But obtaining relations between observables is not as easy as for the amplitudes. Ob- 
servables are constructed, roughly, by squaring amplitudes. Therefore, the observables 
depend on the magnitudes of the invariant amplitudes and also on the strong phases, since 
there is interference between different invariant amplitudes. For example, in the case of 
Eq. (|2,6ip . there are two independent CKM parameters entering the decay amplitudes, 
introducing two weak phases. Each invariant amplitude receives contributions from both 
weak phases. Up to a global phase, one ends up with five independent strong phases 
and six independent real amplitudes, to be fitted with four branching ratios and four CP 
asymmetries [58J. This example shows how a model independent amplitude relation is not 
sufficient to obtain a model independent observable relation. 

Despite the apparent negativity of the discussion above, flavor symmetries have been 
used in phenomenology with the outcome of very interesting analyses and results, helping 
to achieve an understanding on B decays that would have not been possible with dynamical 
approaches alone [5SII5SIEQIEIIEZIIS31I531E51ES!IE 

Flavor symmetry relations will be 
used extensively throughout this thesis, specially for B — > KK modes. 



2 To name a few. 



Chapter 3 

CP Violation in Meson Decays 



The fact that nature is not invariant under parity transformations (P) was an early sur- 
prise. It was of general common sense to think that one could never tell whether a physical 
process was being observed through a mirror or directly. Indeed, electromagnetism and 
strong interactions conserve parity, and in the 1950's experiments in particle physics were 
mainly probing the recently developed theory of QED, and the strong interactions of the 
zoo of hadrons that was just being discovered. The exception was the experiments on 
/3-decay, that was at that time described by Fermi's theory. This theory had just been 
proposed to be universal (e.g, neutron decay and muon decay would have the same cou- 
pling constant), and this was the main issue in weak interactions. 

However, in 1955 the r-6 puzzle raised some confusion. The 9 + would decay in to a 
CP-even two pion state, and the r + into a CP-odd three pion state. Parity invariance 
necessarily implied that these were two different particles with opposite parity, but the 
puzzling point was that they were otherwise identical] same charge, same mass, same 
width. This led to Lee and Yang to point out that while P invariance had been tested in 
strong and electromagnetic interactions, it had never been verified in weak interactions, 
and they proposed several experiments. This proposal was taken seriously by Wu's group, 
who in 1957 provided first evidence of Parity violation studying the angular distribution 
of nuclear /3-decay with polarized Co 60 [67] . Their discovery was simultaneously confirmed 
by Lederman and collaborators with a measurement of the electron asymmetry in pion 
decay [68] • The discovery of P violation immediately solved the t-6 puzzle, a crucial brick 
in the physics of kaons. 

In parallel, the theory of strangeness arose as an explanation to the kaon-hyperon 
problem. These particles were produced very easily in pairs, but they decayed very slowly. 
Gell-Mann and Nishijima proposed the existence of a quantum number (strangeness) that 
was only violated by weak interactions; pair production could then be a strong process, 
but single decay had to be weak. This was the foundation of flavor physics. An immediate 
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consequence was the existence of two different neutral kaons with opposite strangeness. 
The separate identity of these two kaons was established very soon by Gell-Mann and Pais 
in 1955 [69], a discovery that gave birth to the physics of neutral meson mixing. 

The oscillation of neutral mesons quickly opened up a chapter in the history of par- 
ticle physics. The measurement of the oscillation period of neutral kaons through time 
dependent semileptonic decays allowed the extraction of the tiny mass splitting between 
the two physical kaons, AMk/Mk ~ 1CP 15 . Later, studies of regeneration of short-lived 
kaons in matter revealed an excess of regenerated kaons, a result that had an unexpected 
origin. In 1964 an experiment by Christenson and collaborators |70] established that both 
the short-lived and the long-lived components of the neutral kaon system were decaying 
into a CP even ir + Tr~ mode. The clear interpretation of this discovery was the violation 
of CP invariance. 

Today we have much evidence of CP violating phenomena in kaon and Bd meson decays 
|71j . Measurements of CP asymmetries in B s decays are still compatible with zero at the 
2a level, but since the experimental prospects are good at Tevatron, and theoretically 
these asymmetries can be as large as for decays, it is safe to say that it will not take 
long. CP violation has neither been observed in D meson decays, which is very suppressed 
in the SM. However, at the LHC the sensitivity could reach the 1 per mil for D meson CP 
asymmetries, making CP violation in the D sector a powerful probe of New Physics. 

The SM description of CP violation, the CKM mechanism [72], is successful and eco- 
nomic, involving just one CP violating phase. The fits to the Unitarity Triangle are getting 
tight [73(. 174]. and all data is consistent (within experimental and theoretical errors) with 
a single value for its apex. However, the main reason to believe that there must be sources 
of CP violation beyond the SM, is the issue of the baryon asymmetry in the universe. 

The measurement of baryon asymmetry has recently become accurate (about 5% error) 
due to the measurements of the fluctuations of the cosmic microwave radiation background. 
The order of magnitude of the baryon-to-photon ratio is (jib — n B)/ n -/ ~ 10~ 10 . The fact 
that this asymmetry could have been an initial condition at the big bang is excluded, since 
any baryon asymmetry present in the early universe would have been diluted by inflation 
to completely negligible levels. So it is very much accepted that the baryon asymmetry 
must have a dynamical origin. 

The necessary conditions for the dynamical generation of a baryon asymmetry in the 
universe were studied by Andrei Sakharov in 1967 [75]. There are three main conditions 
that have to be fulfilled. First, there must exist baryon number violating interactions (this 
is clear). Second, there must be a departure from thermal equilibrium during a sufficiently 
long period of time. This is because the baryon number violating interactions must occur 
in the forward direction predominantly, so there must be an arrow of time. Third, there 
must be CP violation, since CP invariance would equal the rates of CP-conjugated baryon 
number violating interactions, thus giving a zero net contribution. 

Now, baryon number is broken in the SM by sphalerons at high temperature. The 
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electroweak phase transition, occurring at a temperature T ~ 100 GeV, could provide the 
necessary departure from thermal equilibrium (this is called electroweak baryogenesis). It 
turns out that this phase transition is not strong enough to produce and maintain the 
observed baryon asymmetry, but let's assume it were. Then, the baryon asymmetry can 
be estimated to be (tib — n B)/ n -i ~ 10 -2 • J/T 12 , where J is the Jarlskog determinant 
(a measure of CP violation in the SM). For T ~ 100 GeV, then (rig — n^)/n 7 ~ 10~ 23 , 
so clearly the CP violation in the SM is far too small to account for the observed baryon 
asymmetry. On the other hand, almost any extension of the SM provides new sources of 
CP violation, so the search for physics beyond the SM is much in contact with the study 
of CP violation. 

In general, CP violation is the cause that a process and it's CP-conjugate don't share 
the same rights. This means that it is possible to distinguish objectively between CP- 
conjugated events^. In particular, a difference between the decay rate of a given process, 
A — > B, and that of the CP transformed, A — > B, will signal CP violation. This is 
generally quantified by a so called CP asymmetry, defined as 



The deviation of any CP asymmetry from zero is a measure of CP violation. In this 
chapter we will review the general features of CP violation in meson decays and the 
different types of CP asymmetries that one may encounter. 

3.1 Mixing of neutral mesons 

The phenomenon of neutral meson mixing is a consequence of flavor violation. Therefore, 
the mixing is a weak process, and since the mesons themselves are regarded as asymptotic 
states of the strong Hamiltonian, it is sensible to study the mixing in terms of mesonic 
degrees of freedom. 

The quantum mechanical idea for this phenomenon is very simple. The weak Hamil- 
tonian is not invariant under flavor rotations. This means that flavor eigenstates deviate 
from mass eigenstates and mix under time evolution. This deviation, as already men- 
tioned, is due to weak effects, and this allows to study the oscillations in a perturbative 
fashion. 

Experimentally, oscillations have been observed in all the low-lying neutral meson 
systems. The oscillation of neutral kaons was observed for the first time in 1961 [76]. The 
oscillations of Bj mesons was observed in 1987 \77\ 178] . The observation of B s and D 





1 For example, it is possible to define the electric charge in an absolute sense: positive charge is the 
charge of the lepton more often produced in the semileptonic Kl decays. 
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meson oscillations has taken a much longer time, but they have been observed finally, very 
recently [751150]. 

In this section we present a general description of particle-antiparticle oscillations in 
a quantum mechanical approach. We will see in later chapters that the physics of meson 
oscillations is strongly related to the physics of CP violation. 

3.1.1 The effective Hamiltonian 

Let P° be a neutral meson carrying an internal additive quantum number F^O. For our 
purposes F is a flavor quantum number. Charge conjugation will change the sign of this 
quantum number, so the C-conjugate state P° is a different neutral meson. 

Now, let's say that F is not conserved, that is, the interaction Hamiltonian contains 
(weak) breaking terms with AF ^ 0. Transitions P° <-> P° are then possible as a one-step 
process (through Hfyf =2 ) or through an intermediate state (at second order in H^f =1 ). 
To be specific, we write the full Hamiltonian as 

H = H + H%f= 1 + H%f= 2 , (3.2) 

where Hq is the flavor invariant (strong) Hamiltonian. The states \P°) and \P°) are then 
defined as eigenstates of Hq. As mentioned above, that it actually makes sense to speak 
about P° and P° is due to the fact that H\y is weak. On the other hand, the presence of 
the terms H\y will produce P° — P° mixing. 

Consider now an initial state which is a mixture of a P° and a P°, 

\il>(t = 0)) =c\P°) + c\P°). (3.3) 

At t > this state will evolve in two different ways. First, transitions between P° and P° 
will give rise to unitary rotations in the c — c space (oscillations). Second, the fact that 
these mesons can decay into lighter particles will allow the initial state to evolve outside 
the two meson system: 

\m) = c(t)\P ) + c(t)\P°)+J2cn(t)\n) (3.4) 

n 

where \n) = | 7T7t) , |37r), \ttWi), ... represents any state of any number of particles which does 
not violate symmetry requirement -mainly decay modes of the original mesons. We would 
like to study the time evolution of a general initial state (|3.3[) inside the 2 dimensional 
P° — P° subspace. To that end we split the full Hamiltonian into a part Hq = H&f=o 
and a part H\y = H/^f^q, H = Hq + H\y, and we follow a time-dependent perturbation 
theory formalism treating H\y as a small correction. 

In the interaction picture, the state vectors and operators are defined as follows 

\m)i = j m m))s (3.5) 

T (t) = e iHot O s e- iHot (3.6) 
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The time evolution in this picture defines the evolution operator Ui(t) through 

\m)i = uiimm, (3.7) 

so that the evolution of the coefficients c(t) and c(t) is given exactly by 

c(t) = (P \Ujm(0)) =c(P°\U I (t)\P°) + c(J*|l7 7 (t)|P°) . . 

c(t) = {popjftMO)) = c(P \U I (t)\P )+c(P°\U I (t)\P ) W 

The unitary operator Ui(t) verifies the Schrodinger equation 

id t U T {t) = H w (t) • 17/(4), (3.9) 

where Hw(t) = (Hy/)i- The solution of (|3.9|) with initial condition Ui(0) = 1 can be 
computed perturbatively 

C7/(i) = 1 — z / dt'H w {t')- [ dt' I dt'H w {t')H w {t") + ••• (3.10) 

We are interested only in the time evolution of the projection of the full state \ip{t)) 
on the two-dimensional subspace of the P° — P° system. This evolution is NOT unitary, 
and effectively is described by a Schrddinger-like equation: 

ldt [mj ~ cS \m) ~ V (P°\n eS \p°) {P°\n c{i \P ) J Km) [ ' 

This is an effective description valid at time scales much larger than the typical strong 
interaction scale, and it's called the Wigner-Weisskopf approximation [81j . In particular, 
the effective Hamiltonian H. e ff will not be hermitian. 

By introducing (13. 8p in (|3.1ip . and using (13. 101) . we find an equation for Ti. c s- 

(a\H w (t)\b) -i(a\H w (t) [ dt' H w (t')\b) + ■ ■ ■ 

Jo 

= (a\H eS \b) -i(o|VWeff|m)H / dt' H w (t')\b) + ■ ■ ■ (3.12) 

where a,b,m £ {P°,P }. We can solve this equation order by order in H\y. For our 
purpose here it will be enough to stop at second order: 

He$ = H w + H w \n)(n\H w (v ^ - - inS(m P - E n )) +■■■ (3.13) 

^-f \ m P - E n + le J 
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Here, V stands for the principal part prescription, and the sum runs over all intermediate 
states n. mp and E n are the energies in the center of mass frame defined as Hq\P°) = 
m P \P°), H \P°) = m P \P°) and H \n) = E n \n). 

Up to this point we introduce two hermitian operators M and T defined as 

M = \{H eS + H\ s ) 

r = i(n cS -nl s ). (3.14) 

By inversion of (|3.14p we can write the effective Hamiltonian in terms of M and T, 

1 ( Mn - |r n Mi 2 - |ri2 \ 

HeS = M r = 11 < 3.15 

2 ^ m 2 i - |r 2 i m 22 - §r 22 J 1 ; 

and the explicit matrix elements for M and V can be read off from equations f)3. 13j) and 

in, 



Mn = mp + ^V- 



(n\H^\ P 0)\2 
m P - E n 



0| ITAF^l I „\ /„ I rrAF=l| z30\ 



M 



12 



/ pOl ZJAF=1|„\/„| pjAF=l| pO 



T n = 27r^5(mp-P n )|(n| J ff^= 1 |P°)| 2 
T 22 = 27rY,S(mp-E n )\(n\H^ F = 1 \P )\ 2 

n 

r i2 = r* 1 = 27r^5(mp-P n )(P°|^=»(n|^= 1 |P ) (3.16) 

n 

Assuming CPT as a symmetry of H\y leads to the following relations 

M n = m 22 = M , r n = r 22 = r . (3.17) 



3.1.2 Mass eigenstates: Diagonalizing the effective Hamiltonian 

The mesons P° and P° are charge-conjugated states, and hence they differ only by the 
sign of their internal additive quantum numbers. They are eigenstates of Hq, which 
conserves such quantum numbers, but not of Hy/- Therefore these are not physical states 
(or mass eigenstates) with the corresponding consequence of mixing and decay. There is 
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a feature, however, that decouples the process of mixing with that of decay, and allows 
to treat them separately: oscillations are unitary in the two dimensional subspace, but 
decay is not. A new basis {Pl, Ph} which diagonalizes TC e g defines two fields that do not 
oscillate, just decay. Of course the question arises of whether M and T can be diagonalized 
simultaneously. In general they cannot, so no unitary transformation can diagonalize TC c s, 
and the states {Pl, Ph} will not be orthogonal. The formalism to deal with this issue was 
developed by T.D. Lee and L. Wolfenstein in 1965 [82] . 

To simplify the notation, we parameterize the effective Hamiltonian in terms of three 
complex numbers, A, B and r: 

«- = ( rB B 'a ) <» 8 > 

where CPT has been assumed through (|3.17[) . By comparing (|3.18|) with (|3.15p we see 
that El 

A = M -^r (3.19) 



B = \j(M 12 - l -T 12 ){M* l2 - l -T\ 2 ) (3.20) 



m* — ir* 

1V1 12 2 L 12 



V Mi 2 - §r 12 

The eigenvalues and eigenstates are then 



(3.21) 



H cS \P L ) = {A + B)\P L ) = {M L - % -T L )\P L ) 



2 

H e s\P H ) = (A - B)\P H ) = {M H - ^h)\Ph) (3.22) 



\P L ) = / __ T (|P°)+r|P°)) 



\T 



P H ) = / __ (|P°)-r|P°)) (3.23) 



As mentioned above, the mass eigenstates {Pl,Ph} are n °t necessarily orthogonal. 
This can be easily seen by inspection of (|3.23p . since the non-hermiticity of Ti. e g allows 



2 In the standard literature (see for example |83l 184] ). what we call r is denoted by (q/p), with 



1 r 

P = n , „ - 1 = 
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\r\ ^ 1. However, the amount of non-orthogonality is constrained by the exclusive widths 
of Pl and Pjj- The proof goes as follows. Consider the braket 

(Ph\T\Pl) = (P H \i(H eB - Ht s )\P L ) = [i(T £ + T H ) + i(M L - M H )] (Ph\Pl) (3.24) 

This quantity can be calculated from (|3.16p . 

(P H \T\P L ) = 27r^5(7np - E f )(f\H\P H )*(f\H\P L ) . (3.25) 
/ 

We denote (f\H\P) = T*. Then by means of the Schwartz inequality, |(Pif|r|Px,)| < 
Y, f r{r{f, we find that 



The equation (|3.26p is called the Bell- Steinberg er inequality. Note that if H^f =1 = 0, 
then the states do not decay, only mix (if H§f= 2 + 0, of course). Therefore = 0, but 
still Ml ^ Mjj, so (Ph\Pl) = 0. This would correspond to a conventional unitary mixing. 

3.1.3 Time evolution 

We denote by \P(t)) the state at time t that at t = was a pure |P). The time evolution 
of the flavor states P°,P° is complicated because of the mixing, in the sense that the 
states \P°(t)) and \P°(t)) will be superpositions of P° and P° at t > 0. However, the time 
evolution of the mass eigenstates is quite simple since they diagonalize the effective Hamil- 
tonian. From (|3.1ip and (|3.22j) we find the usual time evolution for decaying stationary 
states: 

\P L ,H(t)) = e-^.HV-K^lPi.tf) (3.27) 
By inversion of (|3.23|) we can find the time evolution of the states P°,P°, 

|P°(i)> = g + (t)\P°) + rg^t)\P ) 

|P°(i)) = g + (t)\P°) + -g^t)\P ) (3.28) 

r 

where, to avoid cumbersome expressions, we have defined 
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an 



AM = M H - M L , AT = T l - T h . 



(3.30) 



We are interested in the time evolution of observables such as branching ratios, decay 
rates, and CP asymmetries. We define Af and Af as the amplitudes for the decay of P° 
and P° into a final state / , i.e., 



A f = (f\H^\P») , A f = (f\H^\P^ 



(3.31) 



The decay rates are proportional to the square of the time dependent decay amplitudes, 
the proportionality factor given by a phase space factor / PS (see the Appendix at the end 
of this chapter). From (|3.98[) . (|3.28[) and (|3.3ip we find the following master equations for 
the time dependent decay rates, 



r(P°(i)->/) = i / PS e~ TLt \A f \ 2 [K+(t) + K-(t) \X f \ 2 + 2Re(L(t) A/) 



1 



r(P (t)^f) = -f PS e-^ t \A f \ 2 



K + (t)\\ f \ 2 + K_(t) + 2Re(L(i)AV 



(3.32) 
(3.33) 



where 



and 



K±(t) = 1 + e Ar * ±2e5 Ar * cos AMt 
L(t) = 1 - e Art + 2ie^ Art sin AMt 

1 A f 



(3.34) 
(3.35) 



3.1.4 CP Violation in the neutral meson system 

Up to now all the discussion about the neutral meson system did not make any reference 
to CP symmetry. All the previous results are independent of whether CP is or isn't a 
good symmetry of the full Hamiltonian. However, it turns out that some of the quantities 
defined so far are closely related to the amount of CP violation. Therefore, the study of 
CP violation is very useful when extracting physical content to the formalism above. On 
the other hand, the study of such systems provides a lot of information about the nature of 
CP violation, to the extent that CP violation can be quantified even from CP conserving 
processes. 



3 It's not clear yet from the construction if AM and AF are positive or negative. It is a matter of the 
choice made in (|3.22|l as to which mass eigenstate should we call \Ph) or \Pl)- We can take by convention 
to call the heaviest one with the subscript "H" - not necessarily as in (|3.22[) - so that AAf > 0. But once 
this convention is taken, it becomes an empirical question whether Ar is positive or negative. 
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The role of CP-violation in the mixing of neutral mesons can be identified by answering 
the following question: How are the results above modified (or simplified) if we impose 
CP as a symmetry? The answer is twofold: 

• We must identify the mass eigenstates (I3.23P with the CP eigenstates. The action 
of CP on the flavor eigenstates i£| 

CP\P°) = e^ p \P°) , CP\P°) = e-^ p \P°) (3.36) 
so the CP eigenstates (defined as CP\P±) = ±P±) are 

|P + ) = -L(|pO) +e ^|pO)) 

|P_) = -L(|P°)- e ^|P°» (3.37) 

Comparing (|3.23p with (|3.37[) we see that in the limit of CP invariance, r = e l ^ p . 
This is a convention dependent quantity, so it is not quite interesting. However, its 
modulus 

\r\ ^ 1 (3.38) 

it's an observable, and its deviation from 1 measures CP violation in the mixing of 
neutral mesons. 

• The amplitudes Aj and Aj defined in (|3.3ip are equal under CP invariance up to 
an unphysical phase: 

Aj = (f\H\P°) = (f\CP ] (CPHCP ] )CP\P°) ^ (/|e if /P" e -^|P°) = e^l'^Af 

(3.39) 

where the phase £f is the analogue to (|3.36|) in the definition of |/). The ratio of 
the moduli of these amplitudes is independent of these phase conventions and its 
difference from unity is a measure of CP violation in the decay process, 

\A f /A f \^l (3.40) 

In the absence of mixing, as in the case of charged meson decays, this is the only 
source of CP violation. 

If /cp is a CP eigenstate, i.e. CP|/cp) = f?/|/cp) with rjf = ±1, then 

%p = (fcp\H\P°) = r) f e^f(fcp\H\P°) = n f e^ A fcp (3.41) 



4 The phase £p is convention dependent and hence unphysical. Because the convention in the election of 
such phase is related to other conventions in the definition of other quantities used in the text, it is useful 
to keep track of it for the moment. 
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and comparing (|3.39j) and (|3.4ip we find that 

This ratio is convention dependent, as contains the phase £p. But this dependence 
is cancelled by r inside the quantity Xj defined in (|3.35p 

A/cp ^/ (3.43) 

This is a physical quantity, whose modulus and phase (or equivalently its modulus 
and its imaginary part) signal CP violation in the decay and mixing processes. 

The results (13.38p . (|3.40p and (|3.43p contain most of the qualitative information on 
how CP is violated in meson decays. 

3.2 Classification of CP-violating effects 

Now that we have seen how CP manifests itself in physical quantities, we are ready to 
make a classification of CP-violating effects. Fortunately, the neutral meson system is rich 
in CP-violating phenomenology, so mostly everything was said and done in section 13.1.41 
There are three main independent ways in which CP is broken in meson decays: 

* CP Violation in decay: \Aj/Af\ ^ 1 
Consider the situation in which no oscillations occur, i.e., 

B = or AM = AT = (3.44) 
Then, following from the master equations (|3.32p and (|3.33|) . we see that 

r(P-/) = f M e- T ^\A f \ 2 

r(P-/) = f PS e- rLt \A f \ 2 (3.45) 

A difference in the CP conjugated decay rates T(P — > f) 7^ T{P — > f) is a signal of 
CP violation, and arises if 

\A f /A f \^l (3.46) 

as advanced in (|3.40p . This type of CP violation is the only possible one in charged meson 
decays, where mixing effects are absent, and hence CP violation in decay is best measured 
in these modes. The relevant asymmetry is then 

A (P± f ^- r(P+ ^ /+) ~ r(P "^ r) - 1 ~ l ^- / ^ +|2 f347^ 

Acp{P n = itF^TW^I " i + l^/-M /+ l 2 (3 7) 
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-k CP Violation in mixing: \r\ ^ 1 

Consider a final state / that can only come from P , not from P°, and also the CP 
analogue, that is 

Af = A f = (3.48) 

These are called flavor-specific decays. Then, the transitions P°(t) — ► / and P°(£) — > / 
are only possible because of mixing; for example P°(t) -> f = P° -A P° -> f. These 
transitions are called "wrong-sign" decays. The wrong-sign decay rates, according to 
(13.321) and (|3.33p . are given by 

r(P°(i)-/) = \f PS e- r ^K.(t)\r\ 2 \A f \ 2 , 



T(P°(t)^f) = \f PS e- T ^K^{t) 



2 

\A f \ 2 . (3.49) 



Consider also that there is no CP violation in decay, such that \Af\ = \Aj\. The 
standard example of this situation is the case of charged-current semileptonic neutral 
meson decays. These are decays P°,P° -tFl, which in the Standard Model verify 
|.(4.;+_x-| = l-^z-xl and Ai~ x = ^i+x = to lowest order in Gp. The CP asymmetry in 
this case is given by 

_ r(p°(t)^f)-r(P°(t)^f) _ i-\r\* 

SL T(pO(t) -» /) + T(P%t) -» /) 1 + |r| 4 1 " ' 



which is nonzero whenever 



|r| + 1 (3.51) 



in agreement with (|3.38p . Note that in this case the CP asymmetry between time- 
dependent decay rates is actually time- independent. These CP asymmetries are suited 
for the extraction of the mixing parameter \r\. 

•k CPV in the interference between mixing and decay: Im(Xf) ^ 

In the case of flavor non-specific decays, that is when a final state / can be reached by 
both P° and P°, two amplitudes interfere in the process: 

A(P°(t) - /) ~ A(P° - /) + A(P° ^P°^f) (3.52) 

Even if neither the decay itself nor the mixing introduce CP violation, the interference 
between these two decay channels can produce a nonzero CP asymmetry. The most 
transparent case arises in the situation in which the final state is a CP eigenstate /cp- In 
the absence of CP violation in mixing and decay, |Aj CP | = 1 holds. Still a deviation from 
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(|3.43p can arise from the phase of A j CP . In order to fix notation we define the mixing angle 
4> M and the decay angle <j) D aM 

c(> M = arg (r*) , 4> D = arg (A f /A f ) (3.53) 

so that 

ar S ( A /cp) = -(4>m+4>d) (3.54) 
In this situation, eqs. (|3.32p and (|3.33p imply 

r(P°(t)^/ C p) = i/p S e- r ^V/cp| 2 f 1 + e Ar * + Re(L(t)e-^ + ^ 



r(P°(t) - /cp) = i/p S e- r ^V/cp| 2 (l + e Art + Re(L(t)e^ + ^))) (3.55) 
The CP asymmetry is then 

r(P°(t) - /cp) - r(P°(t) - /cp) sin ((t> M -H>n) sin AMi 



Ajp(t) 



r(P°(t) - /cp) + r(P°(t) -» /cp) cosh (|ATt) - sinh (±AH) cos (</>m^ d ) 

(3.56) 

This time-dependent CP asymmetry is nonzero only in the presence of P° — P° oscillations 
generating AM ^ and whenever 

Im(A/)/0 (3.57) 

3.3 Strong and Weak phases 

Consider a decay process and its CP conjugate. If CP is not conserved, the two amplitudes 
don't need to be correlated and can be completely different in modulus and phase. Two 
arbitrary complex numbers Af and Aj can always be decomposed in the following way: 



A f = laile^e^ 1 + \a 2 \e i&2 e i4>2 



A f = loile^e - ^ 1 + \a 2 \e i52 e-^ 2 (3.58) 



So the amount of CP violation can be encoded inside the so called CP violating phases 
4>i, which change sign under CP. These phases arise from complex parameters in the 
Lagrangian that contribute to the amplitudes, since the couplings in the Lagrangian appear 
in complex conjugate form in CP conjugate amplitudes (see for instance [El]). In the 
Standard Model complex parameters appear only in the weak sector - through Yukawa 
couplings-, so these CP violating phases are due to weak interaction and therefore are 
usually called weak phases. 



5 This definition for cj> M corresponds to 4>d — 2/3 in B a d — B° d mixing, as will become clear later. 
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The CP conserving phases (<5i — 62 7^ 0) are necessary (together with (f>x — 02 7^ 0) to 
account for a difference in modulus between the two CP conjugate amplitudes. They arise 
even if the Lagrangian is real and generally through a process know as rescattering. This 
is a contribution from a possible final state interaction, with on-shell intermediate states. 
In the Standard Model the main source of CP conserving phases is the strong interaction, 
so these phases are commonly referred to as strong phases. 

Obviously, the weak and strong phases by themselves are not convention independent, 
since they depend (at least) on the choice made for the phases £pj. However, the relative 
strong and weak phases between the different terms in the amplitudes are physical. 

The decomposition (|3.58p is particularly explicit in the SM. As can be seen from the 
AF = 1 SM effective Hamiltonian in eq. (II. 27ft . the AF = 1 amplitudes are the sum 
of a term proportional to A« D) and a term proportional to Ac 0> . This is a consequence 
of CKM unitarity as will become clear in Chapter HI The factors A^c contain all the 
CKM information, so the CP violating phases are all contained in these factors. On the 
other hand, the matrix elements of the operators will contain strong interaction phases, 
as explained in Chapter [2j At the end, the two CP conjugated amplitudes can be written 
in full generality as 

A = \ ( U D) * T + A< D) * P , A = X { U D) T + \ { C D) P , (3.59) 

where T = \T\e it>T and P = \P\e iSp are complex hadronic parameters, called "tree" and 
"penguin" , and 5t,p are strong phases. This representation of the amplitudes is exactly 
of the form of eq. (|3.58|) . We shall go very much into detail with this decomposition in 
the following chapters. 

3.4 CP Violation in B decays 

The discussion about CP violation in the previous sections is quite general and describes 
within the same framework all mesons of different families, such as K 7 s, D's, B's, etc. 
However, the variety of values that nature has chosen for their lifetimes, masses and 
oscillation patterns results in very different phenomenologies. The phenomenology of B 
decays is particularly interesting and contains a number of characteristic features |83| : 

• Nature has provided us with two neutral B meson systems, with a good measure of 
similarities and differences. The mass differences for both systems are [85} [86] 

AM d = 0.507 ± 0.004 ps" 1 (3.60) 
AM S = 17.77 ± 0.12 ps" 1 (3.61) 

• CP violation in mixing can be ignored completely in present analyses of B decays, 
both in the Bj and B s systems. This statement is based on recent measurements of 
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02468 10 12 14 1 2 3 4 

t/T Bd t/T Bs 



Figure 3.1: Ocillation probabilities for Bd and B s mesons. The time axis is given in units of the 
meson lifetimes. Clearly the oscillation of the B s is much faster than the oscillation of the Bd (note 
that they are not drawn to the same scale). 



the CP violating parameter r [85] : 

\r\ d = 1.0033 ±0.0017 (3.62) 
\r\ a = 0.9998 ±0.0046 (3.63) 

• The width difference in the B d system is found to be small. While an initial beam of 
K° and K° mesons is transformed with time into a practically pure Kl beam, this 
does not happen with B d /B d beams in any appreciable way. The width difference 
for the B s system does not share the same fate, and this introduces an interesting 
difference between the two different systems of neutral B mesons. Recent data for 
the lifetime differences are [85]: 

AT d /T d = 0.009 ± 0.037 ± 0.018 (3.64) 
Aiyr, = 0.104+°;™j (3.65) 

After the exposition of Section 13.11 it is easy to understand why it has taken much 
longer to measure AM S , which was just measured last year, than to measure AMj, which 
was measured in the 80's. Of course one of the reasons is that the production of B d 
mesons was possible very early, and the S-factories operating at the T(46 r ) resonance 
have produced millions of B d — B d pairs, while the detailed study of B s mesons has only 
been possible recently at Tevatron. The other reason is that while the B d meson oscillates 
relatively slowly, the oscillations of the B s meson are very fast, oscillating on average 25 
times before decay. Indeed, looking at eq. (|3.28j) . it can be seen that the function g+(t) 
measures the frequency at which a P° or a P° meson turns back into itself. From eq. 
(13.29ft . the probability of finding a P° at time t from an original P° is given by 

V(P°(t) = P°) = \g+(t)\ 2 ~ i e - f *(l±cosAMt) (3.66) 
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Figure 3.2: Amplitude of — B® oscillations as a function of the frequency AM S . These data 
proved that B s mesons oscillate, and measured their mass difference. 



where f is the average lifetime and the approximation is valid if Ar <C T. Then, AM is 
the frequency of the oscillations, so from (|3.60p and (|3.6ip we see that B s mesons oscillate 
about 35 times faster than mesons. Figure I3.ll shows the oscillation probabilities of 
both neutral mesons. Figure I3.2I shows the experimental data revealing the oscillation 
frequency AM S . 



3.4.1 CP Asymmetries for B decays 

The small CP violation in mixing pointed out by (I3.62|) and (|3.63j) tells that r is then a 
pure phase to an excellent approximation (r ~ e~ lr f >M , see (|3.53p ). The time dependent 
CP asymmetry for B^s decays into CP eigenstates adopts now a quite simple form 

_ A d & cos (AM*) + A%1 sin (AMt) 
CP[ ) cosh (Art/2) -^ Ar sinh(Art/2) 1 j 

which defines two quantities, the direct and mixing induced CP asymmetries, of capital 
importance in this thesis: 

- (3-68) 

^ - (3.69) 
1 + |A/| 
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and also 




(3.70) 



These three quantities are not independent, but they are related through 



I A dir 1 2 , | a mix 1 2 i I A 1 2 i 

|w4 CP | + \A CP | + I^IatI = 1 



(3.71) 



Let us take a look at the denominator in (13.6TI) . When measuring the direct and 
mixing induced CP asymmetries, the oscillation period is the relevant time, t osc <~ AM" 1 . 
During an oscillation period the denominator in (|3.67j) is 



And Ar /AM is very small for both Bj s systems. The time dependent CP asymmetry 
for B decays into final CP eigenstates can then be written in the following form 



In the case in which the relevant time is the lifetime, this is also basically true for Bj 
mesons, due to eq. ()3.64p . However, the possibility of a large width difference for the B s 
system (e.g. (|3.65p ). might make possible the measurement of the asymmetries ^Ar- This 
can be important for phenomenology, an example of which will be given in Chapter [6j 

3.4.2 The Bd, s — B^ s mixing angles and the AM^ mass differences 

The experimental facts that AM^ S 3> AI^ s and \r\d S ~ 1, have implications on the 
value of /M12 ■ I n fact, from eqs. ()3.20p and (13.211) one can solve for the real and 
imaginary parts of /M^ in terms of \r\j and AT^^/ AM^ S . In the limit |r|^ s — > 1 

and Ar^s/AM^s — > 0, there is a trivial solution, |r i2 s /-^i2 S | — ► 0. However, this is not 
the only solution. The set of two nonlinear equations always gives two solutions. The first 
is the one close to zero, and the second is the one that makes \T^2 /M^l a huge number, 
which for real data is of order ~ 10 3 . This second solution is not realistic; for example in 
the SM, |ri2/Mi2| ~ Tn\jm\ ~ 10~ 3 . Therefore, data implies that 



cosh (Ar r osc /2) - „4 A r sinh (Ar r osc /2) ~ 1 + O ( 



(3.72) 



Acp(t) = -^cp cos (AMt) + A^f sin (AMt) 



(3.73) 



r 



td,s 
12 



< 1 



(3.74) 



M\ 



12 



When this result is plugged in eqs. (j3.2Q[) and (|3.2ip . one gets 





(3.75) 
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such that 

AM d)S ~ 2|M^ S | (3.76) 
fa s ~ arg(M 1 t) (3.77) 

This means that in order to compute the mixing angles and the mass differences, it is 
enough to compute the mixing parameters ■ Moreover, the expression (|3.77|) is some- 
times taken as the definition of the mixing angle. One should keep in mind, however, that 
the mixing angle, as defined in this chapter, is not a physical quantity because it is sensi- 
tive to unphysical phase redefinitions. Therefore, it is important to work with a consistent 
convention for the weak phases everywhere; then the fact of assigning a numerical value 
to the mixing angle becomes sensible. 

In some particular cases, ona can extract the mixing angle in a very clean way from 
a mixing induced CP asymmetry. Consider, for example, a decay into a CP eigenstate 
which is dominated by a single amplitude, which means that the dominant parts of the 
amplitude have all the same weak phase. Then, working in a convention in which the 
global weak phase is zero, one has 

x = e -i4>M^l = ^ g-tof = Vfe -i4>M > ( 3 _ 78 ) 
A f A f 

and hence Aqp x = — ImA/ — Tjf sin <p M . The prominent example is the case of B d — > 
J/ipK s [57]. Since this decay is dominated by a single amplitude, then A j^ K / A j/ipK s — 
Vj/i>K s = ~~ 1) an d therefore, 

- A^ x (B d J/^K a ) ~ sin fa . (3.79) 

The neglected amplitude is both CKM and a s (rrib) suppressed with respect to the domi- 
nant amplitude, so the corrections to this equation are below the percent (or even the per 
mil) level [Ml W\ ED] . 



3.5 CP Violation in the Standard Model 

In this section we present the SM mechanism for CP violation; a simple and beautiful 
explanation formulated by M. Kobayashi and T. Maskawa in 1972 [72], only 5 years after 
the experimental discovery of CP violation. This formulation required also the existence of 
a third family of quarks (which had not yet been discovered), so together with the positive 
evidence of CP violation, this was a genuine prediction of the existence of the bottom and 
the top quarks. 

Apart from the strong-CP issue -which we will not address here, and in any instance 
is a negligible effect-, the CP violation in the SM comes from the electroweak sector, and 
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in particular from the Yukawa couplings. These couplings are described by the Glashow- 
Weinberg-Salam theory of EW interactions |914 \92\ 193|. 

The GWS theory is a gauge theory based on the SU(2)l <8> U(l)y group. The 'L' in 
SU(2)l indicates that the representations of the fermion fields depend on the chirality: 
Left-handed fermions are doublets and right-handed fermions are singlets. (This is how P 
violation is introduced in the SM.) Then, the kinetic terms in the lagrangian for the gauge 
and fermion fields are 

Ck,f = ^ kL {ifi + gjfeT a + g'$Y kL )^ kL +^ kR {i<P+g'$Y kR )ij kR (3.80) 

C K , G = -\{F a X n 2 - (3-81) 

where A^, are the gauge fields, T a are the generators of SU(2) in the fundamental 
representation, and a sum is understood over the left-handed and right-handed field species 
ki,R- Fav are the usual field-strength tensors and Y k are the hypercharges of the fermion 
fields. Both kinetic terms are invariant under the full gauge group. 

At this point a phenomenological problem arise. First, the chiral character of the gauge 
group spoils the gauge invar iance of the fermion mass terms, so unbroken gauge symmetry 
forbids massive fermions in this theory. Second, gauge invariance forbids also mass terms 
for gauge bosons. The problem is that the observed fermions and gauge bosons that we 
want to describe are actually massive. The solution within the GWS theory is to postulate 
that the gauge symmetry is spontaneously broken. 

The minimal realization of the spontaneous symmetry breaking (SSB) in the GWS 
theory is achieved by means of a single scalar SU (2) doublet. This scalar field <fi, is assumed 
to acquire a vacuum expectation value {4>) = v, breaking spontaneously three independent 
linear combinations of generators and leaving unbroken the fourth combination. The 
unbroken combination of generators gives rise to the corresponding combination of the 
gauge fields that constitute the physical (massless) photon. The three broken combinations 
give rise to the appearance of three goldstone bosons, according to the Nambu-Goldstone 
theorem. These three degrees of freedom become, in an appropriate gauge, longitudinal 
polarization modes for the other three independent combinations of gauge bosons, which 
then acquire a mass. These three combinations can be chosen such that they have diagonal 
couplings with the photon (a definite electric charge). These are the physical and Z 
bosons. The fourth remaining degree of freedom from the original complex scalar doublet 
is a physical real scalar field h, the higgs boson. In this way, the kinetic term for the 
scalar field <f> can be written as 

C KjS =(D^(D^) = ^(d„h) 2 + 

with mjy = vg/2 and mz = V\J g 2 + g' 2 /2 . This mechanism predicts a very definite 



m 2 w W+W-^ + -m 2 z Z^ 



1 + 



(3.82) 
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relation between the masses of the W and Z bosons that was verified precisely at LEP. 
This solves the issue of the massive gauge bosons. 

The kinetic term for the fermions (|3.80j) can then be expanded in terms of the physical 
gauge fields, and the charged and neutral weak currents identified. The electromagnetic 
currents then relate directly hypercharges and electric charges, and the hyper charge assig- 
nation can be made. The rest of the scalar lagrangian consists of an "opposite sign" eft 
mass term, pL(fy(j), and an interaction term, A</> 4 . In terms of physical fields these terms 
lead to interactions among gauge and higgs bossons, and to a higgs boson mass term. The 
parameters \i and A combine into v and run, the higgs mass. 

The issue of the massive fermions is still unsolved, and CP violation has not appeared 
yet in the theory. Here is where the things get interesting. For our purpose we focus on 
the quark fields Q % L = (JJ^D 1 ^) , U R and D R . Having introduced the scalar doublet, gauge 
invariance allows (and hence requires) the following couplings between the fermions and 
the scalar, 

C Y = -AjQE ■<t>D j R -\%e^Qi a -<t>luj i + h.c 

= -^DiD R (v + h)-^UiU R (v + h) + h.c (3.83) 

where a, b are SU(2) indices. The parameters X lJ are arbitrary complex matrices in 
flavor space. This piece of the SM lagrangian provides masses for the quarks, quark-higgs 
interactions, and flavor mixing. Also, because the couplings can be complex, introduces 
the possibility of CP violation. 

Now, to identify the masses, we must rotate the quark fields to get a diagonal mass 
matrix. The squared matrices X u -Xu and A^-A^ are hermitian, so they are diagonalizable 
with real and positive eigenvalues: A u -Ai = V u D^Vu, X u -X u = W u D^Wu, with D u diagonal. 
The same is valid for A^. So X u = V u D u wl and A^ = VdD^Wl. Then we rotate the fields 
such that 

Ul = C'4 , U R = W^u R 

Di = V^d{ , ^ = Wj'4 (3-84) 
In terms of the rotated fields ul,r and oIl,r, the Yukawa lagrangian (|3.83p becomes 

Cy = -mtf L <r R (l + ^) - (l + (3-85) 

where m* = vD™/y/2 are the masses of the quarks. 

This redefinition of the fields will only modify the fermionic part of the lagrangian Ck,f, 
that is, the fermionic currents. In fact, it only modifies the charged currents, because the 
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= -M-fD^ = ^uirV&udi (3.86) 



neutral currents are changed by things like V^V = 1. The charged currents become 

which defines the CKM matrix, V C ^ M = (V^V d )^ . We see that in this mass eigenbasis 
for the quark fields, neutral currents are still diagonal, but the CKM matrix generates 
flavor-changing charged currents. 

The CKM matrix can in principle be any square complex matrix in flavor space, since 
there is absolutely no constraint on the values of the Yukawa couplings X %3 . The only 
constraint is unitarity, since V^km^ckm = (yuV d ){y^V u ) = 1. We are ready now to see 
how CP violation arises in this picture. 

Consider first the case of 2 fermion families. The most general 2x2 unitary matrix 
can be parameterized by one angle and three phases, 

ntf c e'( a+7 ' cos 

However, we may rephase the three of the four quark fields to remove the three phases. 
This means that for 2 fermion families, the CKM matrix can always be chosen to be real, 

^ckm=( COS ° C R Sm t C ) (3.88) 
V — sin v c cos 8 C J 



_ . costf c e amv c e- , , , 



This is the Cabibbo matrix [94], and 9 C is the Cabibbo angle. The conclusion is that with 
only 2 fermion families, the SM cannot account for CP violation. This was the argument 
that led to the prediction of the third family of quarks, once CP violation was observed. 
For three generations of quarks, field rephasing cannot eliminate all the complex phases 
in the 3x3 CKM matrix, but one complex fase is left. Therefore, the SM with 3 families 
predicts CP violation provided this phase is not zero. Another conclusion is that CP can 
only be violated in processes that involve the three families simultaneously, so in the SM, 
CP violation is generally a loop effect. 

For numerical analyses, the most convenient parametrization of the CKM matrix is 
the standard parametrization [95]. However, a parametrization that makes its structure 
more clear is the Wolfenstein parametrization [96], 



v ud 




v ub 


v cd 


V cs 


Vcb 


v td 


Vts 


v tb 



1-iA 2 A AX 3 (p-irj) 

V CKM =\ V cd V cs V cb -A 1-|A 2 AX 2 \+0(X q 

AX 3 (l-p-in) -AX 2 1 

(3.89) 

where A = sin# c ~ 0.22 is treated as an expansion parameter. In this parametrization, 
the complex CP violating phase is 

7 = arg f-^J) = - arg (V ub ) + 0{X A ) = arg (p + i V ) + 0(A 4 ) , (3.90) 
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{p + in) 




Figure 3.3: The unitarity triangle, and the definitions of the angles. 

which up to order A 4 is localized in V u i> = \V u b\e~ t,y and V td = \V td \e~ tl3 , and the rest of 
the entries are real. Here, (3 is 

= ~ arg (~$\ = ~ arg (V td ) + 0(A 4 ) = arg (1 - p + irj) + 0(A 4 ) , (3.91) 

and it is zero whenever 7 is zero. Therefore, in the SM, CP is violated if and only if 
7^0. Experimentally it is found that 7 ~ 60° — 70°. A convenient graphical picture is 
provided by the Unitarity Triangle, which arises from the unitarity of the CKM matrix, 
in particular from the unitarity relation V ud V* b + V cd V* b + V td V* b = (see Fig J3.3j) . CP 
violation can then be quantified by the area of this triangle. SM fits provide already a 
very consistent and constrained Unitarity Triangle, as we saw already in Figj2j 

Let's now take a look at the mixing angles in the SM for B d and B s mesons. We saw 
in the previous section that for B q mesons the given definition for the mixing angle is 
practically equivalent to 4> M = arg(Mi2). We also saw in Section [3.1.11 that M\i contains 
contributions from Hfyf =2 and also from transitions with intermediate on-shell states, 
at second order in Hyf =1 (c./. eq. 13. 16|) . The contribution from intermediate on-shell 
states is mainly non-perturbative. It is an important contribution, for example, in kaon 
mixing. Fortunately, in the Standard Model the mixing of B q mesons is dominated by the 
perturbative AF = 2 box diagrams with a top quark in the loop (see Fig. I3.4j) . 

The CKM structure of these contributions is very simple. Both contributions are 
proportional to the product (V tq V* b ) 2 : oc (V td V* b ) 2 and Afg oc (V ts V* b ) 2 . So it's easy 
to see that the B d -B d mixing angle in the SM is given by 

0f = 2 (3 + 0(A 4 ) . (3.92) 

The mixing angle <j>^ M is zero is at this level of approximation, which means that enters in 
terms which are suppressed by at least A 4 . However, the phase itself is 0{\ 2 ). The angle 
f3 s is defined as 

As - arg (-tl)' (3 ' 93) 
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Figure 3.4: Dominant SM contributions to the B q -B q mixing parameter Mf^ ■ The CKM structure 
has been made explicit. 



Up to C(A 6 ), only Vt s acquires a phase, so Vt s = —\Vt s \e l/3s . In fact (3 S ~ —\ 2 n. Then, 
the B s -B s mixing angle in the SM is given by 

4> S S M = 2/3, + 0(A 6 ) . (3.94) 

This is the way in which the mixing angles are related to the angles of the unitarity 
triangles in the Standard Model. 



Appendix: CP-averaged branching ratio for 2-Body decays 

The differential decay rate dT for the decay of a particle P in the center of mass (CM) 
frame into a final state / is (see, for example, |97j ) 



where i runs over the set of particles in the final state with 4-momenta pi in the CM frame. 
A(P — > f) is the amplitude of the process, and m P is the mass of the initial particle P 
(the initial energy in the CM). For the case of 2 particles in the final state, the 4-momenta 
of the final particles are 

pi = (m P /2 + A , ku) 

p 2 = (m P /2 - A , -k-u) (3.96) 

as described in Fig. 13.51 

We can use the 4-delta function to integrate over the f>2 3-momenta and the modulus 
\pi\ in (|3.95p . This leads to the following general formula: 

^ =^—\A(P^f)\ 2 (3.97) 



dfl J CM 32-7T 2 



771? 
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Figure 3.5: Kinematic description in the CM frame of the decay of a particle P into two final 
particles. Conservation of 4-momenta imposed by the delta function in (|3.95|) implies the final 
state given by (|3.96|) . with A and ku given in this figure. The axis z is chosen arbitrarily. 



If the initial particle P has spin 0, the process is spherically symmetric in the CM frame 
and the amplitude A(P — > f) cannot have an angular dependenc^l. In such a case the 
amplitude factorizes and the integrated decay rate can be written as 

r(P - /) = f PS \A(P - /)| 2 (3.98) 

where / PS is called the phase space factor. This is just the angular integral which depends 
only on kinematic variables: 

= v.^k (3 "> 

Here a s is a symmetry factor, which is a s = 1 if the two final particles are distinguishable 
and a s = 2 if they are identical. The meaning of this factor is that the final state with 9 
and that with 9 + tt are indistinguishable, and must be counted only once. 

In the case of 2 particles with the same mass in the final state, we have mi = m2 = m,f, 
A = 0, and k 2 = (m 2 p — 4my)/4, and the phase space factor is given by 



The branching ratio BR(P — > /) is defined as the partial decay rate T(P — > /) divided 
by the total decay rate T(P — > all). This quantity is naturally normalized to one, since 
the sum of the decay rates of all the exclusive decay modes of P must equal Y{P — > all). 
The total decay rate is the inverse of the lifetime r P . In "untagged" neutral meson decays, 
there is no information of whether the decaying particle is a P° or a P°. In such a case 



6 In the case in which P has spin, but the polarization state is unknown, we must average over the spin 
states and we recover spherical symmetry. 
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the measured branching ratio is an average between the decays P° — > f and P° — > f, 

BR(P° - /) = r(P ° 7/^% ~ ^ = Tp/ps (l ^ |2) (3 - 101) 

where (\A f \ 2 ) = (\A f \ 2 + \A f \ 2 )/2 = {\A f \ 2 + \Af\ 2 )/2, the last equality being true for / a 
CP eigenstate. When the final state is a CP eigenstate then this is called a CP-averaged 
branching ratio. In this thesis we will only deal with CP-averaged branching ratios, so we 
will call them just BR, without the bar BR. Moreover, we define the g PS factor as 

Ops = Tpfps (3.102) 

Some numerical values for these factors, that will be used later, are 

g PS (B% -> vr+vr") = 8.8 x 10 9 GeV" 2 

spsGbS K ° R0 ) = 8 - 6 x 109 GeV ~ 2 

g PS (B° s -> K + K~) = 8.1 x 10 9 GeV" 2 

g PS (B° s -> K°K°) = 8.1 x 10 9 GeV" 2 (3.103) 

for input values of masses and lifetimes in [71] and using the conversion factor 

1 sec = 1.5 19 x 10 24 GeV" 1 (3.104) 

For B decays into two light mesons, the masses of the final mesons are sometimes neglected 
in front of me- If this is done, there are only two g PS factors: g P s{Bd) = 8.8 • 10 9 GeV -2 
and g PS (B s ) = 8.2 • 10 9 GeV~ 2 , which often simplifies the numerics. 
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Chapter 4 

Measuring New Physics in B 
Decays 

Up to now, the theoretical analyses of CP violation in B decays have been mostly focused 
on two complementary subjects. First, many different methods have been proposed for 
extracting in a clean way the CP violating angles a, (3 and 7, within the SM (see e.g. 
[711 \M\ ) . We have now many ways of checking the compatibility of different measurements 
within the SM and detect the presence of NP. Second, many studies have been made of 
observables within specific scenarios of NP in the forward direction; that is: what values 
for these observables does this specific model predict? These second type of studies are 
important because they provide a direct bridge between new data and the reduction of 
the allowed regions in the NP parameter spaces^. In an extreme situation, they may even 
rule out some models (see for example [98]). 

However, the identification of the NP, once this is found, requires NP studies in the 
inverse direction; that is: given this non-standard experimental data, what is the NP that 
describes it? This "inverse problem" is rather complex because of two reasons. First, 
given a specific model, a signature in observable space is generally associated to many 
points in parameter space ( "degeneracies" ) . In the most unfavorable situation this can be 
a huge continuously connected region, meaning that this signature in observable space will 
give absolutely no clue on the basic structure of the underlying theory (for a study of the 
inverse problem in supersymmetry see |99|). Second, the most important information that 
must be extracted from NP signals is related to general characteristics of the NP. Indeed, 
solving the inverse problem in supersymmetry is of no help if nature turns out not to be 
supersymmetric. But model independent analyses of NP are generally unapproachable 
because of the large number of parameters involved (see the discussion in Section [1. 3D . so 

1 We use the plural in "parameter spaces" because each NP scenario has its own parameters, and thus 
its own parameter space. 
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at the end one has to come up with arguments that allow to reduce the number of model 
independent NP parameters. 

A first step in this direction was taken in Ref. [100| . where it was shown that an argu- 
ment concerning new physics strong phases allows to reduce the number of NP parameters 
to a manageable level, and measure them. The knowledge of these model independent NP 
parameters allows then to identify partially the new physics, and to establish some of its 
generic features. 

In this chapter we present an approach, within this framework, to measure the NP 
parameters that appear in B s — > K + K~ amplitudes. This approach requires experimental 
information on branching ratios and CP asymmetries in Bd — > ir + ir~ and B s — > K + K~ , 
and is based on a SU(3) analysis of these modes and the (motivated) assumption that 
any NP contributing to b — > d penguin amplitudes is negligible at the considered level of 
approximation. 

4.1 Parametrization of the SM Amplitudes 

The SM amplitude of a generic two body non-leptonic B decay (B g — > M1M2) is commonly 
parameterized in terms of its CKM structure. We used this fact already in the discussion 
of strong and weak phases in Section I3T31 where it was justified by invoking the form of the 
SM effective Hamiltonian as written in Eq. (jl.27p . However, we have never proved here 
that the Wilson coefficients in Eq. (|1.27j) are CKM independent; in fact, strictly speaking, 
they are not. 

In this chapter we abandon temporarily the effective Hamiltonian description of the 
amplitudes, in favor of a diagrammatic approach that makes this CKM issue more trans- 
parent. 

We will restrict ourselves to a general situation in which the parametrization of Eq. 
( 13.590 holds: decays which contain contributions from tree or penguin diagrams where the 
two up-type quarks produced in the decay (if any) are both either up or charm quarks. 
The leading contributions, which contain a single W propagator, are then all proportional 
to either one of the CKM products V u DV* b , V C DV* b or Vt£>V^, where D is the down-type 
quark in the final state which is not q (see Fig. I4.ip . 

Because in a given decay the CKM elements always enter by pairs, it is customary to 
define the following product of CKM elements, A^ D> = V^,V m . Moreover, the unitarity of 
the CKM matrix implies the triangle relation 

A£» + + \p = (4.1) 

which allows to eliminate the CKM factor from the amplitude. The amplitude is 
then a sum of two terms, proportional to X^ and Xi D) respectively. We define the "tree" 
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Figure 4.1: Leading contributions to B q — > M\Mi which do not involve the spectator quark 
(b — > Dqq). The quark labels stand for q,Q — {u,d,s,c}, U,U' = {u,c} and D = {d, s}. The 
color allowed (a) and color suppressed (b) tree amplitudes are proportional to VubV^. The color 
suppressed and color allowed penguin amplitudes (c) and (d) are proportional to VenV£, where 
£ = {u, c, t} is the quark inside the loop. Note that not all of the diagrams contribute to a given 
decay. 



(^AfiM 2 ) an d the "penguin" (Pm iM2 ) amplitudes as the coefficients of the two terms with 
different CKM structures. Thus the amplitude can be written as 

A(B q - M X M 2 ) = \£» T A \ M2 + \ { c D) 'PS llM2 (4.2) 

Besides the topologies shown in Fig. 14. 1\ there are also leading contributions which 
involve the spectator quark, mainly exchange, annihilation and penguin annihilation di- 
agrams [101] . However, it has been argued that these can be neglected at leading order, 
since they imply the fairly unlikely situation in which the two quarks inside the B meson 
annihilate each other or interact weakly before the decay of the b quark. This is called 
the spectator quark hypothesis, and as any other hypothesis it should be handled with 
care (for example, these contributions could be Cabibbo enhanced). Unlike the penguin 
annihilation topology, the exchange and annihilation contributions do not in general obey 
the CKM structure discussed above, and do not fit in the scheme. However, for charmless 
decays and decays into CP eigenstates they do, and they can be included in the description 
trivially. Since we will be focusing our discussion in two-body decays of neutral B mesons 
(q = d, s) into CP eigenstates, this issue will be of no concern. 

To be specific, the tree and penguin contributions, in terms of the diagrams in Fig. 
14.11 are given by 

Tm x m 2 = T u + C u + (P u — P t ) + (Pu - Pf) 

P£ hM2 = T c + C c + (P C -P t ) + (PC -P t c ) (4.3) 

where X^^Tjj, X^^Cjj, A^ D) *P^ and \g*Pf are respectively the amplitudes (a), (b), (d) 
and (c) in Fig. 14.11 
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Figure 4.2: Exchange (a), Annihilation (b) and Penguin Annihilation (c) contributions, (a) 
and (c) only exist for q = {d, s} and (b) only exists for q = {u, c}. The penguin annihilation 
contribution always fits in the structure (|4.2p . The exchange and annihilation topologies only have 
such CKM structure for U = U' and q = U respectively. For final CP eigenstates this structure is 
automatically satisfied (there is no annihilation contribution in this case). 



The CKM products A[f' have the following numerical values |71j : 



= 0.0038e^ = 
0.00088e^ 



|Ai s) |e-^ 



AW) = -0.0094 = -|A^| 
= 0.04 = |A c s) | 



(4.4) 



where 7 is the CKM weak phase. Note the negative sign of Ac . We denote the relative 
(strong) phase between T^ iM2 and f$ lAfa by 



= arg 



P q 

1 M 1 M 2 



(4.5) 



Now we can give a general useful parameterization of the amplitude in (|4.2p 

\i D) 



A{B g MiM 2 ) = \X ( u D) \\T\e ldT e 47 ± 



(D) 



JO 



\(°)\\T\ e ldT e^±de 



(4.6) 



the plus (minus) sign being for D = s {D = d). We have dropped the labels q and M1M2 
for simplicity, and it will be done so systematically whenever there is no ambiguity. We 
have also defined [64j 

\ ( c D) 



d = d? r _ 



ID) 



p q 

r M 1 M 2 



rpq 



(4.7) 



and St is the strong phase associated to T£ lM2 ■ Note that when (|4.6p is the full amplitude 
(that is, when there are no NP contributions), the phase St is not physical and does not 
appear in the observables. The amplitude for the CP-conjugate process B q — > M\Mi shall 
be obtained from (|4.6p by changing the sign of the weak phase 7. 
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4.2 Measuring the SM parameters 

Assuming that no new contributions come into play, and that the full amplitude is de- 
scribed within SM as in (|4.6j) . it is possible to measure the SM hadronic parameters TJ^ lM2 , 
Pm ± m 2 an d $m x m^- There are three potential experimental measurements that can be car- 
ried out in neutral B decays that allow us to extract this information. These are the 
CP-averaged branching ratio (see the Appendix to Chapter [3]), the direct CP asymmetry 
and the mixing induced CP asymmetry (see section [33]). By virtue of the formulae (|3.10ip . 
(|3.68|) and (|3.69|) one can relate these observables to the theoretical parameters in (|4.6p : 

BR = 5ps |A^ } | 2 |T| 2 (1 + d 2 ± 2(icos7cos#) (4.8) 

Adir _ i 2d sin 7 sing 
" A cp-± 1 + (i2±2dcoS7COS ^ 

ix _ sin(27 + M ) ±2dsm(j + <j) M )cos8 + d 2 sm<j) M 
1 + gt ± 2d cos 7 cos & 

where, again, the plus and minus signs correspond to D = s and D = d, respectively. 

The hadronic parameters |T|, d and 6 are specific of each particular decay. The CKM 
angle 7 and the mixing angle <j) u are, on the contrary, universal theoretical parameters, 
and can be obtained from a multitude of different decays. Therefore, we are left with three 
hadronic parameters (|T|, d and 0) to be obtained from three experimental values (BR, 
and A^ x ) by means of equations M- flCT]) . 

As a matter of illustration and for later usage, consider the decay B® — > ir + ir~ . At the 
quark level (see Fig. I4.3|) it corresponds to a b — > d transition (D = d), so we must choose 
the minus sign in (|4.8p - (|4.10p . For the branching ratio and the CP asymmetries we take 
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Figure 4.4: Allowed regions in the d — 6 plane for B® — > tt + t7 , (a) taking 7 = 61° and (b) 
varying 7 in the range 7 = (61^5)°. The dashed lines show the central value and ±ler deviations 
for Aqp (B° 7T+7T-), and the solid lines the same for A^{B° d 77+77-). 



an average of the Babar and Belle data [1021 11031 1 104] , 

BR{B° d ^ vr+vOcxp = (5.0 ±0.4) x 10~ 6 (4.11) 
-4cp(#° ^ + Ocx P = -0.33 ±0.11 (4.12) 
A^ x (B° d - 7r+7r-) exp = 0.49 ±0.12 (4.13) 

Because .A^p and ^'gp x ' do not depend on the magnitude of the tree (c.f. (|4.9p . (|4.10p ). 
the calculation of the hadronic parameters can be divided in two steps. First we vary the 
asymmetries inside their experimental range to find the allowed region in the d — 9 plane. 
This is illustrated in Fig. I4.4L Second, we vary the branching ratio in its range together 
with the allowed values of d and 9, to find the range for \T\ according to (|4.8]> . The results 
are shown in Table I4~T1 

Note that two different solutions in the d — plane are allowed by a single value of 
the pair (A^p, A™™) (see the two separate regions in Figs. I4.4p . This is due to the 







6 n7T (deg) 


|T^|(10- 6 GeV) 


| -^TTTT I -^7T7T | 


7 = 61° 


(0.18,0.49) 


(107, 150) 


(4.87,6.11) 


(0.07,0.20) 


7 = (6itir 


(0.14,0.61) 


(88, 158) 


(4.60,6.45) 


(0.06,0.25) 



Table 4.1: Allowed ranges for the hadronic parameters in the decay — * 7r + 7r~, calculated 
from the experimental values for the branching ratio and the two CP asymmetries. The values are 
taken from an average of Babar and Belle data. The sensibility on the error of 7 is shown. The 
penguin-to-tree ratio is obtained from d and (|4.7I) . 
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non-linearity of the equations (|4.9p and (|4.10p . However, it should be clear that only one 
solution corresponds to a given pair (T,P), so only one of the two regions in the d — 9 
can be physical. QCDF predicts a value for d n7T around 0.3 [26], so the second solution 
(dim > 3) is likely to be the non-physical one and it's ruled out. This is also confirmed 
independently by an analysis concerning the channel B d -> tt =f K ± [T05] . 



4.3 Parameterization of the New Physics Amplitudes 

Besides the SM contributions, a given decay might also receive contributions from physics 
beyond the SM. Without a specific model for the NP, they cannot be calculated. However, 
it is possible to make a model independent parameterization of these contributions under 
certain general circumstances [lOOj . 

In the context of the low energy effective Hamiltonian, the new contributions will be 
of two types. First, there will be new contributions to the Wilson Coefficients of the SM 
operators. Second, new operators will appear which are not generated by the SM. To this 
end we write the effective Hamiltonian as a sum of the SM and NP contributions: 



(4.14) 



Because the nature of the NP is unknown, TC NP contains all possible operators compatible 
with a given process. As it will become clear below, we are mainly interested on the 
NP contributions to b — > sqq transitions. There are 20 dimension-six operators that 
contribute potentially to such decays, and the general NP effective Hamiltonian in (|4.14p 
can be written as 



njq 



A,B=L.R 



Gf 
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fq,l SalAb/3 q^BQa + fq,2 S^ A bqiBq 

+9q,i Sa^jAbp qpl^lBqa + 5gJ s^jAbqi^Bq 

+h^f SaO-^syAbp qpo-ftvjBqa + K B 2 sa^jAbqa^sq (4.15) 



where 7 H , L = ^(1 ± 75). The coefficients in ()4.15|) will in general contain new (model 
dependent) weak phases and the matrix elements of these operators will introduce new 
(process dependent) strong phases. The strong phases, as mentioned in section [3T3l come 
mainly from rescattering. The important point is that the NP strong phases must come 
from rescattering of NP operators. Let's discuss this in more detail. 

In the SM, the strong phases may be large and must be taken into account because they 
arise from rescattering of tree operators, and although the process of rescattering is size- 
consuming (the rescattered penguin amplitudes are about 5-10% of the tree amplitudes), 
the tree amplitudes are big and the rescattered penguins can be as large as the SM penguin 
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amplitudes. This situation holds even when there are no tree operators that contribute 
directly to the decay, since these can contribute through rescattering. In such cases, the 
rescattered penguin can be of the same order of the largest contributions. 

The NP operators, however, are expected to give direct contributions not larger than 
the SM penguin contributions. This means that the rescattering from such NP operators 
will be at most 5-10% as big as the SM penguins, and at the level of other neglected dia- 
grams. We can therefore argue that these small contributions can be neglected, and since 
these are the main sources of new strong phases, the strong phases can be approximately 
set to zero in the NP amplitude. Thus, without strong phases, the NP amplitude can be 
parameterized in terms of a single effective weak phase: 

{f\HUB) = A^ (4.16) 

Note that there is a potential loop-hole in this argument. We can be sure about the size 
of the largest NP contributions to the decay in consideration, and thus be confident in 
neglecting the strong phases coming from the rescattering of these diagrams. However, the 
possibility remains of the existence of NP diagrams which do not contribute directly but 
through rescattering to the given decay. If they were much larger than the SM penguins 
that do contribute, they could introduce sizeable strong phases in the NP amplitude. 
Other cases which get around this reasoning are when the NP is light, or if there is a 
significant enhancement of certain matrix elements. 

4.4 Hints of New Physics in B decays 

The question arises of whether we should expect sizeable contributions from physics beyond 
the SM, or if these, even if they exist, are too small to be observed at present. A second 
question would be what transitions are likely to be appreciably affected by NP. There 
are several indications that may help to answer these questions. For example, the CP 
asymmetry in b — > sqq modes (q = u, d, s) is found to differ from that in b — > ccs decays by 
2.6cr (they are expected to be approximately equal in the SM) [85 11106] . In addition, some 
B — > ttK measurements disagree with SM expectations [1071 1108] , although the so-called 
B -> ttK puzzle pIU Q3U1 [ml QUI [H2] has been reduced pHl Ell QUI EBj. One also 
sees a discrepancy with the SM in triple-product asymmetries in B — > 4>K* |117tlll8|lll9| . 
and in the polarization measurements of B -> <f>K* [12Q1CE21] and B -> pK* [7X1 fT22| fT23j . 
These discrepancies are (almost) all not yet statistically significant, being in the l~2a 
range. However, if these hints are taken together, the statistical significance increases. 
Furthermore, they are intriguing since they all point to New Physics in b — > s transitions. 
For this reason it is interesting to consider the effect of NP on B decays dominated by the 
quark-level b — > s transition, and assume that no NP affects b — > d transitions. 



4-5 Measuring the New Physics parameters through Flavor Symmetries 



81 



4.5 Measuring the New Physics parameters through Flavor 
Symmetries 

For the sake of clarity we begin this section summarizing the basic ideas discussed so far. 

The amplitudes for B decays that at the quark level are governed by the transition 
b — > Dqq, can be written as 

A(b - dqq) = \>tf\\T\ (e* - de w ) (4.17) 
A(b sqq) = \X^\\T'\e i5T (V 7 + d'e i6 '^ + A q t?** (4.18) 

These parameterizations rely on several conditions: 

• SM diagrams which do not obey the CKM structure in (|4.2p can be neglected with 
respect to the penguin amplitudes. 

• Potential NP contributions to b — > d decays can be neglected with respect to the SM 
penguin amplitudes. 

• There are no NP operators large enough to produce a sizeable contribution to the 
b — > s decay through rescattering. 

In section 14.21 we showed how to measure the SM hadronic parameters in (|4.6p from 
measurements on the branching ratio and the CP-asymmetries. The method was applied to 
obtain the hadronic parameters of the decay B^ — > ■k + tt~ , that according to the previous 
conditions it is described entirely within SM as in (|4.17p . The same could be done, in 
principle, with a b — > s decay. Taking the amplitude in (|4.18|) . the branching ratio and 
the asymmetries can be expressed in terms of the SM hadronic parameters and the NP 
parameters. The expressions are much more complicated in this case: 

BR = BR SM • (1 + B NP ) (4.19) 

ASM I -T)NP 

where the SM functions BR SM , A b d f r and A s ™ ix are those in (^8]) - (^T0|) taking the plus 
sign and D = s, and the NP functions are given by 

fi NP = -^-{z 2 + 2z[cos(<5> q -7)cos5 T + dcos<S> q cos(6 + 5 T )]} (4.22) 
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Figure 4.5: Leading contributions to the decay B® — > K + K~, (a) tree and (b) penguin contribu- 
tions. These diagrams are related directly to those of — > tt + it~ (sec Fig. I4.3[) by U-spin (mainly 
a permutation d <-> s). 



£> NP = _ {2z[sin($ 5 -7)sin«5 T + dsin$,sin(0 + <yT)]} (4.23) 

X NP = ^{z 2 sin (2$ 9 + S ) + 2z[sin ($, + 0, + 7 ) cos<5 T + dcos (0 + *r) sin (*, + S )] 
+ [cos (2 7 + <f) s s M ) + 2d cos (0) cos (7 + (f) s s M ) + (i 2 cos <^ M ] sin<fyf } (4.24) 

The B® — B® mixing angle 4> s can also suffer from the same NP contributions. We take 
this into account by adding a new contribution^], (j) s = + 50^ p . In order to simplify 
the expressions we have defined 

A = 1 + 2<i cos 7 cos 6> + (i 2 , z = A q /\\ { U D) T\ . (4.25) 

Note that in the limit A q , 5(ff s p -> the NP functions B NP , V NP and M NP go to zero and 
the observables in (|4.19p - (|4.2ip reduce to the SM functions. 

The CKM angle 7 is considered an input parameter and it is obtained from the SM 
fit [HlEll, 7 = (6Fts) - The 

mixing angle <j) s (including its NP contribution) can also be 
measured independently (see Chapter E]). However, three measured observables are not 
enough in this case, since we have six parameters to determine, \T\, d, 6, 5t, A q and <& q . 

Fortunately, flavor symmetries allow us in certain situations to obtain the parameters 
\T\, d and 6 from a different decay, related to the first by a symmetry such as U-spin. If 
this can be done, then 5t, A q and <& q can be obtained from (|4.19|) - (|4.21|) , This strategy 
has been used in the literature for the SU(3)-related decays B® — ► ir + 7r~ and B® — > K + K~ 
[64" 1 11241 H25| . although in pursue of slightly different goals. 

In [M] (see also [1261 1127] ). the relationship between these decays was used to extract 
the CKM angles (3 and 7 in the SM. Assuming no NP in B® — > K + K~ , and obtaining the 

2 Note that 4> M in (|4.10|l should be understood in A s ™ ix as <j> M = <^ M 
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mixing phase 4>l M independently, (3 (or = 2/3) and 7 can be obtained from A.Qp(Bg — > 
7r + 7r _ ) and ^Qp^-Bg — > 7r + 7r~) once the parameters d and have been obtained from 
A^(B° S -> #+#-) and -4^(5° K + K~). There is, however, two-fold and four-fold 
ambiguities in the determination of 7 and /3, respectively. 

In [124] . the branching ratio for B® — > K + K~ was computed from (|4.8|) . also within 
the framework of SM, once the hadronic parameters d, 9 and \T\ were obtained from 
B® — ► 7r + 7r~. The prediction was compared with the experimental value by CDF, looking 
for a clear discrepancy that would signal NP. Due to the present large experimental errors 
no discrepancy can be drawn, but once this errors get squeezed, this method will certainly 
reveal if there is need for a nonzero value of A u . 

The first time the NP amplitude for B® — > K + K~ was included in the scheme was in 
[125] . Although the lack of experimental data on the CP asymmetries for B® — ► K + K~ 
makes it impossible to apply the method to extract the NP parameters, predictions can be 
computed for the B® — > K + K~ observables once a NP model has been used to compute 
the NP amplitude. This was done in refs. [1281 [36] and will be explained in Chapter [7[ 

Up to this point, let's explain the method in somewhat more detail. The decays 
B® — > 7r + 7r~ and B® — > K + K~ are described at leading order by the diagrams in Figs. 14.31 
and 14.51 respectively. In the U-spin limit the contributions are equal except for the CKM 
elements. This means that T KK = T wn and P KK = P wn hold up SU(3) breaking effects. 
The SU (3) breaking is accounted for in the deviations from unity of the ratios 

_ T KK P KK T KK j^g I A nn\ 

r T = Tfr- > t pt = p /T = & ■ (4.26) 

7Y7V +7Y7T / 7T7T 

An estimate of this U-spin breaking within QCDF can be found in |129j . The SU(3) 
breaking in r T is expected to be much larger than that in r PT since both factorizable and 
non-factorizable terms are present in the former, while factorizable contributions cancel in 
the later. A more ambitious calculation of the ratio \r T \ using QCD sum rules [UJ gives 

\r T \ = 1.52±g;i| (4.27) 

A fair estimation for the SU(3) breaking parameter r PT is still lacking. In section H~6l we 
will make a quick comment on this issue. For the moment it will be assumed that this is 
not more than ±20% in magnitude, 

£ = 1.0 ±0.2 (4.28) 

and neglect the U-spin breaking in its phase (AO = 0), whose impact has been shown to 
be very small [127] . 

With these values for the U-spin breaking parameters and the results for {T^nl and 
IPtt-k/T^I in Table [O] it is possible to calculate the hadronic parameters for B® — > K + K~ 
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(10" 6 GeV) 


\Prk /T KK | 


7 = 61° 

e = i 


(8.57, 10.75) 


(0.07,0.20) 


7 = 61° 
£ = 1±0.2 


(8.57,10.75) 


(0.06,0.24) 


7 = (61^)° 

e = i 


(8.10,11.35) 


(0.06,0.25) 



Table 4.2: Allowed ranges for the hadronic parameters in the decay B® — > K + K~ . These are 
obtained from Table EP1 by U-spin correspondence as explained in the text. We take the central 
value in (|4.27p for \tt\- The impact of the uncertainty in the U-spin breaking parameter £ and 
CKM-angle 7 is shown. 





(xlO 6 ) 


ds SM 


ASM 
dir KK 


ASM 
mix KK 


7 = 61° 
$ = 1 


(6.4,42.6) 


(1.2,9.3) 


(0.15,0.45) 


(-0.32,-0.10) 


7 = 61° 
£ = 1±0.2 


(4.2,61.9) 


(0.8,13.5) 


(0.12,0.56) 


(-0.38,-0.09) 


7 = (61^)° 

e = i 


(5.0,60.7) 


(0.9,13.2) 


(0.08,0.58) 


(-0.34, +0.08) 



Table 4.3: SM predictions for the branching ratio and CP-asymmetries of B® — » K + K . 

using (|4,26p . These are shown in Table 14.21 Once the hadronic parameters have been 
calculated, we can give SM predictions for the observable^! [1281 136] - These predictions 
are given in Table l4~3l (see also [105] for comparison), and the correlations between them 
are illustrated in Figure 14.61 for different values of the U-spin breaking parameter r PT . 

These results can be used as a consistency check when comparing with experiment, or 
can be used together with experimental data to extract the NP parameters A u , $ n and 
5<r- Up to date the only available piece of data on B® — ► K + K~ is the value of given by 
CDF [130] . 

BR(B S -» K + K~) cxp = (24.4 ± 1.4 ± 4.6) • 10~ 6 (4.29) 

which fits fairly good in the SM prediction, and implies a value for B NP compatible with 
zero, 

B NP e (-0.54,3.56) (4.30) 

3 Regarding BRf^, it is sometimes more useful to present the ratio of branching ratios of B® — > K + K~ 
and B° d -> tt+tt": R 3 d = BR(B° S -> K+K-)/BR{B° d -> tt+tt") [T25]. 
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Figure 4.6: Correlations between the observables A s d f r KK - BR S K ™ and A s d f r KK — A^- x KK , for 
7 = 61° and £ = 1, £ = 0.9 and £ = 1.1. 



but with large errors that still allow for the presence of NP. 



4.6 A glimpse further 

The method described above to measure the NP parameters would rule out automatically 
all those NP models that fit the conditions listed in Section 14.51 but do not reproduce the 
measured values for these parameters. To that end one should compute the contributions 
coming from each particular model to be tested, which can be quite a tedious job. There 
are, however, several generic features of the NP flavor structure that could be identified at 
once by the power of flavor symmetries. All NP models which do not share these features 
would be ruled out. 

Consider the decay B® — > K°K°. Only the diagrams (c) in Fig. 14.11 contribute to this 
decay, so only the penguin P^o^o can be related to P^+ K - by U-spin, and not the tree 

^k°k ( wrncn does not contain the "pure" tree (a)). However, because |Au | <C |Ac |, the 
amplitude for this process can be approximately written as 

A(B° S K°K°) = X^P^ko + A d e m * (4.31) 

Then measurements on B® —> K°K° allow to extract the NP parameters A d and &d- 

Consider the case in which the NP is isospin conserving. This means that whatever 
this NP is, it should predict A d = A u and &d = <£ u . Now, if the measured values for the 
NP parameters do not fit this rule, then all isospin conserving NP models can be discarded. 

The main advantage of using experimental measurements on a particular decay to- 
gether with flavor symmetries to extract information in other channels, is that the hadronic 
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effects are taken into account in a model independent way. The price we pay is the large 
uncertainties coming from the SU(3) breaking effects (up to 30%), which are still not under 
control. 

A different approach is to use a particular theoretical framework to calculate the matrix 
elements, but the inclusion of hadronic effects always leads to model dependence. It would 
be a step forward if both approaches could be combined to reach more reliable predictions. 

An attempt in this direction has been made in [35] using QCDF. According to the 
discussion at the end of Section 12.1.41 and as will be explained in detail in Chapter (5J the 
quantity 



From (|4.32p and the measured values for the branching ratio and the direct CP-asymmetry 
for this decay, one can obtain its hadronic parameters T^. ^. and P^oro- QCDF predicts 
a relationship between these and the hadronic parameters for B® — > K + K~ , which is now 
free of most of the model dependence, and thus much more reliable than the pure QCDF 
computation. 

As an example, one can obtain d s K+K _ in this way, and together with from 
Table 14.11 obtain a value for the U-spin breaking parameter £. This is described in detail 
in Chapter [5j The result as given in ref. [35] is 




i = 0.81 ±0.35 



(4.33) 



This value can then be used in the analysis of Section 14.51 instead of the poorly justified 
one taken above. However, one should consider carefully up to what point these two 
approaches are independent when using the value in eq. (I4.33j) for £. 



Chapter 5 



Exploring B —> KK Decays with 
Flavor Symmetries and QCDF 

Along the past chapters of this thesis we have seen continuously the difficulty and the 
phenomenological importance of the computation of the hadronic parameters ( "tree" and 
"penguin") in two body non-leptonic B decays. We have seen that this type of decays 
contain a huge amount of information that can be used to understand and test the SM 
and to look and measure NP; but we have also seen that in order to do that the effects of 
strong interactions have to be controlled in some way. 

Moreover, the increasing precision that is being achieved experimentally, and the op- 
timistic prospects for the future coming from the successes of the CDF collaboration, the 
starting of LHC, and the possibility of a Super-B, urges the development of theoretical 
tools improving on the treatment of hadronic uncertainties. 

We have already seen several ways in which hadronic parameters in two-body non- 
leptonic B decays can be extracted. First, in Section 14.21 we have seen how these can 
be extracted directly from data. The argument is that three hadronic parameters can 
be extracted from three independent observables, which for neutral B decays are the 
branching ratio and the two components of the time dependent CP asymmetry. One just 
has to invert the equations (|4.8p - (|4.10p and solve for the hadronic parameters. For the case 
of Bd — > 7r + 7r~, the results were given in Table I4~T1 One should be careful, though, when 
using this procedure in NP analyses, since this is only correct when the decay at hand is 
not affected by the NP. In this thesis we bypass this potential problem by assuming that 
any NP affecting b — > d transitions is negligible at the considered level of accuracy. 

A second approach that may be used to extract the hadronic parameters is QCDF 
as outlined in Sections 12.1.31 and 12.1.41 The amplitude of a B decay in QCDF is given 
by the matrix elements in eq. (|2.43[) . Therefore, the hadronic parameters are given by 
T MlM2 = (M^TX + T^B) and P MlM2 = (MiM 2 |T| + T£\B) (see Section ECU) • For 
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Theory 


S2 S3 S4 


|7W| (10- 6 GeV) 


7 7+0.04 
' • ' -0.06 


6.3 7.8 6.8 


| -^7T7T /-^7T7T | 


u - oz -0.09 


0.49 0.37 0.48 



Table 5.1: Hadronic parameters for B c i — >• tt + tt from QCDF. S2, S3 and S4 correspond to 
different scenarios for the inputs as described in ref . [26] . 

the Bd — > tt + tt~ example, the hadronic parameters predicted by QCDF are shown in 
Table [5TT1 126| . The corresponding QCDF prediction for the B^ — > 7r + 7r~ branching ratio 
is 

BR(B d - tt+Oqod* = 8.9±|:2 -il -H * 10 ~ 6 (5- 1 ) 

which gives a central value which is high compared to the experimental number (14. 111) , 
but with tremendous uncertainties. Clearly this is presently a problematic decay for the 
QCDF approach. 

An alternative way to extract the hadronic parameters with the aim of giving predic- 
tions for the observables, is the flavor symmetry analysis described in Section 14.51 The 
process consists in establishing SU (3) relations between the hadronic parameters such that 
the knowledge of some hadronic parameters (that are extracted from data) allows to know 
the other hadronic parameters on which we are interested. Then, after a sensible inclu- 
sion of SU(3) breaking effects, one can calculate observables. The analysis in Chapter [4] 
relating Bj — ► ir + ir~ and B s — > K + K~ gave the predictions plotted in Fig. 14.61 

A comparative plot between the predictions for B s — ► K + K~ from QCDF and flavor 
symmetry is shown in Fig. 15.11 As can be seen, the theoretical predictions have uncer- 
tainties much above the experimental precision. Moreover, this decay will probably be 
measured with better accuracy in Tevatron shortly, including the direct CP asymmetry. 
In LHCb both observables will be measured with very high accuracy. Therefore, it is 
mandatory to improve substantially the theoretical predictions if the experimental efforts 
want to be exploited. But, how can this be done? To answer this question it is necessary 
to understand what are the advantages and limitations of each of these approaches. 

The flavor-symmetry approach has several positive features. It is a model-independent 
approach based on symmetries of the QCD lagrangian alone, as described in Section 12.21 
It also includes all the hadronic long-distance effects naturally, since they are contained in- 
side the hadronic parameters that are extracted from data. However, the major drawback 
is the inclusion of SU (3) breaking corrections, which are obviously completely unacces- 
sible within the framework. These symmetry breaking effects depend on the dynamics 
and computing them requires the use of non-perturbative methods. Alternatively, one 
might just assume that these effects are smaller than a certain percentage, and allow for 
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-0.2 0.2 0.4 0.6 



A dir (B° 3 k+k-) 

Figure 5.1: Comparison between the predictions for B s — > K + K~ from QCD factorization and 
flavor symmetry. The green box shows the QCDF prediction [26j (errors summed in quadrature), 
the gray regions show the predictions using B c i — ► 7r + 7r~ and flavor symmetry for different amounts 
of SU(3) breaking (c.f. Chapter 3]), and the blue horizontal line is the experimental value for the 
branching ratio [131] . 

arbitrary breaking of this size. However, even in the cases in which this can be done 
model-independently, the associated uncertainties that are generated preclude any preci- 
sion in the predictions. Moreover, the predictions are also affected by large uncertainties 
related to the fact that the used data has often large experimental errors. Therefore, the 
flavor-symmetry approach requires of two ingredients: more accurate experimental data, 
and more reliable estimates of SU(3) breaking. 

On the other hand, QCD factorization is a purely theoretical framework based on 
QCD and the heavy-quark expansion. In principle its predictions are unlimited up to a 
particular order in a s and AQ CD /m;,. This means that it can predict all the hadronic pa- 
rameters (not only those for which a flavor-symmetry relation can be found) provided the 
non-perturbative input of form factors, decay constants and LCDAs is given. However, as 
discussed at the end of Section 12.1.41 there are some difficulties with non-factorizable 1 / m& 
corrections that are numerically relevant because they are chirally enhanced. When includ- 
ing these corrections, the infrared divergencies must be subtracted in a model-dependent 
way, and the predictive power is partially lost. Moreover, as was already commented in 
Section 12.1.41 in QCDF the strong phases are perturbatively calculable, and therefore 
appear first at order a s , which means that they are in general predicted to be small. 
Phenomenologically, however, there are certain decays (for example B — > irir) in which 
the strong phases could be large, and not accountable by a 0(a s ) correction. The QCD 
factorization approach should then be provided with a prescription to control the phe- 
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nomenologically relevant power-suppressed infrared divergencies, and with a mechanism 
that generates large strong phases whenever they are there. 

In this chapter we present an approach that improves (at a phenomenological level) 
on the weak points mentioned above. This approach was proposed in [35] and since then 
it has been applied in [361 E3 EH1 [39] . 



5.1 The theoretical input: A 

We consider now the following quantity defined as the difference between tree and penguin 
contributions, 

A=T-P. (5.2) 

This quantity is a hadronic, process-dependent, intrinsically non-perturbative object, and 
thus difficult to compute theoretically. Such hadronic quantities are usually either ex- 
tracted from data or computed using some factorization-based approach, according to the 
long discussion given above. In the latter case, A could suffer from the aforementioned 
phenomenological problems of factorization. 

However, for a certain class of decays, T and P share the same long-distance dynamics: 
the difference comes from the (u or c) quark running in the loop, which is dominated 
by short-distance physics [35]. Indeed, in such decays, A = T — P is not affected by 
the breakdown of factorization that affects annihilation and hard-spectator contributions, 
and it can be computed in a well-controlled way leading to safer predictions and smaller 
uncertainties. 

Consider first the effect of the ^-operators in eq. (|2.44|) . The contribution to A will 
come from — T|- The current-current pieces will be just ~ ai(- • • ) + a 2 {- ■ • ) since they 
only appear in Tf. They will be proportional to the hard spectator functions H\ and H 2 
and they will be affected by the already familiar IR divergencies. However, the rest of 
the pieces will be ~ (af — af). So looking at eq. (|2.46p we see that they will just receive 
contributions from the penguin contributions, proportional to the functions (P" — P£). 
These functions are well behaved, and dominated by short-distance physics. At order 
a s , and neglecting electromagnetic corrections, penguin contractions are only present for 
i = 4,6, so the only contribution to A (besides the a\ and a 2 part) will be proportional 
to [PI - PI) ± r x (Pg - P 6 C ). Specifically, (see ref. [26]) we have 

P:(M 2 ) - P!(M 2 ) = -^± Cl [G M2 (0)-G M2 (ml/ml) 

P 6 «(M 2 ) - P 6 C (M 2 ) = Ci [Gm 2 (0) - G M2 {ml/ml)\ , (5.3) 

where G(x) and G{x) are certain penguin functions. The i?-operators in Eq. (|2.47p contain 
also the pieces proportional to b\ and 62 that enter and not T^, and these pieces will 
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fi af of Fq^ K m s (2GeV) m c /m b 

(m b /2,2m b ) 0.2 ± 0.2 0.1 ±0.3 0.34 ± 0.05 103 ± 20 MeV 0.30 ± 0.06 



Table 5.2: Inputs used in the computation of A d KK . 





H af 




m s (2GeV) 


m c /m b 


\A d 1 


19.7% 4.1% 0.9% 


15.8% 


5.7% 


53.7% 



Table 5.3: Relative impact of the input uncertainties on the error of A' 



introduce the IR divergencies from weak annihilation. The rest of the pieces, however, as 
can be seen looking at the specific form of the coefficients hi [26], cancel completely in the 
difference bf — b\. 

So now consider penguin-mediated decays that do not receive contributions from a\ t 2 
or b\$.- For these decays, the quantity A is at 0(a s ) free from dangerous IR divergencies 
in QCDF. This makes it a quantity that can be computed with acceptable accuracy and 
can be used as a reliable theoretical input. In general, A is given by 

A MlM2 = A MlM2 ^-C 1 [G M2 {ml/ 'ml) - G M , 2 (0)] , (5.4) 

where G M2 = G M , 2 ± r^ 2 G M2 , with the plus or the minus sign depending on whether Mi 
is a pseudoscalar or a vector meson. 

In our analysis of B —> KK modes, we are mainly interested in Aj (K corresponding to 
B(i — ► K°K°. A careful computation gives 

\A d KK | = (2.96 ± 0.97) x 10" 7 . (5.5) 

Only the absolute value is physical, because an unphysical global phase in the amplitude 
would rotate simultaneously P and T, rotating A accordingly. In particular, in the for- 
mulae for the hadronic parameters given in the following section, A can be chosen to be 
real, introducing only a global phase in the amplitude. From now on we will always take 
it as a real quantity. The values of the inputs used in the numerical evaluation of A^ K , 
and the impact of each one of them on the error, are given in tables 15.21 and 15.31 

5.2 Tree and Penguin contributions 

We now describe how to obtain the hadronic parameters T and P from A and experimental 
data. The derivation will be completely general, and is valid for any decay, independently 
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of whether the computation of A is theoretically clean or not. We begin writing the two 
self-conjugated amplitudes in terms of tree and penguin contributions, 



A = \ U D >T + A C D) *P , A = \ u d) T + a { 1 P ) P 
Now we put T = P — A and we square the amplitudes, 



(5.6) 



\A\< 



\A\< 



n (D)* X (L>)*, 2 
A c -+- A u 



n (D) .(D), 2 
A c -r A u 



P + 
P + 



.(D)* 



AD) 



■A 



(5.7) 



.(£>) ,(£>) 



A 



But the squared amplitudes are directly related to observables, 



\A\ 2 
\A\ 2 



BR{l + A dil )/g PS 
BR{1 - A dir )/g PS 



(5.8) 



where gps is the usual phase-space factor (c.f. Eq. (|6.5|) ). So two equations can be written 
that relate observables with P and A, 



BR(1 + A dil )/g PS 

I ~r Au | 

BR(1 - A dir )/g PS 



P + 



A u 



.(£»)* .(D)* 
Ac T Au 



P + 



A U 



X P) , >(£>) 
Ac "T Au 



A 



A 



(5.9) 



These are the equations for two circles in the complex P plane, whose solutions are the 
two points of intersection. This will result in a two-fold ambiguity in the determination of 
P and T (see Fig. I5.2p . Before writing down the analytical solutions, notice that in order 
for solutions to exist, the separation between the centers of these circles must be smaller 
than the sum of the radii but bigger than the difference. This translates into a consistency 
condition between BR, A d { r and A: 



< 



BR > 




(5.10) 



where BR = BR/gps , and IZd is a specific combination of CKM factors, TZd 

-o 



%c^\/Jc 2 D) (see below). The approximation for the upper bound on \A d \ r \ holds up 
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Figure 5.2: Eqs. |5J|) can be written as i?i = \P - Cf | and R 2 = \P - C^\, where the R t 
depend only on data and CKM elements, and the Cf^ depend only on A and CKM elements. 
These are the equations for two circles with centers in and radii Ri in the complex P plane. 
Their crossing points are the solutions for P. The condition for the existence of solutions is 
\Ri - #2 1 < |Cf -Cf\<\Ri+ R 2 \, which translates into Eqs. (j5riU)l . 



to very small corrections in the usual situation A < O(10~ 7 ) and BR ~ O(10 -6 ). The 
first condition in eq. (|5.10p turns out to be highly nontrivial. For example, Fig. 15.31 shows 
the allowed values for the direct CP asymmetry of B^ — ► K°K° in terms of its branch- 
ing ratio. It can be seen that for the present data of the branching ratio, the direct CP 
asymmetry is predicted to be less than about ~ 20%. 
The hadronic quantities P and T are then given by 



Im[P] 



BRA diT 
2c^A 



Re[p] = - C [ D) A± 



(D) A\ BR 



-im[py- v + 



(D) ' JO) 



T = P + A, (5.11) 
where the coefficients cf^ are again some specific combinations of CKM elements (see 



Table \52 

Equations (|5.1ip allow to extract the hadronic parameters T and P from experimental 
data on BR and Ada, information on sides of the unitarity triangle and the weak phase 7, 
and the theoretical value for A. This method is also powerful because if no experimental 
information is available for ^4diD one can j us t vary Adir over its allowed range in eq. (|5.10p . 
So in fact T and P can be extracted from BR, A and CKM elements. This is just 
the case for B d -» K°K°: while there is a preliminary experimental number for the 
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1 2 _ 3 4 

-» iT°ir ) x 10 6 



Figure 5.3: Allowed values for the direct CP asymmetry of — > if ^ in terms of its branching 
ratio, according to the value of £± d KK - The vertical blue band corresponds to the experimental 
branching ratio, which sets an upper bound on |^4di r | slightly above 0.2. 



c o 


M 
c i 


M 


n d 


-3.15 • i(r 5 


-0.034 


6.93 • 10~ 5 


7.58 • 10~ 3 


w 

c o 


c {s) 




n s 


3.11 • 1(T 5 


0.011 


1.63 • 10" 3 


1.54- 10" 3 



Table 5.4: Numerical values for the coefficients c- and TZd for 7 = 62°. 
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direct CP asymmetry, it is still too uncertain. Therefore, in our analysis we will assume 
that AdhiBd — ► K°K°) is unknown, and we will give results for different values of this 
asymmetry between —0.2 and 0.2, according to Fig. 15.31 We will also see that this range 
can be still made smaller by looking at data on Bj — > ir + ir~, but this will not be taken 
into account in our analysis. 

Finally, we would like to mention that this procedure is also useful to extract hadronic 
parameters from charged B decays, were the CP observables are just two: branching ratio 
and CP asymmetry. An example is the decay B + — ► ir + (j), for which the very small value 
~ 10~ 8 indicates a very small CP asymmetry, as long as the branching ratio is not 
too small (at present there is only an upper bound). 



5.3 Flavor symmetries and QCDF 

Using [/-spin symmetry, we can relate the two penguin-mediated decays — > KqKq and 
B s — > KoKq, as exemplified in fig. 15.41 Let us stress that we work with the operators of the 
effective Hamiltonian: internal loops have already been integrated out to yield four-quark 
operators, so that the internal loop of the u-penguin is not affected by [/-spin rotations, 
[/-spin breaking should be much smaller here than usual: it does not affect final-state 
interaction since both decays involve the same outgoing state, and it shows up mainly in 
power-suppressed effects. This is confirmed by QCDF: 



P 



sO 



fP 



(10 



1 + 



(A d KK /P d0 ){Sai + 5p c 3 + 25f3i}] = f P d0 (1 + 5 ps0 ) , 



T so = fT d0 [l + (A d KK /T d0 ){5a% + 5(3% + 25/?f}] = fT d0 (1 + 5 T s ) , (5.12) 

where we define the [/-spin breaking differences 5a? = a?— a? (id. for (3). The superscripts 
identify the channel and the bar denotes quantities for decays with a spectator s quark. 
Apart from the factorizable ratio : 



/ 



Ai 



Ad 

KK 



B S ^K 



M 2 Bd F»« 



(0) 



B d ^K 



(0) 



0.94 ±0.20 



(5.13) 



which should be computed on the lattice, [/-spin breaking arises through l/m^-suppressed 
contributions in which most long-distance contributions have cancelled out. 

First, the hard-spectator scattering (5a) probes the difference between B^- and B s - 
distribution amplitudes: 



5a\ = a s C F C^/N 2 c x 5\ B x [{x)% + r«{x) K X%] 
5X B = B s KK M B J{A s KK \ Bs )-B d KK M Bd /{A d KK \ Bi ) 



(5.14) 
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Figure 5.4: Diagrams contributing to (a) Bd 
J7-spin transformation. 



K°K° and (b) B s K°K° related through a 



B q KK = fBqfliGp I 'y/2, {x)k and Mb /An are first and first inverse moments of K and -Bq 
distribution amplitudes [2(1: . respectively. 5Xb is expected small, since the dynamics of 
the heavy-light meson in the limit nib — > oo should vary little from Bj and B s . Second, the 
annihilation contributions {5(3) contain a [/-spin breaking part when the gluon is emitted 
from the light quark in the .B^s-rrieson (this effect from A\ and A\ defined in [26] is 
neglected in the QCDF model for annihilation terms). 
Taking the hadronic parameters in [26], we obtain 



|<5 P .o| < 0.03 , |5 Ts o| < 0.03 . 



(5.15) 



KK 



fAj CK up to l/m^-suppressed corrections, 



These relations yield also the constraint Aj 
relating observables in both decays. 

Relations exist between B^ — > KqKq and B s — > K + K~ as well. A combination of 
[/-spin and isospin rotations leads from the penguin contribution in B^ — ► KqKq to that 
in B s — ► KqKq, then to B s — ► K + K~, up to electroweak corrections (it corresponds to 
fig. 15.41 up to replacing d — > u in the right-hand diagram). On the other hand, there are no 
such relations between tree contributions, since B s — > K + K~ contains tree contributions 
which have no counterpart in the penguin-mediated decay B^ — > KqKq. This is seen in 
QCDF as well: 

p s± = fp do^ + {A d j pd o ) ^ {§a c + m + 26p c^ = fP do {1 + 6ps±) ; (5 . 16) 

1 + ^-{h + 8a% + Sffi + 2^4}] = fT d0 + A s KK ai (1 + 8 T ,±) . 



From QCDF, we obtain the following bounds: 

\5 P s±\ < 0.02 , \S T s±\ < 0.04 



(5.17) 



The latter bound shows that the flavor-symmetry breaking corrections are smaller than 
T dG /{A d KK ai) = O(10%). Fortunately, T s± is strongly CKM suppressed in B s -» K + K~ 
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so that the uncertainty on its QCDF determination will affect the branching ratio and 
CP-asymmetries only marginally. 

Finally, these relations between and B s hadronic parameters are affected by elec- 
troweak penguins, small in the SM but potentially enhanced by NP effects. 

5.4 SM predictions for B s — > KK and SU(3) breaking 

In this section we present the results for the branching ratios and CP asymmetries of 
B s — ► K°K° and B s — > K + K~ . The knowledge of the B s — > K + K~ hadronic parameters, 
which require data from ^ -> i^if (and not S d — > 7r + 7r as in Chapter 0|) allows us 
to give also predictions for the SU(3) breaking parameters \r?\ and £ used in Chapter 0] 
relating B s -> K + K~ to £ d 

Let us outline briefly the basic steps in the process. 

• The Bd — > K°K° hadronic parameters P d0 and T d0 are obtained from Eqs. (|5,11|) . 
with the value of A c j XK given in Eq. (|5,5p and the experimental value for the B^ — > 
K°K° branching ratio [132] 

BR(B d ir K°) cxp = (1.08 ± 0.30) • 10^ 6 . (5.18) 

Since the B& — > K°K° direct CP asymmetry (A^) is still uncertain, we vary its 
value between —0.2 and 0.2 according to the discussion in Section [5.21 (see Fig. 15. 3p . 

• The hadronic parameters P s0 , T s0 , P s± and T s± for B s K°K° and B s K+K~ 
are obtained as described in Section [573| in particular from Eqs. (|5.12p . (|5.16p . (|5.15p 
and (|5.17p . In order to account for well-behaved short-distance l/m^-suppressed cor- 
rections not yet included, we stretch the SU (3) breaking bounds to more conservative 
ranges: 

|<5pso|, |<5t"°I; l<5p s ±l> |<5t s± I — 0-0^ • (5.19) 
This 8's are complex quantities and their phases are varied from to 27r. 

• The observables are then obtained from the ranges of the hadronic parameters. 

There is a further comment before presenting the results. We already mentioned in 
Section 15.21 that the extraction of the hadronic parameters results in a two- fold ambiguity 
(see Fig. 15. 2p . We lift this ambiguity on the basis that the rejected predictions would be 
in conflict with the Bj — > 7r + 7r _ analysis made in Chapter HI First, this solution requires a 
way too large [/-spin breaking in arg(P/T). Second, a prediction of the rejected solution 
is that A 1 ^. < 0, which is clearly in contradiction with the Bj — > ir + ir~ analysis, that 
predicts it positive (see Fig. \A.Qb . 
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BR S " x 10 6 


Af ir x !0 2 


A s ° ix x 10 2 


Kr = -0.2 


20.8 ± 6.7 ± 4.6 ± 1.6 


0.9 ±0.3 ±0.1 


-0.2 ±0.8 ±0.1 


Af ir = -0.1 


20.6 ± 7.5 ± 4.6 ± 1.6 


0.4 ±0.3 ±0.0 


-0.7 ±0.8 ±0.1 


A d d i = 


20.5 ± 7.5 ± 4.5 ± 1.5 


0.0 ±0.3 ±0.0 


-0.8 ±0.6 ±0.1 


Ajl = 0-1 


20.6 ±7.6 ±4.5 ±1.6 


-0.4 ±0.3 ±0.0 


-0.7 ±0.8 ±0.1 


Af ir =0.2 


20.8 ± 7.0 ± 4.6 ± 1.8 


-0.9 ±0.3 ±0.1 


-0.2 ±0.7 ±0.2 



Table 5.5: Observables for B s — > A A as functions of the direct asymmetry Aiir(-Bd — > A' A" ) 
within the SM. The errors for the branching ratios are split in: (1) Uncertainties of all inputs 
except / and 7, (2) Uncertainty of /, (3) Uncertainty of 7. For the asymmetries, the second error 
is due to 7. The preferred values for Af9 r arc highlighted (see the text). 

In Table [531 we present the results for the B s — > K°K° observables. An average over 
all the values for A^ r gives a single prediction for the branching ratio, 

BR{B S -» K°K°) = (21 ± 8 ± 5 ± 2) x 10" 6 , (5.20) 

and the following ranges for the CP asymmetries, 

- 0.011 < Af ir < 0.011 , (5.21) 
-0.015 < Af iix < 0.005 . (5.22) 

In Table 15.61 we show the results for the B s — > K + K~ observables. Again, a global 
average gives the following prediction for the branching ratio, 

BR(B S -» K + K~) = (24 ± 9 ± 5 ± 2) x 10~ 6 , (5.23) 

and the following ranges for the CP asymmetries, 

- 0.22 < A s d f r < 0.49 , (5.24) 
-0.55 < A s n f ix < 0.40 . (5.25) 

The information from the B^ — > tt + tt~ analysis can also be used to chose some preferred 
values for A^ r between those considered in Tables 15.51 and 15.61 Indeed, as can be seen 
from Table 15.61 there is a sign anticorrelation between A^ v and ^4^ x - The analysis in 
Chapter d] clearly predicts negative values for A s ^ ix , which correspond, roughly, to positive 
values of A^ T . All these correlations between all these B — > KK and B — > irn modes 
turn out to be extremely powerful. In the tables we have highlighted the preferred results 
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Si? s± x 10 6 


Alt x 10 2 


A S mix X 10 2 


Kr = -0-2 


24.9 ±9.3 ±5.0 ±2.3 


20.1 ± 15. 7± 2.6 


26. 2± 11.5 ± 1.9 


*fir = -0-1 


22.1 ± 8.5 ± 4.7 ± 1.7 


33.9 ±13.0 ±3.3 


8.2 ± 15.1 ±0.2 


^?r = 


20.1 ±7.3 ±4.4 ±1.2 


35.3 ±12.8 ±2.8 


-9.1 ±15.3 ±1.1 


^?r = 0-l 


18.4 ±6.8 ±4.2 ±0.9 


27.3 ±27.9 ±2.1 


-26.4 ± 18.4 ± 1.8 


Kr = 0-2 


17.5 ±6.5 ±4.1 ±0.9 


0.5 ±33.6 ±1.4 


-38.9 ±14.8 ±2.8 



Table 5.6: Observables for B s — > K + K~ as functions of the direct asymmetry A^^d — > K°K ) 
within the SM. The errors mean the same as in Table l5.5I Thc preferred values for Af? r are high- 
lighted (see the text). 

(corresponding to > 0). With these restrictions the global averages given above get 
tightened, for example, 

BR(B S -> K + K~) -> (20 ± 7 ± 4 ± 1) x 10~ 6 , (5.26) 

The prediction for the branching ratio in Eq. f|5.23j) or (|5.26p can be already compared 
with the experimental value |130| . 

BR{B S -> K + K~) cxp = (24.4 ± 1.4 ± 4.6) x 10~ 6 , (5.27) 

and fits nicely within the errors. Notice that the experimental errors are quite small. 
It is interesting to note that before this experimental number was released, the current 
experimental value was ~ 40 x 10~ 6 . In [35], the discrepancy between the theoretical 
results derived here and that high experimental number for the branching ratio led to the 
statement that [35] ". . . the data suggest a departure from the SM, to be further checked 
experimentally" , so these results were genuine predictions, and indeed, the 'further exper- 
imental check' fell right on top of the theoretical number. 

The predictions for the B s — » K + K~ observables can be compared to the other deter- 
minations (c.f. Fig. I5.ip . This is shown in Fig. 15.51 where it can be easily seen that the 
improvement is substantial, specially in the prediction of the branching ratio. 

Finally, we comment on how the results derived for the hadronic parameters can be 
used to predict the SU(3) breaking parameters \tt\ and £ connecting B^ — ► 7r + 7r _ to 
B s — ► K + K~ , introduced in Chapter HI The point is that, once the B& — > 7r + 7r~ hadronic 
parameters T*^ and are extracted from data as in Chapter HI and the B s — > K + K~ 
hadronic parameters T s=t and P Sziz are obtained from B^ — ► K°K° as described here, the 
SU (3) breaking factors that connect them can be obtained. 
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Figure 5.5: Comparison between the determinations of the B s — > K + K~ branching ratio and 
direct CP asymmetry from flavor, QCDF and the method described in this chapter (we refer to 
it as DMV). This plot is an extension of Fig. 15.11 and shows how the new predictions improve 
considerably the situation posed at the beginning of this chapter. 

To that end it is more useful to write the — > 7r + 7r - hadronic parameters obtained 
in Section [4.21 not as ranges, but as numbers with errors: 

\T%±\ = (5.48 ±0.42) x 10" 6 , 
\P^/T^\ = 0.13 ±0.05. (5.28) 

Then Eqs. (I4.26D give the values for the SU (3) breaking parameters. In Table [5771 we show 
the derived B s — > K + K~ hadronic parameters and the corresponding predictions for the 
SU(3) breaking parameters. 





\T S± \ x 10 6 


\P S± /T S± \ 


arg {P s± /T s± ) 


M 


e 


Kr = -0-2 


13.0 ±2.4 


0.10 ±0.04 


(38 ±28)° 


2.4 ±0.5 


0.74 ±0.40 


Af ir = -0.1 


12.3 ±2.6 


0.10 ±0.04 


(77 ±27)° 


2.2 ±0.5 


0.78 ±0.41 


AZ = 


11.7±2.6 


0.11 ±0.04 


(105 ±25)° 


2.1 ±0.5 


0.82 ±0.46 


AZ = M 


11.2 ±2.6 


0.11 ±0.05 


(134 ±38)° 


2.0 ±0.5 


0.86 ±0.51 


AZ = 0-2 


10.8 ±2.4 


0.12 ±0.06 


(179 ±30)° 


2.0 ±0.5 


0.90 ±0.56 



Table 5.7: Hadronic parameters for B s — > K + K and [/-spin breaking parameters |rx| = 
\T S± /T^\ and ^ = \P S± /T S± \/\P£/T£\, that relate B s -» K+K~ and B d -» it+tt-. 
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Penguin-Mediated B —> VV Decays 
and the B s — B s Mixing Angle 

The method presented in the previous chapter is powerful for many reasons. One of the 
reasons is that the theoretical input (A) is quite clean, and it is the minimal information 
that one needs to obtain the hadronic parameters T and P without any information on the 
mixing-induced CP asymmetry of neutral B decays. This property was used in Chapter 
[5] because there is still no experimental information on A^^B^ — ► K°K°). However, 
because this asymmetry is the one that contains the information on the mixing angle 
(c.f. eq. (|3.69p ). this method can be used to extract the mixing angle from data, without 
the necessity of neglecting a part of the amplitude (see Section [3.4.2p . or doing a 577(3) 
analysis with some other modes, which probably requires the input of other mixing angles. 

In this chapter we use the method developed in Chapter [5] to provide three alternative 
strategies to extract the B s — B s mixing angle (fi s from non-leptonic B decays into pairs 
of vector mesons. In particular we focus on B^ g — > K*° K*° , (pK*° , <jxf> decays, for which 
the extraction of A is clean, according to the discussion in Section 15.11 

However, the theoretical study of B decays with vector-vector final states is compli- 
cated. The main reason is the fact that transverse amplitudes do not factorize at leading 
order. In fact, naively one expects the transverse amplitudes to be suppressed in the heavy 
quark limit with respect to the longitudinal one. This works quite well, for example, for 
B — ► pp, where the longitudinal polarization fraction is ~ 95%. However, it does not work 
for some penguin mediated decays such as B — ► (pK*°, where the measured polarization 
fraction is ~ 50%. This has been called the polarization puzzle. It is therefore important 
to estimate properly these contributions when used; but because theoretically this is a 
challenge, we will concentrate on longitudinal observables alone. 
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6.1 Longitudinal observables in B —> VV modes 



The amplitude for a B meson decaying into 2 vector mesons can be written as 



A(B -> ViV 2 ) 



— 4— ( e i -Pb)^2 -PB) 



B 



-(e* 
2 {€l 
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(PB ■ e*)(ps • 4) 



+ 



{ * 



(PB-e*)(pB 



+ 



m. 



%e lLvp<j^\ ^2 PlP2 

2pi • P2 



A, 



A. 



(6.1) 



where Aq + correspond to the amplitudes for longitudinal and transversely polarized final 
vector mesons. It is also customary to use the basis -Ao,||,±) where A»± = (A+ ± A-)y/2. 

The vector mesons in the final state decay typically into pairs of pseudoscalar particles. 
A full angular analysis of vector-vector modes provides the following set of observables: 
three polarization fractions /n, f± and fu (only two of them are independent) and their 
CP-conjugate counterparts ib,_i_,||> two phases <p±,\\ (again, together with (j)± : \\), a total CP- 
averaged branching ratio BR, and a total direct CP-asymmetry A<±\ r . The polarization 
fractions are defined as 



/o,. 



I An 



W + \A A 



+ \A 



k. 



|A | 2 + |Aj 



+ All 



(6.2) 



A full angular analysis is available for Bj — > (f>K*° from BaBar and Belle [118|. I121j , and 
the same type of analysis is expected for B& — > K*°K*°. 

We will focus in this chapter on observables for the longitudinal polarization (£>i? long , 
Ad°ir g ; A^fx an< ^ A^p g ), where only A occurs. These observables, free from the positive 
and negative helicity components, can be predicted with a much better accuracy. Indeed 
the negative-helicity (positive-helicity) component of the amplitude is l/m^-suppressed 
(1/m^-suppressed) because of the nature of the interactions involved (left-handed weak 



interaction, helicity-conserving strong interaction at high energies) \29\ I27j. This suppres- 
sion makes longitudinal observables better behaved and easier to compute than transverse 
ones. 

Some decay channels exhibit the 1 /m^-suppression of transverse amplitudes in a very 
striking way : the longitudinal polarization is very close to 1, e.g. fi — 97% for B — > 
p + p~ . In such cases, the full observables (where Aq is replaced by the sum A = Aq + 
A_ + A+.) coincide with the longitudinal ones to a high degree of accuracy. On the 
other hand, for penguin dominated AS* = 1 decays, fi can be as low as ~ 50%, so that 
the transverse amplitudes (or ± helicity amplitudes) contribute significantly to the full 
observables. Therefore, one must determine whether purely longitudinal observables can 
be extracted from experimental measurements. 
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Figure 6.1: Definition of the various angles in the decay B — > V\{— > Pi-P2)V2(— > Q1Q2). 

We start from the normalized partial decay rate of B — > V1V2, where the two vector 
mesons go subsequently into pairs of pseudoscalar mesons. It can be written [133] 

d?T 9 1 

Fd cos 6\d cos e 2 d(j) 8n \A \ 2 + |yl||| 2 + \ A ± \ 2 

\A \ 2 cos 2 61 cos 2 9 2 + |^||| 2 ^sin 2 61 sin 2 9 2 cos 2 (j> 

+ \Aj_\ - sin 2 0± sin 2 6 2 sin 2 <f> + Re [AJ A|| ] — ^ sin 26»i sin 209 cos d> 
2 2y 2 



+ ImLAgAJ — -= sin 20i sin 20 2 sin + Im[A|j — sin 2 X sin 2 2 sin 2</> 
2y 2 2 



(6.3) 



where (0i , 62, <j>) are angles introduced to describe the kinematics of the decay B — > V± V2 
followed by V\ — > P1P2 and V2 — > Q\Q2- 0\ is the angle of one of the V\ decay products in 
the rest frame of V\ relative to the motion of V\ in the rest frame of the I?-meson (same 
for 02 with V2). 4> is the angle between the two planes formed by the decay products of V\ 
and V2 respectively (see Fig. 16.11 for a representation of the angles) . 

There are different ways to perform the angular integrations in order to extract the 
purely longitudinal component from the differential decay rate. A first option consists in 
computing moments of cos0i (or equivalently cos 02 ) '■ 



r '° ng = / a 7TT- a a a (l cos2 ^ ~ I) dcos0 1 dco S 02d<j ) = g P s\A o \ 2 /T B (6.4) 

J dcos0\dcos02d(p \2 2 J 



'id cos V2 C ° 1 2, 
where gps is the product of phase-space and lifetime factors 



9PS = J^ps\/l M B ~ (™i + m 2 )2][M 2 - (mi -m 2 ) 2 ] (6.5) 

B 
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A second possibility amounts to performing asymmetric integrations over one angle [134j 

r-l r r-2n J3p 

= / dc os0i / dcos9 2 / #- — =g PS \A \ 2 /r B (6.6) 



i 



T J d cos ^idcos 02i 



with 

11 r 71 "/ 3 5 /•27r/3 ii 
dcos0 2 = I — / -- / +— / I (-sme 2 )de 2 (6.7) 



T \ y ./0 y Jtt/3 9 ./27r/3, 

In the same way we can obtain the CP-conjugate T lons (B q — > f) from the corresponding 
CP-conjugate distribution, leading to the CP-averaged branching ratio of the longitudinal 
component 

I A I 2 I A 1 2 

Belong = ZE. (r^(B° q -, /) + r lon S(S g ° - /")) = g P5 ' o| +' o| (6.8) 

where j4o is the CP-conjugate amplitude of Aq. 

If we include the dependence on time in the above expressions, B-B mixing modifies 
the expressions. We will focus on CP-eigenstates / CP in the final state K*°K*° and <fi<fi, 
as well as cpK* with a subsequent decay of K*° into a CP-eigenstate (K s ir° or K^ir ). 

The time evolution of these observables is obtained by considering the time dependence 
of .Ao(t) |135j . Inserting this time dependence one arrives at the usual expression for the 
longitudinal component of the time-dependent CP-asymmetry: 

.long,. r^(B° q (t) -> f CP ) - T^(B° q (t) -> f CP ) = A 1 ^ cos (AMt) + A 1 ^ sin (AMt) 

CP I ) r iong (jB (t) ^ /cp) + r iong(B0(t) - / cp ) cosh (AR/2) _ ^ng ^ (AR/2) 

(6.9) 

where the direct and mixing-induced CP asymmetries are defined by: 



long o Re(e-^M^ ) 
Ar " l^ol 2 + lAol 2 



d long _ 1^01 ~ 1^01 A long _ 9 ^Hg r '"^Q^0j / R , m 

- |Ao| 2 + |A | 2 ' mix ~ lA^ + IA)! 2 1 j 

together with the asymmetry related to the width difference : 

(6.11) 

(f>M is the mixing angle and Ar = T H — T L . ??/ is the CP eigenvalue of the final state / 
(±1): r] K *ox*o = rj<j}<j) = 1) whereas rj K ,o^ = 1 if K*° decays into K s ir° and —1 if it decays 
into Kltt°. In the latter case, the contribution from the strong process K*° — > Kir is 
the same for both B and B decays and it cancels in the time-dependent CP-asymmetry 
Eq. (|6.10p . which depends only on the amplitudes Aq and Aq. 

Finally, if the direct CP-asymmetries of all three helicity components are negligible, 
the longitudinal branching ratio can be estimated very easily from: BR lon z = BR total f . 



6.2 Inputs 
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BR long (B d -> K*°K*°) x 10 6 



Figure 6.2: Allowed region on the BR lon s - A l ^ s plane for B d -> K*°K*°, according to the value 
of A£, K ,. 

6.2 Inputs 

In this section we compute the values of the quantities A for the decays B djS — ► if* A'* , 
-Brf.s - ► <pK*° and -B s — ► 00. Because the final states can appear in three different polar- 
ization states, a different A should be associated to each polarization. However, since we 
are only interested in longitudinal polarizations, we only give the results for this case. 
For B dtS -> K*°K*° we obtain, 

|A£* K *| = A d £ K . x^Ci \G K * (8 e ) -G K * (0)| = (1.85 ± 0.79) x lO" 7 GeV (6.12) 

47TiV c 

|AJ^.| = Af.^^^d \G K *(s c ) - G K *(0)\ = (1.62 ±0.69) x 1(T 7 GeV (6.13) 

47TiV c 

where Gy = Gy—r^Gy are the usual penguin functions and ^viV are ^ ne n& i ve factoriza- 
tion factors combining decay constants and form factors (see Sections 12.1.21 and 15. ip . The 
numerical values of the used inputs are given in Table 16. 1\ which are taken from Ref . [27] . 
The contributions to each error from the various sources are detailed in Table 16.21 

The conditions in eqs. (|5.10p derived in Chapter [5] can be applied to the B d -> K*°K*° 
longitudinal branching ratio and CP asymmetry given the value for A^., x * in eq. f|6. 12j) . 
The allowed region in the observable space is shown in Fig. 16.21 

For the A, as well as for the other quantities computed in this chapter (as well as in 
Chapter [5]) , we quote as the central value the value obtained from taking the central value 
of the inputs. To estimate the error, we vary one by one each of the inputs, compute the 
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m c (m b ) 


f B 




^B ; Ag s 


ai(K*) 


1.3 ±0.2 


0.21 ±0.02 


0.24 ±0.02 


0.35 ±0.15 


0.06 ±0.06 






fx- 


/^.(2GeV) 


aB^k' 




0.1 ±0.2 


0.218 ±0.004 


0.175 ±0.025 


0.39 ±0.06 


0.33 ±0.05 



U 


#(2GeV) 






0.221 ±0.003 


0.175 ±0.025 


0.38i°:J° 


0.0 ±0.3 



Table 6.1: Input parameters required in QCD factorisation to compute the quantities A's and S's 
described in the text. The masses and decay constants are given in GeV. 

difference with the central value, then add in quadrature the resulting uncertainties. The 
main sources of uncertainties are the scale of factorization fi, the mass of the charm quark 
m c , and the form factor Aq^ k *. 

In a similar way, we can compute the corresponding longitudinal A for the decay modes 
B d:S -> <j>K*° and B s -> 



= (1.02 ± 1.11) x 1(T 7 GeV (6.14) 
= (1.16 ± 1.05) x 10~ 7 GeV (6.15) 
(2.06 ± 2.24) x 10" 7 GeV (6.16) 



In the following Sections we show how to apply the results of this section and Chapter 
[5]to the longitudinal contribution of penguin-dominated B — ► VV modes. We will see that 
they can be used to extract the B s — B s mixing angle and some longitudinal observables 
like branching ratios and time-dependent CP asymmetries within the Standard Model. 
In particular, we outline three different strategies to determine the B s — B s mixing angle 
(in the SM and beyond). Indeed, concerning New Physics we will see that under the 
assumption of no significant New Physics affecting the amplitude, while Strategy II can 
detect the presence of New Physics by comparing the obtained 4> s with 4> S S M = 2(3 S , 
Strategy I and III can not only detect New Physics but allow also for the extraction of 4> s 
even in the presence of New Physics in the mixing. 



*+k>\ = A^^C.IG^-G^m- 
\ A U\ = ^^C 1 \G,(s c )-G,(0)\ = 



6.3 First strategy to extract (j) s : Bounding T/P 
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m c 


^0 


fx* 


M 


ai(K*) 


a 2 (K*) 


ai(K*) 




Ad 

K* K" 


37.3% 


13.2% 


0.2% 


44.2% 


0.1% 


4.6% 


0.1% 


0.3% 




37.5% 


12.9% 


0.2% 


44.4% 


0.1% 


4.7% 


0.1% 


0.3% 





m c 


AB^K* 
^0 


Ik- 


/^.(2GeV) 


a 


a x {K*) 


Ad 


44.2% 


2.0% 






52.3% 




A s 


35.0% 




0.1% 


0.7% 


58.2% 


0.7% 


A s 


44.1% 








52.3% 







a 2 (A'*) 




,5^0 


#(2GeV) 


a 2 (0) 




Ad 








0.4% 


0.7% 


0.3% 




0.1% 


0.1% 


5.0% 


0% 


0% 


0% 


A s 






2.1% 


0.4% 


0.7% 


0.3% 



Table 6.2: Relative contributions from the inputs to the errors in A for the various decays. 



6.3 First strategy to extract (j) s : Bounding T/P 



The b — > s penguin-dominated decays like B s — ► are in principle clean modes to 

extract the mixing angle (f> s . In this section and those following, <p s refers to the same 
mixing angle that will be measured, for instance, in the mixing induced CP asymmetry 
of B s — > ifj(j) including possible New Physics contributions in the mixing. When focusing 
only on SM we will use the notation 4> s = 2(3 S . 

In an expansion in powers of A« '/Ac , the amplitude for the decay B s — > K*°K*° is 
given by: 



A l Z!(B s - #* u iT u ) ~ sin</> s + 2 



*0 T^*0\ 



A 



Re 



K*K* 



P, 



K*K* 



sin 7 cos i 



(6.17) 



In order to determine the accuracy of this relation, we must assess the size of the CKM- 
suppressed hadronic contribution T. Notice that this relation is valid even in presence of 
New Physics in the mixing. In the SM, one can derive from the Wolfenstein parametri- 
sation that eq. (|6.17p is of order A 2 (with A = V us ), and both pieces shown on the 
r.h.s of eq. (|6.17p are of this same order. However, despite the smallness of the ratio 

(s) (s) 

\X£>/\?>\ = 0.044, a significant value of the hadronic ratio Re(T/P) could spoil the po- 



108 Penguin- Mediated B — > VV Decays and the B s — B s Mixing Angle 



tentially safe extraction of sin0 s . The deviation from sin</> s is: 
AS(B S -► K*°K*°) = 2 



Re ( K s * K * ) sin 7 cos S (6.18) 

™K*K* ■ 



We want to set bounds on Re(T/P), which can be related to the inputs: 
Re ( 1) = b. = t + Re fA\ 1 + Re(F) I A 



P7 \ P J \P J Re(P) 2 + Im(P) 

Eqs. (|5.1ip show that the maximum of Re(T/P) is reached for »4.^ g = together with 
the positive branch for Re(P). The following bound is obtained 

Re (^) < 1 + (" c i S) + a/-( c o S) M S) ) 2 + (V4 S) )^P/A 2 ) 1 (6.20) 



where the lower bound for PP long and the upper bound for A must be used. In a similar 

along 

dir 



way, the minimum of Re(T/P) occurs for A^^ e = 0, for the negative branch of Eq. (|5,lip 



for the solution of Re(P) 

Re (p) > 1 + (-cS s) - \/-(4 S) M S) ) 2 + (l/4 S) )^/A 2 ) _1 (6.21) 

where the lower bound for BR lon& and the upper bound for A must be used once again. 
As a conclusion, we obtain a range for Re(T/P) from two inputs: the branching ratio 
Bi? long (B s K*°K*°) and A^.*., given in Eq. lRH5j) . 

Using Eq. (|6. 18|> . these upper and lower bounds on Re(T/P) are converted into a 
bound on the pollution AS(B S -> i^*°K*°). The latter is plotted as a function of the 
longitudinal BR lon ^(B s -> K*°K*°) in FigESJ 

Once a measurement of «4™jlf (P s — ► iC* .R"* ) is available, upper and lower bounds for 
(p s are easily obtained. For instance, if we take as a lower bound for the branching ratio 
BPj on z(B s -> K*°K*°) > 5 x 1(T 6 , Fig. Ogives 0.03 < AS{B S -» K*°K*°) < 0.06. In 
the case of a moderately large branching ratio BR lon z(B s -> iT* ^* ) ~ (30 - 40) x 10~ 6 , 
the bounds get sharper, with 0.04 < A5(P S -> iP ^* ) < 0.05 and 

(A l °^(B s ^K*°K*°)-0M) <sin<^< (A^(B S ^ K*°K*°)-0M) (6.22) 

The same strategy can be applied to B s — > (j)K*° and B s — ► cjxf> decays 

• Take the experimental value for the longitudinal branching ratio BR long (once avail- 
able), and the theoretical value for A from Eqs. (16.150 or (I6.16p . 



6.4 Second strategy : Measuring CP asymmetries and BR 
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0.055 

* 

° 0.05 
| 0.045 
S. 0.04 
< 0.035 
0.03 




10 20 30 40 50 

BR lon z(B s -> K*°K*°) x 10 6 



Figure 6.3: Absolute bounds on AS(B S — > K *°K *°) as a function of the longitudinal branching- 
ratio BR lon z(B s -» K*°K*°). 



• Apply Eqs. (f6T20|) and ([6T21]) to constrain the range of Re(T/P). 

• Derive the allowed range for AS according to the equivalent of (|6.18|) 

• From the measured value of -4 m ^f , determine <fi s from 

(A^-AS max ) <sin0 s < (A l ^-AS min ) (6.23) 

A weak mixing angle (f) s different from (p^ M would signal the presence of New Physics. 

Interestingly, if the longitudinal direct CP asymmetry becomes available and happens 
to be inconsistent with zero, the bounds for Re(T/P) in Eq. (|6.20p and (|6.2ip can be 
tightened. Eq. (|6.19p can be exploited to derive expressions similar to Eq. (|6.2U|) and ([6.2ip 
with a non-vanishing ij™ g , leading to stronger bounds on Re(T/P) and consequently on 
sin 4> s . 

6.4 Second strategy : Measuring CP asymmetries and BR 

In this section, we show how we can extract mixing angles and related CKM phases in a 
clean way from experimental data, the length of two sides of the unitarity triangle and the 
theoretical quantity A. The only theoretical requirement is that the decay must allow for 
a safe way of computing A. The approach is general in the same sense as in the previous 
section, since it can be applied to any B decay into two pseudoscalars or vectors. But it 
yields different results for the four groups of decays: 
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1. B<i decay through a 6 — > <i process, e.g. i?^ 

2. B s decay through a 6 — > s process, e.g. _B S - 



3. Bd decay through a 6 
a CP eigenstate) 

4. i? s decay through a 6 
a CP eigenstate) 



s process, e.g. Bj 
d process, e.g. B s 



-> (j)K*° (with a subsequent decay into 
■> <pK*° (with a subsequent decay into 



As far as weak interactions are concerned, the difference between B^ and B s decays consists 
in the mixing angle, whereas b — > d and b — > s processes differ through the CKM elements 
^u?c , where D = d or s. 

In the case of a B^ meson decaying through a b — > D process (D = d,s), we can extract 
the angles a [57j and j3 from the identities: 



sin 2 a 



sin 2 /3 



BR 



2|A^P| 
BR 



A| 



1 



1-(A 



dir J 



(-^mix) 



2IA 



IAI 



dir J 



(A 



(6.24) 
(6.25) 



In the case of a B s meson decaying through a b — > D process (D = d,s), we can extract 
the angles (5 S [136] and 7, assuming no New Physics in the decay, from the following 
expressions: 



sin 2 p s 



sin 2 (f3 s + 7) 



BR 



2|Ai D) | 2 |A| 2 

BR 
2|a1 D) | 2 |A| 2 



1- A/1" (Alir) 2 - (Anix) 2 



1- a/1" (^dir) 2 " (^mix) 2 



(6.26) 
(6.27) 



If the obtained (3 S differs from its SM value, this would signal the presence of New Physics. 
Notice that this strategy is obtained by combining the definition of A with the unitarity 
of the CKM matrix, so it is designed to work only in the context of the SM. Consequently 
the previous expressions should be understood as a way of testing the SM. This is an 
important difference with Strategies I and III where one can obtain a value for the weak 
mixing phase also in the presence of New Physics in the mixing (but not in the decay) . 

While the previous equations are quite general (they can be used for B — > PP de- 
cays), it is understood that BR and Ajirmix refer to the longitudinal branching ratio and 
longitudinal CP-asymmetries, respectively, when they are applied to B — > VV decays. 

Eq. (|6.26p provides a new way to perform a consistency test for the SM value of | sin/9 s | 



from the measurements of A ^ f(B s — ► K* U K 



K*°K*°) and BR lon z(B< 
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K*°K ). The same strategy can be applied to B s — > (pK*° and B s — > cfxj) using the 
corresponding sum rules. This sum rule offers several advantages : it is independent of 
CKM angles, and all the hadronic input is concentrated on a single well-controlled quantity 
A. 

Note that all these equations depend actually on the corresponding branching ratio 
and *4^p g . The asymmetry A 1 ^? is indeed related to the direct and mixing-induced CP- 
asymmetries through the equality (^4^ s ) + {A l £^) + (A 1 ^?) = 1. It was already noticed 
in [127] in the context of B s -> K+K~ and in [W] in the context of B -> J/rpK* , D* + D* 
decays that it is possible to extract A 1 ^? directly from the "untagged" rate: 

T lons (B s (t) VV) + T lons (B~ s (t) -» VV) oc R^* 1 + R^e^ 1 (6.28) 
If the time dependence of both exponentials can be separated, one obtains 

A l Z & {B s - VV) = (6-29) 

The branching ratio and -4.^p g are thus the only required observables to extract (3 S through 
this method, which offers the advantage of concentrating in A all the hadronic input needed 
to bound the tree-to-penguin ratio. 

6.5 Third strategy : Relating B s -> K*°K*° and B d -> K*°K*° 

Once an angular analysis of B d — > K*°K*° is performed, it is possible to extract the CP- 
averaged branching ratio corresponding to the longitudinal helicity final state. Eqs. (|5.1ip 
can be used to extract the hadronic parameters, if one assumes that no New Physics 
contributes in an appreciable way. If flavor symmetries are sufficiently accurate for this 
particular process, this estimate can be converted into a fairly precise determination of 
hadronic parameters for the b — > s channel B s — > K*°K*°. We have seen in Chapter [5] 
that for Bd t s - > KK modes a [/-spin analysis combined with QCD factorization leads to 
tight constraints on the ratio of the tree contributions to both decay modes, as well as that 
for the penguins. In this section we show how to relate Bd — > K*°K*° and B s — ► K*°K*° 
decay modes following the same approach. 

We define the parameters 5^* K * and 5^* K , as 

Plc-K- = fP&* K *(l + €* K *), T> K . K . = fT*. K .(l + (6.30) 

where the factor / is given by 



ml A*^ K * 



m%A$- 



>K* 



0.88 ±0.19 (6.31) 
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~ 2 Rc(P? { f K ,) (xlO- 6 ) 



Figure 6.4: Values for P^* K , corresponding to a real Ak*k*, as a function of the longitudinal 
observablcs BR°(B d -> A"* /?* ) and Adir°(B d -> K*°K*°). The contours from black to gray 
correspond to branching ratios from 10~ 6 to 1CP 5 in steps of 2.5 x 10~ 6 . The contours from red 
to yellow correspond to values of the CP asymmetry from to its maximum allowed value in 
steps of 30% of its maximum value. The contours from red to pink are the same for negative CP 
asymmetry. 



We compute |<5 A ,'. K *| using QCDF. These parameters are affected by the model depen- 
dent treatment of annihilation and spectator-scattering contributions, so the results should 
be considered as an estimate. A significant part of long-distance dynamics is common to 
both decays, and we find the following upper bounds 

\&k*k*\ < °- 12 , \5 T k*kA < 0-15 (6.32) 

where the largest contribution comes from the lower value of Xb- 

We could in principle apply the same strategy to B^ s — ► <pK*°, but the corresponding 
5's are much larger. Indeed, the computation leads to corrections up to 8^* ~ 50%. This 
shows that [/-spin symmetry cannot be expected to hold at a high accuracy for any pair 
of flavour-related processes. 

#(*)#(*) ff er a mucri 

more interesting potential than other 
final states such as (j)K*°. Moreover, we cannot perform a similar analysis for (jxfi since 
Bd ~* ^ is a pure weak-annihilation process, contrary to B s — > (fxp mediated through 
penguins. Therefore we focus on the precise B s — ► K*°K*° modes in the remaining 
part of this section. Notice that the large hadronic uncertainties affecting B s — > (fxfr and 
B s — > (pK*° have no impact when we use these modes in the strategies described in Sections 
16.31 and 16.41 since we exploited a quantity A where they cancel out. 

Once the hadronic parameters P^*^* and T^, K , have been obtained from Eq. (|6.30|) . 
one can give predictions for the B s — * K*°K*° observables. Note that the branching 
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BR long (Bd -» K*°K*°) x 10' 



Figure 6.5: Longitudinal branching ratio for B s — > X* ^* in terms of the longitudinal — ► 
if* ^* branching ratio. The light-shaded area corresponds to the uncertainty on the ratio of 
form factors /, whereas the dark-shaded area comes from varying the various hadronic inputs. 



ratio BR long (B ( i — > K*°K*°) is an experimental input in this analysis, and this piece of 
information is not available yet. The result for the branching ratio of B s — > K*°K*° 
is given in terms of the Bd — ► K*°K*° branching ratio in Fig l6.51 Once the branching 
ratio of B d — ► [ s measured one can use this plot to find the SM prediction for 

BR long (B s K*°K*°). 

The ratio of branching ratios BR long (B s -► K*°K* )/ BR long (B d -» K*°K*°) and the 
asymmetries turn out to be quite insensitive to the exact value of BR long (B ( i — > K*°K*°) 
as long as BR} ong {B ( i — > > 5 x 10 -7 . The numerical values are summarised in 

Table E3 

Under the standard assumption that New Physics contribution to b —> d penguins is 
negligible, and since the experimental input comes entirely from B d — ► K*°K*° (a b — > d 
penguin), the results given in Table [6751 are SM predictions. In presence of New Physics 
in b — > s penguins the full prediction can be obtained by adding to the SM piece extra 
contributions to the amplitude and weak mixing angle as explained in [125} 1100], [138j. 

One may also use this as a strategy to extract the mixing angle <f) s . If one assumes no 
New Physics in the decay B s -> K*°K*°, this method relates directly A^(B S -> K*°K*°) 
and0 s . Fig l6.6l shows 

^Zx ( B * K*°K*°) vs. 4> s . Once this asymmetry is measured, this 
plot can be used as a way to extract (fr s , and this result can be compared to the one found 
in tree decays such as B — > DK. A disagreement would point out New Physics. Moreover, 
it is possible to distinguish whether New Physics affects the decay or the mixing itself, 
by confronting BR long (B s K*°K*°) and A l ^ g (B s -» K*°K*°) with the SM predictions 
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f BR lon z(B s -> K*°K*°y 


) =17 ±6 

SM 


\BR lo ^(B d -> K*°K*°) / 


A l ™ s (B s - K*°K*°) SM = 


0.000 ± 0.014 


A^(B S - ir*°^*°) SM = 


0.004 ±0.018 



Table 6.3: Results for the longitudinal obscrvables related to B s — > if* A"* according to Sec. 16.51 
These are predictions for the SM contributions under the standard assumption of no New Physics in 
b -> d transition. We used 0f M = 2/3 s = -2° for A^ , and we assumed BR lon s(B d -> K*°K*°) > 
5 x 10~ 7 . The quoted uncertainty includes the errors associated to all input parameters including 
the variation of 7 inside the range 56° < 7 < 68° [73| [74] . 

given in Table [6731 If the predictions for the branching ratio and the direct CP asymmetry 
agree with experiment, but the <f) s extracted from A^^(B S —> x*°K*°) differs from <f)f M , 
this will be a clear indication of New Physics in B s — B s mixing. An interesting comparison 
will be allowed between the value for <j) s obtained here and the measurement of <j) s from 
the mixing induced CP-asymmetry of B s — > DK decay |139[ 1140] . 

6.6 Discussion 

The increasing list of measured non-leptonic two-body Bj- and -E^-decays provides many 
tests of the CKM mechanism of CP violation in the Standard Model. Of particular interest 
is the determination of angles through time-dependent CP-asymmetries. For instance (j) s , 
related to B s — B s mixing, should be constrained : it is tiny in the Standard Model, and can 
be measured through many penguin-dominated decays. However, for such determination 
to be valid, one must assess the size of the various hadronic quantities involved as precisely 
as possible. 

In this chapter we have applied ideas presented in Chapter [5] for B4 s — ► KK to 
vector- vector modes mediated through penguins : B^^ — > K*°K*°, (pep and cpK*° (with 
the condition that K*° decays into a definite-CP eigenstate). In order to combine flavor 
symmetries with QCD factorization, we have restricted our analysis to longitudinal ob- 
servables, which are under better theoretical control. These observables have been related 
to the angular analysis performed experimentally in Sec. 16.11 Penguin-mediated modes of- 
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Figure 6.6: Relation between A l ^{B s -> K*°K*°) and the B s - B s mixing angle S . We 
assumed BR lons (B d — > K*°K*°) > 5 x 10~ 7 and 7 = 62°. A measurement of this asymmetry leads 
to a prediction for <fi s , which includes hadronic pollution and SU(3) breaking effects, according to 
Scc.Id31 





Strategy 1 


Strategy 2 


Strategy 3 


Inputs 


BR lons (B s -> K* a K*°) 
A S K , K «, 7 


BR lons (B s -> K*°K*°) 
A^iBs -> l-r°A* ) 
-^Cf ( B ° -» A'* A* ) 


BR lons (B d -> K*°K*°) 
Ax-.jf., <5t, <Sp, 7 


Outputs 


0s 


1 sin/3 s |, 7 


BR lons (B s -> A* A* ) SJU 
•4i™ g ( B « A*°A*°) SJf 

0s 


Advantages 


Applies also to 
B s -> 4>K*° and 
B s -* 00 


Applies also to 
B s -> 0A'*° and 
S s -► 00 


It can be easily gen- 
eralized to include 
New Physics in the 
decay and mixing. 


Limitations 


It assumes no New 
Physics in b — > s de- 
cay. 


It assumes no New 
Physics in b — > s de- 
cay. 


Does not apply to 
B s 0if*° or 
£? s — > 00 because 
5t,p are big. 



Table 6.4: Comparison between the three strategies for B s — ► AT* AT* . 
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fer the very interesting feature that the difference between tree and penguin contributions 
A = T — P should be dominated by short-distance physics. It can be computed fairly 
accurately using QCD factorization, and it can be used to determine tree and penguin 
contributions from observables as explained in Chapter This theoretical piece of infor- 
mation is used to relate CP-asymmetries of B ds — > K*°K*°, (jxf> and (j)K*° to CKM angles 
according to different strategies. For illustration, we have focused on Bd lS — ► K*°K*°, 
where all three strategies apply. 

In Sec. I6.3t we have proposed to use A = T—P to put stringent bounds on the pollution 
due to hadronic uncertainties. Indeed, even though the ratio |A« /Ac | = 0.044 is small, 
a large value of the hadronic quantity Re(T/P) could spoil the naively safe extraction of 
sin0 s from the mixed asymmetry of B s — > K*°K*°. This strategy to control the pollution 
can be applied to all penguin-mediated processes of interest here. A similar analysis for 
sin 2(3 from B d — > 4>Ks can be found in [39] . 

In Sec. EH we have suggested a second approach, using A^f(B s — ► K*°K*°), A l ^ g (B s 
K*°K*°) and BR long (B s — ► K*°K*°) to extract |sin/3 s |. In principle, one can also use an 
alternative set of experimental quantities : the branching ratio together with a direct mea- 
surement of the longitudinal untagged rate. The sum rule needed for the B s — > K*°K*° 
is independent of the CKM angle 7 and the input on hadronic dynamics is limited to a 
single well-controlled quantity: A S K * K *. This strategy can also be applied to extract S 
from B s — > (j)K*° and B s — ► cfxj) using the corresponding sum rule. 

In Sec. 16.51 we proposed a last method to determine <fi s , by relying on the predic- 
tion of the mixing induced CP-asymmetry A l ° n ^(B s -> K*°K*°) as a function of the 
BR lons (B d -> K*°K*°) and the theoretical input A d K * K *. In this strategy, tree pollution 
is controlled using the hadronic information from flavour symmetry and QCD factorisa- 
tion. The outcome of our analysis is presented in Fig. 16.61 This strategy requires data 
on B d -> K*°K*° and on the mixing-induced CP-asymmetry A^(B S -> K*°K*°). The 
input from B s decay is therefore minimal : A l ™^(B s — > K*°K*°), while all other inputs 
can be obtained from £>-factories. 

A comparison among the three different strategies discussed in this chapter is given in 
Table 16.41 where the needed inputs are enumerated as well as the predicted observables 
and the range of validity. 

If both hadronic machines and super- B factories |141j running at T(5S I ) provide enough 
information on i? s -decays, it will be interesting to compare the determination from <f> s 
following those methods, which rely on penguin- mediated decays, with the value obtained 
from tree processes like B s — > DK. Differences between the values obtained through these 
two procedures would provide a clear hint of physics beyond the Standard Model. In 
such a situation, the different methods presented in this letter would yield very useful 
cross-checks for the penguin-dominated vector modes. 



Chapter 7 

Supersymmetric Contributions to 



In Section B~5l we saw that the NP parameters for the decay B° — ► K + K~ can be extracted 
from the measured values for the branching ratio and the CP-asymmetries, once the SM 
hadronic parameters are known. These hadronic parameters can be obtained from other 
decays related by flavor symmetries to B° — > K + K~ , such as B^ — > -k + -k~ or B^ — > K°K°, 
provided that these are not affected by NP. In Chapter |4] we obtained the B° — > K + K~ 
SM hadronic parameters from — > 7r + 7r~ using flavor symmetry, and in Chapter [5] we 
obtained the SM hadronic parameters of B° s -» isT+if- and £° -» i^if from B° -»■ K°K° 
using an interesting combination of flavor symmetry and QCDF. 

Experimental values for the CP-asymmetries in B° — ► K + K~ are still not available, so 
the method of Chapter U] cannot be applied yet to measure the NP parameters *A U and *P U . 
However, it is always interesting to study which existing NP models can generate sizeable 
contributions to these parameters, since once A u and <& u are measured, these NP models 
will be the most sensitive to this new data. If they generate too big NP parameters, their 
parameter space will be drastically constrained. If the NP parameters are indeed measured 
to be large, then these models will prevail over more conservative ones. 

On the other hand, being able to introduce hadronic effects in a fairly reliable way 
allows to give full predictions for observables in the presence of NP. Of course, no extension 
of the SM will be given full support until the SM predictions (and in particular those 
predictions derived in Chapters [Hand [5]) are found inconsistent with experiment. However, 
due to the place of honor that supersymmetry holds within the NP models in the present, 
it is worth to study in detail how can it modify the predictions for these observables. For 
this reason, in this chapter we compute the B° — ► K + K~ and B° — ► K°K° observables in 
a generic SUSY model. 

Naively, one would guess that all NP contributions to A u and <& u are suppressed by 
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My^/M^ p , where M NP ~ 1 TeV, and are therefore small. However, SUSY contributions 
involving squark-gluino loops can be important since they involve the strong coupling 
constant a s . This contributions are thus proportional to a s /M^ p and can compete with 
the SM contributions which are of order a/M^, since (a s /a)(Myy/M^ p ) ~ 1. Indeed, 
these are the dominant effects, and are the only ones which are considered below. As we 
will see, they can account for an important enhancement of the CP asymmetries. 
The NP amplitude for a B — > f decay is given generically by (see (|4,16p ) 



The procedure to follow is the one motivated in Chapter [TJ 

• We consider all operators generated at the heavy scale, taken to be M\y. We compute 
the SUSY contributions to the coefficients of these operators. 

• Using the renormalization group, we run the operator coefficients down to the fac- 
torization scale [i. Operator mixing takes place in this step. 

• We compute the matrix elements of the various operators at uib- This allows us to 
calculate A u , ® u , A d and 

Once we have computed the values of the NP parameters, we use the SM hadronic param- 
eters calculated in Chapter [5] together with equations (|4.19p - ([4.24p to give full predictions 
for the observables. 

7.1 The Effective Hamiltonian in the presence of NP 

We begin by listing all the operators which are generated by the New Physics at the 
heavy scale. In general, NP can induce many new operators with Dirac structures that 
cannot fit in the SM operator basis presented in Section PL.2i Therefore it is mandatory to 
generalize the set of effective operators beyond the one generated by the SM. Computing 
the contributions in which we are interested, we will come up with a set of operators whose 
basis we choose as [57] : 




(7.1) 




(b a S a ) V -A{qf3q(3)v+A 
{b a S a )v-A{q~f3qf3)v-A 
(b a q a ) V -A{q0Si3)v+A 




(b a Sf3)v-A(q/3qa)v+A 
(baSp)v-A(q~l3qa)v-A 
(j>aq/3)v-A(qpS a ) V+ A 



(7.2) 
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and the corresponding effective Hamiltonian considered here is the following 



-l/NP 
^cff 



G F 



E 

i,q=u,d 



4^) of + 4^) of 



+ C &g (n) Q 8 g + Csg(fi) Qs 



(7.3) 



Despite the fact that, at the quark level, B° — > K + K~ is b — > suu, (i-quark operators 
must be included above since they mix with the u-quark operators upon renormalization 
to mi,. The same is true for B° — > K°K° and the it-quark operators. The operators Of 
and Qs g are obtained from Of and Qs g by chirality flipping, i.e, V-A <-> V+A. 

The above list of operators includes new-physics contributions to electroweak-penguin 
operators. As we will see, these effects can be significant. This shows that, although the 
SM electroweak-penguin contributions to B° s — > K + K~ are negligible, the same does not 
hold for the NP. 



7.2 Wilson Coefficients 



In order to calculate the relevant SUSY contributions to the Wilson coefficients in (17. 3D 
we must compute the corresponding diagrams in the full theory. These, as stated above, 
are the contributions involving squark-gluino loops. Specifically, we have to deal with box 
and penguin diagrams (see Figs. 17.11 and 17. 2p . The contributions from the box diagrams 
to the Wilson coefficients at the matching scale are [57] 





a 2 


c l,box — 


2V2G F m 2 - 


r d 


a 2 


c l,box — 


2V2G F m 2 - 


r u 


a 2 


c 2,box — 


2V2G F m 2 g 


r d 


a 2 


c 2,box — 


2V2G F m 2 - 


n u — 


a 2 


c 3,box — 


2V2G F m 2 - 


r d 


*l 


c 3,box — 


2V2G F m 2 - 



2V2G F m 2 ^ T ib L y T is L ( T 'juT T 'ju [ 18 F ( X di^ X ^j9) lg G ( X dig' Xa 



r l 



2V2G F m 2 - ^ L ^* r ^ 1 '^uT^u '[Q F ( x dip x ufg) + 6 G ( x dig> x ^9j 



2V2G F m 2 ^ i ^* r S I '^ r S ii )* r S ii [6 i?1 ^ d ^' :C ^' ) + Q G( - X ^9 



' x djg) 



d u u r 5 

0^ib L ) ^is L (Tju) Fju ~ g F ( X dig^ X Ujg) + 7^G(x^ i -,Xu j g / 



1 

36 



+ (r^)*r?/(rf/)^ 
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Figure 7.1: Gluino box contributions to b — > sqq. 



°4,box 



L 4,box 



r d 

u 5,box 



c 6,box 



'5, box 



digi X v-jg. 



{cr 



1 



ib L T r is L ( r jd L )* r jd [^^dig^djg) + ^G^dig^djg) 



+ (r^)*rf/(rf/)*rf/ 



F ( x d^ x d rg ) + ^ G ( x d^ x d rg )}} 



a 



2V2G F m 2 - 



1 5 

Fix t - ,it -) G(it - 



ft 



IY'Dl \*-pDR f-pDR\*-pDL 
\ L ib ) 1 js V 1 id J 1 id 



L18 



id' 



2^ GFm 2( r ib L )* r ?( r S fl )* r S L _6 F ( X dig> X d j9 ) + ^^dig^djg) 



-6,box 



9 





(7.4) 



The coefficients Cg box and Cg box are zero since their structure arises from the diagrams 
(a) and (b), but these diagrams don't contribute to u-type coefficients, as mentioned above. 
The coefficients c? box are just like (|7.4p with the replacement L «-> R. Also Xig = mf/m|, 
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Figure 7.2: Gluino-squark penguin contributions to b — > sgg. 



and the loop functions are given by 

x In x y In y 



F(x,y) 



(x - y)(x - l) 2 (y - x)(y - l) 2 (x - l)(y - 1) ' 



. x 2 lnx W 2 lnu 1 ,„ „. 

y) = I I (7 5) 

(x-y)(x-l) 2 (y-x)(y-l) 2 (x - l)(y - 1) . 

The contributions to the Wilson coefficients from the penguin diagrams are given by 

~ 9 r q „2 



„<1 _ „<? _ C 2,peng _ C 4,peng _ « s /t^l \*-p g l ^ " 

■S 



C l,peng - C 3,pcng ~ 3 " ~ '3 ~ 2 y/2G F m? 1.2^*^ + Q^^^is), 



(7.6) 

where the loop functions in this case are 

1 (1 + 2x)lnx 

{X) ~ 2(1 -x) + 6(1 -x) 2 ' 
ll-7x + 2x 2 In a; 

B{X) = ~ 18(1 -x)3 ~ 3(l-x)4. (? - 7) 

Again, the coefficients of the operators with opposite chirality are obtained via the 
substitution L <-> ii. 

It is possible to cast this Wilson coefficients into a more convenient form. To that end 
we will deal with the two 6x6 squark mass matrices M u and M d , 

(Mil M&\ 
M q2 = (7.8) 

\M& M&J 

These matrices are hermitian and so are Mil and M RW while M LR = M RL- These 

are the 

mass matrices in the interaction basis for the squark fields, {uf,dj} and are diagonalized 
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by the rotation matrices T D and T u that link this basis to the mass eigenbasis {ui,dj}, 



d 1 = {T D )U 



(7.9) 



and similar for up- type squarks. 

It turns out that the entries of the left-right down-squark mass submatrix are sup- 
pressed (at least) by a factor mb/M SVSY [57], and we neglect it. Generically, SUSY can 
induce large FCNC's. Experimental bounds - from neutral meson mixings, b — > 57, etc- 
provide strong constraints on the SUSY parameter space. In particular, Bd — Bb and K—K 
mixing allow us to decouple the down squark from the strange and bottom squarks. The 
same can be said for the up squark based on D-D mixing, although this doesn't directly 
decouple the up from the top squarks |57| . 

These simplifications imply, first, that M*^ — 0. Second, the d-type rotation matrices 
can be parametrized by two mixing angles Or^l and two weak phases Sl,r, 



( 1 









cos 9l 







sin9 L e- i5L 



sin f^e 
cos 6l 



\ 



\ 





cos 9r 
sin 



i5 R 





sin 9 R e iS 
cos 8r 



(7.10) 



The same approximation is made for T u , although we just shall need the first and fourth 
rows and no angles arise from the u-type sector in our computations. We choose the 
mixing angles to take values \6l,r\ < 45°, which defines the mass eigenstates s and b as 
the ones more closely aligned to the s and b quarks respectively. 

By substituting the simplified rotation matrices into (j7.4j) and (17. 6h we find 

a 2 sm29 L e l5L \ 1 , 5 „. . I., . 2„, 



4V2G F m 2 § 

~( X b L g ~~ ^ x S L g 

a 2 s sm26 L e i5L 
4V2G F m 2 ~ 



lo 



3 2 

) - 2 A i X b L g) ~ ^B( X - bL g 



-(x 



J s L gj 



a 2 s sm26 L e tdL 

4V2G F m 2 g 
-(x- b 



-F(x 



b^gi x qls 



+ -^ G ( X brJi X 



L91 QL9 . 



+ g B ( X h9 



b L g 



Xz 



LB) > 



r 1 



a 2 s sm26 L e ldL 
AV2G F m 2 g 

~( X f>L9 ~ 



r^{ X b L g1 X q L g) + Y^{ X b L g1 X 9Ls) 



; A ( x ~b 



(7.11) 
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with the loop functions F,G,A,B given by (|7.5p and (|7.7p . The expressions for the coef- 
ficients cf are obtained from those in (|7.1ip via the exchange L <-> R, and the coefficients 
c| 6 and c| 6 vanish. In addition Xq ig - = rri~./rn~, where (q = d,u) is the mass of the 
j-th squark mass eigenstate. 

Using the same approximations, the SUSY contribution to the WC of the chromomag- 
netic operator is given by 

x 2o VSUSY _ 8 «|sin2gxe^ [--SUSY/ \ {'l \~| r719 x 

A t ^C 89 -- i/s " (fc ~ (7.12) 

where the loop function is 

.SUSY/ x -H + 51x-21x 2 -19x 3 + 6x (-l + 9s) log(x) 
/§ {X) ~ 72(-l + x) 4 (7 - 13) 



7.3 RG evolution of the Wilson Coefficients 

Once the Wilson coefficients have been calculated at the perturbative scale Miy, the next 
step is to compute their renormalization-group running, including operator mixing, down 
to the factorization scale. This scale, as mentioned in Section is taken to be m&. The 
QCD evolution of the Wilson coefficients is given by (|1.18p . 

C(/x) = U 5 (fi,M w )C(M w ) (7.14) 

where U^{pL,M\y) is the evolution matrix, given at leading order by (|1.17j) . For our pur- 
poses it is sufficient to work at leading order in qJ3 and neglecting electromagnetic cor- 
rections (zeroth order in a). 

The easiest way to perform the running is to translate our basis (|7.2p to the usual basis 
of 12 AB = 1 SM operators in eq. (|1.28p . The coefficients Cj in (jl.27p are related to the 
NP c?'s through 



r u - 
c l — 


-At (C 5 + C 7 ) 


r d _ 
c l — 


-At (C 5 - 


fC 7 ) 


c 2 — 


-A t (C 6 + C 8 ) 


„rf _ 
c 2 — 


-A* (C 6 - 


2^8) 


r u - 
C 3 — 


-X t (C 3 + Cg) 


c 3 — 


-At (C 3 - 


2C9) 


c 4 — 


— Aj (C4 + C10) 


„rf — 
c 4 — 


-A 4 (C 4 - 


2^10 



(7.15) 



Note that the c^'g are zero at the Mw scale in our case. We take them to be zero also at 
the mb scale since the electroweak combination c§ 6 (mb) — cf q^i,) is at LO a function only 



1 The inclusion of next-to-leading-order corrections in a s in the anomalous dimension matrix has an 
impact of less than 10%. 
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of Cg'g(MvK), and the QCD combination (cf 6 (m;,) + 2c5 6 (m{>)) is mostly dominated by the 
same combination at Myy, taking into account that all NP penguin Wilson coefficients are 
of similar size. 

The running is then performed over the operators C% by means of eq. (|1.17p and the 
anomalous dimension matrix in Table [TTl 

7.4 Hadronic matrix elements for B° — ► K + K & B° — > K°K° 



The final step in the program is to compute the hadronic matrix elements of the oper- 
ators in Eq. ([73]) for B° -> K + K~ and B° s -> K°K°. These are calculated using the 
naive factorization approach as explained in Section 12.11 plus the contribution from the 
chromomagnetic operator. 

We consider first the contribution from Of, 0%, Of and O4 to B° — > K + K~ . From 
eq. (|2.2ip we have 



A(B° -»■ K + K-) 1)2 ,3,4 = X p 3) (a4 + r«a 6 + a 10 + rKa 8 )A KK 



p=u,c 



-A (s) A 



C 4 + 



c 3 



+ [c w + 



c 9 



+ r5 C 8 + 



A, 



(7.16) 



where we have used the unitarity relation (|4.ip and the relationship between the two 
different sets of Wilson coefficients (|7.15p . In the same way, for B° s — ► K°K°, one finds 



A(B S ° -> K°K°) 1>2i3A = A KK 



c2 + ^ +rC cS + -i 



(7.17) 



The contributions from the operators 0|, Or,; ^3 and 0| are the same with the substi- 
tutions c\ — ► — cf. Finally, the factorized matrix element ^4^^ (c.f. eq. (|2.32p ) and the 
chiral factor (c.f. eq. (I2.34p ) are given by 



&m%f K Fr K 



2m 2 K 



(7.18) 



The contribution from the chromomagnetic operator can be obtained from the order 
a s QCDF formulae as found, for instance, in [26]. In reference to the QCDF coefficients 
in eq. (|2.46p . the coefficient CfJ contributes only to the penguin contractions P%(K) and 
Pq(K). These only enter in a\{KK) and cLq(KK), so at the end the contribution to the 
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amplitude is 



A(B° 8 - KK)t 



CfOs 

4vriV, 



-2C; 



eff 



1 



X 



■$ K (x) - 2Cf g r« 



(7.19) 



The integral over the kaon light-cone distribution amplitude can be expanded in terms 
of Gegenbauer polynomials as in eq. (|2.36p . which gives integral = 3 + 2>af + Salf. 
However, because we are working in naive factorization, we may just take the asymptotic 
form for the distribution amplitudes and set af — ► 0. 

Then, the NP amplitudes for B° s -> K + K~ and B° s 
read 



K K in naive factorization 



A u e 



(K + K-\H»£\B°) 



A, 



c 4 



(K°K°\H^\B°) 



1 + 



A, 



-d 



c 4 + 



+ < 



4 + 



1 + 



(7.20) 



where the barred coefficients are q = c, — q and CojJ = Coj? — (71^. 



7.5 NP contributions to 5° — B° mixing 



As mentioned in Section I4.5|, the NP will in general affect also the mixing, and not only 
the decay amplitude. As the NP correction to this angle - 5(j)^ p - is necessary in order 
to compute the mixing induced CP asymmetry (|4.2ip . we must consider the SUSY con- 
tributions to B° — Bg mixing, M]^ p . Adding this new contribution, the mixing angle 4> s is 
given by (pT/Tj) 



arg(M 1 s 2 M +Mf 2 p ) 



arg 



Mtf 



1 + 



'12 



ji rsn 



arg(M 1 s 2 M )+arg 1 + 



(7.21) 



so that the NP contribution to the mixing angle is 



Mf 2 p 
( 1 + T7k 



M 



12 



(7.22) 
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These contributions, in our case, come from the squark-gluino box diagrams in Figure 
17.31 These diagrams are calculated like the ones in Figure [7TT1 (see for instance |142| ) and 
match to the following AB = 2 effective Hamiltonian, 



(7.23) 



i=l i=l 

where the operators are defined as 

0\ = {sLj^,b L )(s L ^b L ) , 4 = (s R b L )(s L b R ) 

2 = (s R b L )(s R b L ) , Oj = 

3 = (4^)(4^) > o t = Oi\ L ^ R . 



(7.24) 



When only LL or RR mixing is present in our scenario, only C\ or C\ are generated. When 
both LL and RR mixings coexist, the operators O4 and O5 contribute. However, in our 
scenario the operators O2 and O3 are not generated at all. The Wilson coefficients at the 
matching scale My/ are given by 
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(7.25) 
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Figure 7.3: Gluino box contributions to B° — B° s mixing. 



where the loop functions F and G are given in (|7.5j) . The Wilson coefficients at the scale 
m b are obtained applying the RG procedure. The two loop running gives |143j 



Ci(m 6 ) = 0.848 Ci(A%) , C 4 {m b ) = 2.395C 4 (M W ) + 0A85C 5 {M W ) , 
Ci(m b ) = 0M8C 1 (M W ) , C 5 (m b ) = 0.061 C A {M W ) + 0.904 C 5 (M W ) . 

The calculation of the matrix elements in terms of the B-parameters gives 

2 



(7.26) 
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M Bs 



1 ~ 1 

-Bi{m b ){Ci{m b ) + Ci{m b )) + -B 4 (m b ) 
3 4 \ fJ-b + m t 



C 



4 (m 6 j 



4 B5 <™'» Ufa 



C 5 (m b ) 



where the B-parameters can be found in [144| and their central values are 

Bi{m b ) =0.86, B A {m b ) = 1.17, B 5 (m b ) = 1.94. 
In the SM, only Cf M is non-zero, and it is given by 



cr(Mw) 



G F M W x (s) 2 q ( , 



where x t = mf/M^r and Sq is [8] 



Sn(x) 



c 4 — 12x 3 + 15x 2 — 4x + 6x 3 In 
4(x - l) 3 



X 



(7.27) 



(7.28) 



(7.29) 



(7.30) 



128 Su^er symmetric Contributions to B° — > K + K and B° — > K°K° 



Finally we can write 

Cf M (m b ) + 4^1(7^6) \m b + m s J Cf M (m b ) 

l B 5 (m b ) ( M Bs \ 2 CI^(m b ) \ 

4 5iN \m b + mj C?"(m b ) J 1 " ' 

We have now all the ingredients to calculate the NP amplitude and the B° — ► ETA' ob- 
servables in terms of the set of SUSY parameters. 



arg 1 + 



+ 



7.6 Additional Bounds on the SUSY parameter space 

Let us address now the issue of the input parameter space. To that end we begin summa- 
rizing the basic points discussed so far. 

• The final aim is to see how the space of allowed values for the B® — ► K + K~ and 
Bg — > K°K° observables within a general SUSY model is increased with respect to 
that of the SM alone. 

• The expressions for full observables in the presence of NP are given in (|4.19p - (|4.24j) 
and (14.8p - (|4.10p . These depend on the following set of parameters: |T XJr |, d KK , 9 KK , 
7, Al s) , <t> s s M , 5^ p , A u and * u . 

• The CKM parameters 7, A{f and 0f M = 2/3 s are measured independently. 

• The SM hadronic parameters |T XK |, d KK and 6 KK are obtained from B® — ► K°K° 
as described in Chapter [5l paying especial attention to the correlations, since these 
parameters must be varied together with the scanning of the NP parameter space. 

• The NP parameters A u , A d and $^ are the magnitudes and weak phases of 
the NP amplitudes. This amplitudes have been calculated in this chapter within 
supersymmetry, as a function of the following SUSY parameters: 

— The masses of the squarks and the gluino. We take = nig = mg = 
m~ H = 250 GeV and 250 GeV < mu R ,m SLR < 1000 GeV. 

— The mixing and weak angles. These are taken unconstrained: — 7r/4 < 9l.r < 
7r/4 (see Section [7^2]) and —it < 5l,r < tt. 

The Wilson coefficients are sensitive to the s — b mass splitting. They vanish for nig = 
m~ b and grow when the splitting is large. We therefore expect these contributions to 
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be most important for large values of mj (keeping m~ h fixed). In the same way, NP 
effects in b — > dqq transitions depend on the difference m^ — m^. By setting mg = 
we ensure that b — > d decays get no such contributions, which is consistent with the 
discussion in Sec. 17.21 A difference between A u e® u and A d e^ d is only possible in the 
presence of a nonzero u — d mass splitting. Without it there are no contributions 
to isospin-violating operators. However, this mass splitting must be very small in 
the left-handed sector due to SU(2)l invariance. We therefore set rriu L = m ( i L , but 
allow for a significant mass splitting in the right-handed sector. 

A quick look at the form of the Wilson coefficients in (|7.f 1|) shows that the NP ampli- 
tude (|7.20p has the following structure, 

A u e iS> u = (jFsin20 L e^ - ^sin 26 r^) A (7.32) 

where J- and J- are functions of the masses of the squarks and gluino. They are related 
by a change L «-> R, and for rriq L = rriq R they are equal (J- = T). Clearly, a way of 
maximizing the NP amplitude is to set 6r = Ql = ir/4 (maximal s — b mixing) . In such a 
situation we have, 

A H = 2F\A\\8m(8 L -$ u )\ (7.33) 

which shows that if 8l r are unconstrained, then A u and & u are not correlated and can 
be maximized independently. 

However, the range in the SUSY parameter space has to be constrained to fit to several 
existing experimental bounds. In the following we show what these bounds are and how 
they constrain the various ranges taken for the input parameters. 



7.6.1 Bounds from B — > ixK 

The same SUSY contributions to B° — > KK will also affect B — > irK, since they both 
share the same quark-level decay of the b quark. In particular, there will be effects on the 
quantities R* and Acp(tt + K°) [26], whose definitions and measured values are [85] 

R - BR{B + ^v + K»)+BR(B- 

R * = 2[BR(B+ *°K+) + BR(B~ -> ^K~)] ~ L °° ± °-° 8 (7 ' 34) 

, , , n , BR(B + ir + K°) - BR(B~ -> ir-R ) 

Acp(k + K ) = i— ■ ■ — ^ =r-£ = -0.020 ± 0.034 7.35 

y ' BR(B + — > 7T + K°) + BR(B~ — > tt~K°) y ' 



In order to incorporate these two constraints we follow the approach in [57] where QCDF 
is used, except for the strong phase related to ^cp which we take as a free parameter. 
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7.6.2 Bounds from B — > X s 7 

The bounds from B — > X s j are usually important constraints on any FCNC NP analysis, 
and this is no exception here. For the branching ratio of B — ► X s ^f we take [85] 



BR(B -» Y s7 ) = (3.55 ± 0.26) x 10" 4 (7.36) 

and we consider a 2-cr range, allowing for various theoretical uncertainties. As we shall 
see, even a wide range of this sort has a considerable impact in our analysis. This bounds 
are incorporated following [145] . 

7.6.3 Bounds from B° — B° mixing 

Data from B° — B° mixing also provides a potential constraint in the SUSY parameter 
space. In particular, the recently measured mass difference AM S in eq. (|3.6ip [86], together 
with the SM fit [731 [74] gives 

AM S \ . . 

(0.81 ±0.19) , (7.37) 



AM SM , 

s / exp 

which constrains the SUSY parameter space through the modulus in (|7.3ip . since from 
(I3.76|) we have 



AM S _ 2\Mfg + M$\ 



Ai\/| M 2|Afg 4 | 



1 + 



(7.38) 



This bound also turns out to be important. 



There are other traditionally important constraints, like B° s — > that are very 

sensitive to other SUSY parameters, mostly tan/3 and tua- However, for low tan/3 and 
values of above 200 GeV, they have no effect on our region in the SUSY parameter 
space. 

Taking into account the various constraints, the contributions from LL and RR mixing 
have been analyzed. AM S is the strongest constraint, and it is the relevant one when 
considering only LL or RR mixing separately. In particular, it has a large impact on the 
phases $ u and <&d- I n the case of B° — ► K + K~, LL mixing gives the largest contribution 
to the amplitude, more than twice that of RR. However, in the case of B° s — > K°K° both 
contributions are similar in size. 

When both LL and RR mixings are allowed simultaneously, the constraints on the 
SUSY parameter space are changed. In this case new operators for B° — B° mixing are 
generated, so that the effect is not simply the combination of the two separate contributions 
(for instance, see Ref. [146] ). We find that (i) now BR{B — > X s j) is also important, not 
only AM S , and (ii) the global effect of the constraints is weaker. The upshot is that 
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Figure 7.4: SUSY contribution to the NP amplitudes A v e i<s> » (left) and A d e^ d (right) in the 
scenario with simultaneous LL and RR mixings. The dark regions correspond to the variation of the 
SUSY parameters over the considered parameter space. The light regions satisfy the experimental 
bounds, including the recent measurement of AM,. 



there is a certain enhancement of the NP amplitudes when both LL and RR mixings are 
combined. In particular, the weak phases <I> U and <J>d are not so strongly constrained as 
when either LL or RR mixing is taken to vanish. 

7.7 Results: NP amplitudes and Observables 

In Fig. 17.41 we show the allowed ranges for A u e l ^ u and A d e l ^ d in the scenario with si- 
multaneous LL and RR mixings. The dark regions correspond to the values that these 
amplitudes take when varying the parameter space over the initial ranges. The light 
regions show how these values are reduced by the existing experimental constraints men- 
tioned above. There are two important remarks. First, we see that the above constraints 
do indeed greatly reduce the allowed SUSY parameter space. Second, even so, the effect 
on A u e l ® u and A d e l ® d can be significant. 

At this stage we can identify what are the effects of the various constraints in reducing 
the SUSY parameter space. The bound from B — > X s ^/ affects only the left-handed 
sector. In particular, for large mg, the regions with \0l\ > 10°, \5l\ < 60° and \9l\ > 10°, 
\5l — tt| < 60° are excluded. The bound from AM S is much stronger: when rag > 400 GeV, 
any values of 16*^^1 ^ 5° are excluded, as well as those regions in which 5l + 5r ~ 



132 Supersymmetric Contributions to B° — > K + K and B° — > K°K' 



-3tt/2, -vr/2, tt/2, 3tt/2J£J After these bounds are imposed on the parameter space, the 
constraints from B — > ttK have very little effect on the regions in Fig. 17.41 

Note that the allowed region for A u is much larger than that for A d , by approximately 
a factor of 3. In the isospin limit, these should be equal, so this factor of 3 is a measure of 
isospin breaking in this NP scenario. In particular, for mn R = 250 GeV (zero ur-cLr mass 
splitting), the values of A u reduce to those for A d . 

We now examine the effect of these contributions on the observables. By adding the 
SUSY contributions to the SM amplitudes, it is possible to compute the branching ratio 
and the CP asymmetries in the presence of SUSY. Fig. l7.5l shows the allowed values for the 
B° — ► K + K~ observables, for three different values of A^, compared with the predictions 
of the SM and with the recent experimental value for the B° — > K + K~ branching ratio 
reported by the CDF collaboration [130j : 

BR(B° S -> K + K~) CKp = (24.4 ± 1.4 ± 4.6) x 1(T 6 . (7.39) 

The agreement between the CDF measurement and the prediction of the SM in Ref. [35] 
erases any discrepancy between experiment and the SM. This branching ratio will now be 
an important future constraint. The branching ratio within SUSY should not deviate much 
from the SM prediction so as not to generate any disagreement with data. Indeed, Fig. 17.51 
shows that the impact of SUSY on the branching ratio of B° — > K + K~ is practically 
negligible. Interestingly, for positive values of A^ r (preferred region), the SM predicts a 
smaller value for BR(B° — > K + K~), but it is now compatible with the new data. Still, it 
is in this case that SUSY shows a larger deviation in the correct direction. 

A completely different picture arises for the CP asymmetries. The results for the 
direct CP asymmetry reveal that SUSY can have an impact. This is not surprising: 
SUSY introduces a term in the total amplitude which is of the same order of magnitude 
as that of the SM and carries a weak phase that is not constrained. The mixing-induced 
CP asymmetry gets affected in a more dramatic way. The interpretation is that the SUSY 
contribution to the mixing angle (j) s can be large (in fact it can take all values between — ir 
and 7r), while in the SM it is tiny: 0f A/ — —2°. Any experimental measurement falling 
inside the dark area in the plots, but outside the dashed rectangle, would not only signal 
NP but clearly could be accommodated by supersymmetry. 

Fig. 17.61 shows the results for B° — ► K°K°. Although the branching ratio is little 
changed in the presence of SUSY, the enhancement of the CP asymmetries due to the 
inclusion of the SUSY contributions is in this case even more important. The reason is 
that, within the SM, the CP asymmetries are much smaller in B° — > K°K° than they are 
for B° — ► K + K~ , because of the absence of the tree diagram. Thus the impact of SUSY 
is much greater. This is evident by looking at the three lower plots in Fig. 17.61 where 



For further details of the AM a constraint on this parameter space, see Ref. [146] . 
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Figure 7.5: Predictions, in the form of scatter plots, for the correlations between BR(B° — > 
K+K-) - A dh (B° -► A'+A"-) (up) and A mix (.B° -> K + K~) - A dh (B° -> K+K~) (down) in 
the presence of SUSY, for (a) Af lr = -0.1, (b) ^9. = and (c) Af ir = 0.1. The dashed 
rectangles correspond to the SM predictions. The horizontal band shows the experimental value 
for BR(B° S -► K+K-) at la. 




(a) (b) fc) 



Figure 7.6: Predictions, in the form of scatter plots, for the correlations between BR(B° — > 
K°K°) - A dil (B° s -> K°K°) (up) and A mix (B° s -> A°A°) - A dir (B° -> A A ) (down) in the 
presence of SUSY, for (a) Af ir = -0.1, (b) = and (c) Af lr =0.1. The dashed rectangles 
correspond to the SM predictions. These are quite small in the three lower plots, so they are 
indicated by a circle. 



134 Supersymmetric Contributions to B° — > K + K and B° — > K°K' 



the tiny rectangles corresponding to the SM predictions can hardly be observed. We have 
drawn a circle around them to indicate their position. 

These are a good illustration of the general scenario discussed in the introduction. 
While the branching ratios in this case are relatively insensitive to supersymmetry, the 
direct and mixing-induced CP asymmetries of these decays are greatly affected. Thus, 
these CP asymmetries are the observables to focus on in order to observe NP, particularly 
SUSY, while the branching ratio of B° — ► K + K~ can become an important constraint on 
models beyond the SM other than SUSY. 

7.8 Summary of the results 

In this chapter we have considered the branching ratios and CP asymmetries for the 
decays B° s — ► K + K~ and B° — ► K°K° in a supersymmetric (SUSY) model, focusing on 
the dominant gluino-squark contributions [57] . The determination of the SM contributions 
has been taken from the combination of B% — ► K°K° and QCD factorization described in 
Chapter [5] [35]. In ref. |128j the SM contributions were taken from — ► tt + it~, as was 
done in Chapter HI [125] . but the uncertainties are much bigger. Here we have described the 
analysis of ref. [36] • We have also included the constraints coming from BR(B — > X s -y), 
B — > ttK and AM S , and we find the following results. 

• The new-physics (NP) amplitudes are A u e^ (B° s -> K + K~) and A d e^ d (B° -> 
K°K°). We find that both can get significant contributions from SUSY. In the 
isospin limit, these quantities are equal. However, our calculations show that, for 
the region of parameters considered, in SUSY there can be a difference of up to a 
factor of 3 between the NP amplitudes. This indicates the possible level of isospin 
breaking in this type of theory. In particular, in the SUSY model considered here, 
large isospin violation is possible when there is large mass splitting in ur-oIr. 

• The branching ratio BR(B° — ► K + K~) is very little affected by SUSY. At most, the 
SM prediction can be increased by 15% for Af ir = 0.1. In fact, SUSY can somewhat 
improve the already good agreement between the SM prediction and the new precise 
CDF measurement pJ30]. The impact of SUSY on BR(B° -»■ K°K°) is even smaller, 
reflecting the reduced allowed region for A d e l ^ d as compared to A u e l ® u . 

• The situation is very different for the CP asymmetries; the size of the effect depends 
strongly on the decay and the type of asymmetry. For B° — ► K + K~ , the direct CP 
asymmetry within SUSY is in the range -0.1 < A dir (B° s -> K + K~) SUSY < 0.7 for 
-0.1 < Af ir < 0.1. Depending on the value of Af ir , it may be possible to disentangle 
the SUSY contribution from that of the SM. This is due to the competition between 
the tree and the NP amplitudes for each value of A^ r . As for A m i x (B° — > K + K~), 
its value can vary all the way from —1 to +1, signaling a large impact from SUSY. 
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• Turning to B° s — > K K , the CP asymmetries are particularly promising. This 
decay is dominated by the penguin amplitude in the SM, and so the direct CP 
asymmetry is strongly suppressed: it is predicted to be at most of the order of 1%. 
However, in the presence of SUSY, the direct CP asymmetry can be 10 times larger. 
The mixing-induced CP asymmetry is also predicted to be very small in the SM. 
However, A m i x (Bg — > K°K°) SUSY covers the entire range, and so this asymmetry 
can be large in the presence of SUSY. 
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Part III 

Conclusions 



Chapter 8 

Concluding Remarks 



Non-leptonic weak decays of B mesons are suitable for the combined study of different 
aspects in the phenomenology of particle physics. We have seen their relevance in studies 
of CP violation, strong interactions, heavy quark physics, flavor physics and physics beyond 
the SM. All these topics are interconnected and related to the search for new physics; a 
search that will most probably succeed in the next few years, but with consequences that 
will surely require theoretical and experimental work for a much longer time. 

This is a very interesting moment in the history of particle physics, in which we expect 
to finally break the TeV wall and go beyond a theory of elementary particles formulated 
forty years ago. The outcome can be a completely new understanding of the world we 
live in, leaving aside the period of speculation and uncertainty of the last decades. A 
joint effort of the whole particle physics community has brought us here, and the same 
effort will be necessary to make this breakthrough possible. It is an impressive world-wide 
project that will have a historic result. 

The direct search for new particles is a crucial part of this program. The main objective 
of the LHC, with its experiments ATLAS and CMS, is to find the elusive higgs particle, 
the last requirement of the SM. However, a more important task will be the search for 
other completely new particles that might help to understand all the puzzles that the SM 
leaves unanswered. In fact, from the theoretical point of view, the observation of new 
particles might be much more valuable than the establishment of the higgs' existence. 

Despite the importance of the direct searches, the correct identification of the new 
physics and its properties, and the understanding of its theoretical consequences, requires 
also from indirect analyses. Some of the major open questions in particle physics that 
might find a partial answer in the next future experiments are related to flavor physics. 
B physics has been the center of attention in the experiments of the last decade, with a 
huge amount of data collected at the B factories and at Tevatron, and will keep a place of 
honor in the future, with the B physics experiment LHCb at the LHC and a future Super 
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B. The reason for this is that B physics is a powerful probe of physics beyond the SM in a 
completely complementary way to direct searches. In fact, all the B physics experiments 
carried out so far have put stringent bounds on the nature of the new physics, so thus far 
the indirect searches are on the head of the list. And whether or not the direct production 
of new particles takes place, B physics will still be an important way to test the SM and 
look for new physics. 

However, in order to look for tiny new physics effects, the SM predictions must be 
under strict control. In this thesis we have focused on non-leptonic B decays. These type 
of decays are very difficult to study because of their sensitivity to long distance strong 
interaction effects. So in order to look for new physics in these decays one must find a way 
to deal with hadronic phenomena leading to predictions of acceptable accuracy. Indeed, 
the claim of a NP signal in non-leptonic B decays will have to rely on a solid inclusion of 
long distance QCD contributions. This also means that the study of these decays might 
have a double prize: understanding the NP, and understanding QCD. We have already 
obtained benefits from both sides. 

Towards this end, we have presented analyses of two body non-leptonic B decays within 
two approaches. The first one is the flavor symmetry approach, which takes advantage of 
the approximate SU(3) symmetry of the QCD lagrangian to establish relations between 
hadronic quantities for different decay modes, in such a way that predictions can be made 
without the need of direct computations of these hadronic quantities. One can then use 
data on some sort to predict the values of other observables. In Chapter|3]we have used this 
approach to extract the SM hadronic parameters in B s — > K + K~ decays from Bj — ► 7r + 7r~ 
data. This has provided predictions for B s — > K + K~ observables within the SM that can 
be compared with experiment. A deviation between the theoretical predictions and the 
experimental data will be a signal of NP. A useful general parametrization of the NP 
contributions to this decay has been given, together with a method to measure the new 
physics parameters. Therefore, if NP was found in this decay, the measurement of these NP 
parameters would provide interesting information on the nature of this NP, as for example 
if the NP is isospin violating, and so on. The limitations of the SU (3) analysis have been 
discussed in detail, in particular the problem of estimating the symmetry breaking effects, 
that constitute the main source of uncertainties. 

A second approach that has been discussed here is the QCD factorization method. 
This provides a tool to compute matrix elements from first principles, based on a heavy 
quark expansion. Therefore this approach provides genuine QCD predictions. However, 
some of the phenomenological difficulties of this approach are those non-factorizable power 
suppressed contributions that are numerically important for realistic b quark masses. In 
Chapter [5] we have proposed a method to get around this difficulty. A combination of 
QCDF with flavor symmetry relations has been used to study B s — > KK decays, giving 
predictions for the observables with uncertainties reduced by more than a 50%. This 
method has also been used in Chapter [6] for B s — ► VV decays, and used to propose 
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three strategies to measure the B s — B s mixing angle. This mixing angle is predicted to 
be extremely small within the SM, and it is therefore an excellent place to look for NP. 
However, the claim of a NP contribution to this mixing angle cannot be easily disentangled 
from a NP contribution to the decay amplitudes in these type of decays, and we have 
discussed the ways in which these strategies could help to make this distinction. 

This proposal is clearly an improvement on the flavor symmetry and the QCDF ap- 
proaches, giving more reliable results with smaller uncertainties. In fact, the results 
presented for B s — > K + K~ were predictions that were confirmed very nicely -maybe 
unfortunately- by the CDF collaboration shortly after. We would like to make some 
further comments concerning this method. The applicability of the approach has to be 
checked individually for each mode: it can only be applied to those decays for which A 
receives no contributions from annihilation or hard spectator scattering graphs, such as 

— > tt + 4>, etc, as discussed in Section \5. 11 The predictions 
derived in this way include most of the long distance physics, which is contained inside 
the experimental input. The used theoretical input is minimal, and is the most reliable 
input that QCDF can offer, free from the troublesome IR-divergencies. Moreover, the 
theoretical error is under control and is likely to be reduced in the near future due to, for 
example, the fast progress that is taking place in lattice simulations. 

There is, however, an honest criticism due to the fact that some long distance effects 
that are controversial in QCDF, are also absent here. The prominent one is the contribu- 
tion from the charming penguins, which could easily account for a significant discrepancy 
between theory and experiment that would not be due to NP. 

Finally, in Chapter [7] we have studied B s — ► KK decays in supersymmetry. Using 
the SM analyses performed in the previous chapters, and computing the most relevant 
contributions from SUSY to the decay amplitudes and to the mixing, we have studied the 
departure from SM expectations that SUSY could provide in the branching ratios and CP 
asymmetries. The conclusion is that the enhancements are important, which means that 
moderately precise experimental measurements on these observables, if consistent with the 
SM predictions, will provide useful constraints on the SUSY parameter space, especially 
on the CP violating phases. 

To be able foresee what is the most relevant direction that future work on these topics 
should follow requires a considerable scientific maturity, which I lack. However, there are 
some general aspects concerning this issue that are of a general opinion. 

On the QCD side, an all-order proof of factorization within QCDF would be appealing, 
but maybe not too useful for practical purposes. A further development of the SCET 
formulation will probably put factorization theorems on a much more firm basis, specially if 
the issue of the long distance charm loops finds a phenomenologically satisfactory solution. 
The final aim is to have a completely systematic way to compute QCD corrections to 
matrix elements up to any desired accuracy (even if it means that the actual computation 
is extremely difficult and lengthy). Moreover, a similar progress should be expected from 
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the lattice community, which is certainly conceivable. 

On the NP side, the progress should be faster, since immediate future events will 
have a drastic influence on any work in this direction. The SM fits are almost completely 
systematized, and any tension between theory and experiment can be dragged to the origin 
with ease. Probably a similar systematization should be built for the largest possible of 
phenomenologically relevant NP models that exist in the literature, so that the explosion 
of new phenomena that might occur in the first five years at the LHC can be used to 
accept and discard models in the blink of an eye. 

However, the true model of particle physics at the TeV scale might have nothing to 
do with any of the NP models that have been proposed up to now. Therefore, model- 
independent studies of NP are necessary in order to be sure that we are not missing 
anything new, or old. 
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